TOUGH4 User Manual
  • Quick Entry to Keywords for Data Input
  • 1️⃣INTRODUCTION
    • About TOUGH
    • TOUGH Development History
    • TOUGH4 Implementation
    • Scope and Methodology
  • 2️⃣WHAT IS NEW IN TOUGH4
  • 3️⃣CODE COMPILATION AND INSTALLATION
    • Setup for Compilation
    • Code Compilation
      • 1. Compilation of TOUGH4 using Visual Studio
      • 2. Compilation of TOUGH4 on Linux-like platform
    • Installation
    • Running the Executable for Simulations
  • 4️⃣GOVERNING EQUATIONS
    • Mass-Balance Equation
    • Accumulation Terms
    • Flux Terms
    • Sink and Source Terms
    • Virtual Node Well Treatment
    • Semi-Analytical Conductive Heat Exchange
    • Drift Model
    • Non-Darcy Flow
  • 5️⃣NUMERICAL METHOD
    • Space and Time Discretization
    • Interface Weighting Schemes
    • Initial and Boundary Conditions
      • Initial Conditions and Restarting
      • Neumann Boundary Conditions
      • Dirichlet Boundary Conditions
      • Atmospheric Boundary Conditions
      • Constant Temperature Boundary Conditions
    • Parallel computing schemes
    • Linear Solvers
    • Python Functions
      • Relative Permeability
      • Capillary Pressure
      • Initial Condition Calculation
      • Fetching Output Data
      • Fetching Thermophysical Property Data From NIST Webbook
      • Coupling With Third-Party Software
  • 6️⃣SOFTWARE ARCHITECTURE
    • Program Design
    • Data Structure
    • Linear Equation Setup
  • 7️⃣PROCESS MODELING
    • EOS1
    • EOS2
    • EOS3
    • EOS4
    • EOS6
    • EOS7
    • EOS9
    • ECO2
    • EWASG
    • TMVOC
    • Tracers/Decay Chain
    • Biodegradation Reaction
    • Wellbore Flow
    • Non-Darcy Flow
    • Enhanced Coal Bed Methane
  • 8️⃣PREPARATION OF MODEL INPUT
    • Input Formatting
    • Keywords and Input Data
      • TITLE
      • BIODG
      • CBMDA
      • CHEMP
      • COFT
      • CONNE
      • COUPL
      • DIFFU
      • ELEME
      • ENDCY
      • ENDFI
      • FLAC
      • FNIST
      • FOFT
      • FORCH
      • GASES
      • GENER
      • GOFT
      • HYSTE
      • INCON
      • INDOM
      • MESHM
      • MODDE
      • MOMOP
      • MULTI
      • OUTPU
      • PARAM
      • ROCKS
      • ROFT
      • RPCAP
      • SELEC
      • SOLVR
      • SPAVA
      • TIMBC
      • TIMES
      • TRACR
      • WELLB
    • Inputs for Initial Conditions
      • EOS1
      • EOS2
      • EOS3
      • EOS4
      • EOS6
      • EOS7
      • EOS9
      • ECO2
      • EWASG
      • TMVOC
    • Geometry Data
      • General Concepts
      • MESHMaker
      • Multiple-continuum processing
    • Inputs for MESHMaker
      • Generation of radially symmetric grids
        • RADII
        • EQUID
        • LOGAR
        • LAYER
      • Generation of rectilinear grids
      • MINC processing for fractured media
    • Adjustment of Computing Parameters at Run-time
  • 9️⃣OUTPUTS
  • 🔟VALIDATION AND APPLICATION EXAMPLES
    • EOS1
      • Problem 1 - Code Demonstration
      • Problem 2 - Heat Sweep in a Vertical Fracture (rvf)
      • Problem 3 - Five-spot Geothermal Production/Injection (rfp)
      • Problem 4 - Coupled Wellbore Flow (r1q)
      • Problem 5 - Five-Spot Geothermal Production/Injection under extremely high temperature
    • EOS2
      • Problem 1 -Five-spot Geothermal Production/Injection (rfp)
    • EOS3
      • Problem 1 - Code Demonstration (eos3p1)
      • Problem 2 - 1D TH Problem with Heating and Gas Source (by Guanlong Guo)
      • Problem 3 - Heat Pipe in Cylindrical Geometry (rhp)
      • Problem 4 - 3D Thermal Consolidation Test, Coupling with FLAC3D Simulator (by Guanlong Guo)
    • EOS4
      • Problem 1 - Code Demonstration (eos4p1)
      • Problem 2 - Heat Pipe in Cylindrical Geometry (rhp)
    • EOS6
      • Problem 1-Validation with EOS2
      • Problem 2-Noble Gas Transport
    • EOS7
      • Problem 1-Multiphase and Nonisothermal Processes in a System with Variable Salinity (rf1)
      • Problem 2-Thermal and Tracer Diffusion (EOS7R/rdif7)
      • Problem 3-Contamination of an Aquifer from VOC Vapors in the Vadose Zone (EOS7R/rdica)
      • Problem 4-Density, Viscosity, Solubility, and Enthalpy of Real Gas Mixtures (EOS7C/SAM7C1)
      • Problem 5-CO2 Injection into a Depleted Gas Reservoir (EOS7C2/SAM7C2)
      • Problem 6- CO2 Injection into a Saturated System (EOS7C/SAM7C3)
      • Problem 7-Density, Viscosity, and Enthalpy of Real Gas Mixtures (EOS7CA/SAM7CA1)
      • Problem 8-CO2 Injection into a Shallow Vadose Zone (EOS7CA/SAM7CA2)
      • Problem 9-Non-Isothermal Compressed Air Energy Storage in Reservoir (by Julien Mouli-Castillo)
    • EOS9
      • Page 1
    • ECO2
      • Problem 1-Demonstration of Initialization Options (ECO2N/rtab)
      • Problem 2-Radial Flow from a CO2 Injection Well (ECO2N/rcc3)
      • Problem 3-CO2 Discharge Along a Fault Zone (ECO2N/r1dv)
      • Problem 4-CO2 Injection into a 2-D Layered Brine Formation (ECO2N/rtp7)
      • Problem 5-Upflow of CO2 along a Deep Fault Zone (ECO2M/r1d)
      • Problem 6-Migration of a CO2 Plume in a Sloping Aquifer, Intersected by a Fault (ECO2M/rwaf)
      • Problem 7-GCS/GHE with a double-porosity reservoir (Case6_50kg_DP/ECO2NV2)
    • EWASG
      • Problem 1 - Brine Density Calculation (dnh)
      • Problem 2 - Production from a Geothermal Reservoir with Hypersaline Brine and CO2 (rhbc)
    • TMVOC
      • Problem 1-Initialization of Different Phase Conditions (r7c)
      • Problem 2-1-D Buckley-Leverett Flow (rblm)
      • Problem 3-Diffusion of components (rdif2)
      • Problem 4-Steam Displacement of a NAPL in a Laboratory Column (rtcem)
      • Problem 5-Steam Displacement of a Benzene-Toluene Mixture in a Laboratory Column (rbt)
      • Problem 6 -Air Displacement of a NAPL from a Laboratory Column (rad)
      • Problem 7-NAPL Spill in the Unsaturated Zone (r2dl)
    • T4.Well
      • Problem 1-Steady-state two-phase flow upward
      • Problem 2-Non-isothermal CO2 flow through a wellbore initially full of water
  • CONCLUSION REMARKS
  • REFERENCES
  • ACKNOWLEDGEMENT
  • Appendix
    • ☑️A: RELATIVE PERMEABILITY FUNCTIONS
      • IRP=1 Linear function
      • IRP=2 Power function
      • IRP=3 Corey's curves
      • IRP=4 Grant's curve
      • IRP=5 Perfectly mobile
      • IRP=6 Fatt and Klikoff function
      • IRP=7 van Genuchten-Mualem Model
      • IRP=8 Verma function
      • IRP=10 Modified Brooks-Corey Model
      • IRP=11 Modified van Genuchten Model
      • IRP=12 Regular hysteresis
      • IRP=13 Simple hysteresis
      • IRP=31 Three phase perfectly mobile
      • IRP=32 Modified Stone's first 3-phase method
      • IRP=33 Three-phase Parker's function
      • IRP=34 Alternative Stone 3-phase
      • IRP=35 Power-law function
      • IRP=36 Faust for two-phase Buckley-Leverett problem
      • IRP=37 Another alternative to Stone function
      • IRP=40 Table lookup
      • IRP=41 User-Defined relative permeability function
    • ☑️B: CAPILLARY PRESSURE FUNCTIONS
      • ICP=1 Linear function
      • ICP=2 Function of Pickens
      • ICP=3 TRUST capillary pressure
      • ICP=4 Milly’s function
      • ICP=6 Leverett’s function
      • ICP=7 van Genuchten function
      • ICP=8 No capillary pressure
      • ICP=10 Modified Brooks-Corey Model
      • ICP=11 Modified van Genuchten Model
      • ICP=12 Regular hysteresis
      • ICP=13 Simple hysteresis
      • ICP=31 Parker et al 3-phase function
      • ICP=32 Parker 3-phase function, alternative 1
      • ICP=33 Parker 3-phase function, alternative 2
      • ICP=34 Parker 3-phase function, alternative 3
      • ICP=40 Table lookup
      • ICP=41 User-Defined capillary pressure function
    • ☑️C: ADDITIONAL PROGRAM OPTIONS
    • ☑️D: DESCRIPTION OF FRACTURED FLOW
      • Multiple Continuum Approaches
      • Active Fracture Modle
Powered by GitBook
On this page
  1. PROCESS MODELING

TMVOC

PreviousEWASGNextTracers/Decay Chain

Last updated 1 year ago

  1. Description

TMVOC was originally developed by Karsten Pruess and Alfredo Battistelli. It was designed for simulating the flow of multicomponent mixtures of volatile organic chemicals (VOCs), such as crude oil, gasoline or diesel, and organic solvents, in variably saturated media. It can be used to analyze the fate and transport of non-aqueous phase liquids (NAPLs) in the vadose zone as well as below the water table. The NAPL may consist of a general multicomponent mixture of organic fluids. In addition, one or several non-condensible gases may be present. Any and all phase compositions in a gas-water-NAPL system are treated, including single, two-phase, and three-phase conditions. Flows can be non-isothermal, and may involve advective, diffusive, phase-partitioning, and sorptive processes.

In the TMVOC formulation, the multiphase system is assumed to be composed of water, non-condensible gases (NCGs), water-soluble volatile organic chemicals (VOCs), and optional tracers. The number and nature of NCGs and VOCs can be specified by the user. There are no intrinsic limitations to the number of NCGs or VOCs in the TMVOC formulation; Current default total number of components are 20 (ComponentNum=20). NCGs may be selected by the user from a data bank provided in TMVOC; currently available choices include O2, N2, CO2, CH4, H2, ethane, ethylene, acetylene, hydrogensulfide, ammonia and air (a pseudo-component treated with properties averaged from N2 and O2). Thermophysical property data for VOCs must be provided by the user. The fluid components may partition (volatilize and/or dissolve) among gas, aqueous, and NAPL phases. Any combination of the three phases may be present, and phases may appear and disappear in the course of a simulation. In addition, VOCs may be adsorbed by the porous medium, and may biodegrade according to a simple half-life model or using a more complex biodegradation reaction model.

TMVOC allows flow system in any one of the three phases (Gas, aqueous, oil/NAPL phase) or combination of them. Each phase flows in response to pressure and gravitational forces according to a multiphase extension of Darcy’s law, which includes effects of relative permeability and capillary pressure between the phases. Transport of the mass components may also occur by molecular diffusion in all phases. Multiphase diffusion is treated in a fully-coupled manner that can cope with diffusion of phase-partitioning components under conditions of variable phase saturations.

Depending on thermodynamic conditions and relative abundance of the different components, the fluids may exist in seven different phase combinations, as shown in Figure 23. Arrows denote routes for appearance or disappearance of phases that are checked after each update of thermodynamic conditions during the Newton-Raphson iteration process.

  1. Specifications

A summary of TMVOC specifications is given in Table 15. The default number of mass components is 2 (water and one gas). The maximum number of components is 20. The flow system can be in any one of the seven phase conditions. The component numbering is based on the number of gases, hydrocarbons, and tracers in the model. The first component is always water; the second component is gas1; the third component is gas2; the first hydrocarbon is number (1+NumGases+1) component; the second hydrocarbon is number (1+NumGases+2) component; the first tracer is number (1+NumGases+NumHyC+1) component; the second tracer is number (1+NumGases+NumHyC+2) component. The selection of primary variables uses following rules:

(1) The first primary variable is always the phase pressure (if multiple phases, it is the lighter phase pressure).

(2) The last primary variable is always temperature.

(3) For single phase conditions, primary variables from second are mole fractions of the gases, mole fractions of the hydrocarbon, and mole fraction of tracers, respectively.

(3) For two phase conditions, the second primary variable is the saturation of lighter phase. Primary variables from third are mole fractions (in the lighter phase) of the gases (starting from second gas), mole fractions of the hydrocarbon, and mole fraction of tracers, respectively.

(4) For three phase condition, the rule (3) is applied except the primary variable for mole fraction of last hydrocarbon is replaced by water saturation ( SaS_aSa​) .

Table 15 Summary for TMVOC

Specification
Parameters

Components

(1) Water

(2) Gases (1, ....., NumGas) (2-NumGas are optional)

(3) Hydrocarbons (1, ......, NumHyC) (optional)

(3-7) Tracer1-Tracer5 (optional)

Phase condition and its state name and index

(1) Gas, GAS

(2) Aqueous, AQU

(3) Oil, OIL

(4) Aqueous-Gas, G_A

(5) Oil-Gas, G_O

(6) Aqueous-Oil, O_A

(7) Three Phase, GOA

Primary variables

Optional process modeling

Molecular diffusion, Wellbore simulation, Biodegradation reactions, and non-isothermal simulation.

  1. Specific input requirements

IE(402) selection of air solubility model

=0 use model by D'AMORE AND TRUESDELL (1988), CRAMER (1982)

=1 use constant air solubility, HC=1.0×10−10HC=1.0 \times 10^{-10}HC=1.0×10−10

IE(403) selection of component fraction unit in the simulation. The initial conditions must be consistent with the selection, and the source/sinks are also consistent with the selection. If mass fraction is selected, injection or production rate must be in kilograms. Otherwise, the injection or production rate must be in moles

=0 use mole fraction (default).

=1 use mass fraction.

Details for assumptions, physical processes, fluid properties, and phase partitioning for TMVOC module can be found in the original .

See

TMVOC simulation requires definition gas, hydrocarbon, and tracer components. They are defined through keywords "", "" and " ", respectively. Details can be found from the inputs for these keywords. If no gas is defined in a simulation, the default gas "AIR" will be included as the gas component. TMVOC does not require other specific input. The following inputs may be required, but are optional.

7️⃣
TMVOC user manual
GASES
CHEMP
TRACR
Table 28
Figure 23. Phase compositions and phase changes considered by TMVOC. The phase designations are: g - gas, w - aqueous, n - NAPL