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. PREPARATION OF MODEL INPUT
  2. Keywords and Input Data

CHEMP

CHEMP provides data for calculating the thermophysical properties of the NAPL/chemical. The units of many of these constants are not standard metric units, and care must be taken to ensure that the appropriate units are used. Most of the data used in this block can be taken from Appendix A of Reid et al. (1987), where the same units are used as here.

Record CHEMP.1

Free format, or Format(I5)

NumHyC

NumHyC number of organic chemicals for which data are to be entered. NumHyC must be <19.

Record CHEMP.2

Free format, or Format(A20)

HcNames

HcNames name of the organic chemical.

Record CHEMP.3

Free format for 5 parameters, or Format (5E10.4)

TCRITM, PCRITM, ZCRITM, OMEGAM, DIPOLMM

TCRITM chemical critical temperature, K.

PCRITM chemical critical pressure, bar (1 bar = 10510^5105 Pa).

ZCRITM chemical critical compressibility.

OMEGAM Pitzer's acentric factor for the chemical.

DIPOLMM chemical dipole moment, debyes.

Record CHEMP.4

Free format for 5 parameters, or Format (5E10.4)

TBOILM, VPAM, VPBM, VPCM, VPDM

TBOILM chemical normal boiling point, K.

VPAM chemical vapor pressure constant from Reid et al. (1987). If VPAM ≠ 0 use Wagner correlation (Reid et al., 1987), for VPAM = 0 use Antoine Correlation.

VPBM chemical vapor pressure constant from Reid et al. (1987).

VPCM chemical vapor pressure constant from Reid et al. (1987).

VPDM chemical vapor pressure constant from Reid et al. (1987).

Record CHEMP.5

Free format for 5 parameters, or Format (5E10.4)

AMWTM, CPAM, CPBM, CPCM, CPDDM

AMWTM chemical molecular weight, g/mole.

CPAM chemical ideal gas heat capacity constant from Reid et al. (1987).

CPBM chemical ideal gas heat capacity constant from Reid et al. (1987)

CPCM chemical ideal gas heat capacity constant from Reid et al. (1987)

CPDDM chemical ideal gas heat capacity constant from Reid et al. (1987)

Record CHEMP.6

Free format for 5 parameters, or Format (5E10.4)

RHOREFM, TDENREF, DIFV0M, TDIFREF, TEXPOM

RHOREFM reference NAPL (liquid) density, kg/m3.

TDENREF reference temperature for NAPL density, K.

DIFV0M reference binary diffusivity of VOC in air, m2/s.

TDIFREF reference temperature for gas diffusivity, K.

TEXPOM exponent for calculation of chemical diffusivity.

Record CHEMP.7

Free format for 5 parameters, or Format (5E10.4)

VLOAM, VLOBM, VLOCM, VLODM, VOLCRITM

VLOAM liquid NAPL viscosity constant from Reid et al. (1987).

VLOBM liquid NAPI viscosity constant from Reid et al. (1987).

VLOCM liquid NAPL viscosity constant from Reid et al. (1987). If VLOAM and VLOBM = 0, VLOCM is reference NAPL viscosity in units of cP (1 cP = 10-3 Pa-s).

VLODM liquid NAPL viscosity constant from Reid et al. (1987). If VLOAM and VLOBM = 0, VLODM is reference temperature for NAPL viscosity in units of Kelvin.

VOLCRITM chemical critical volume, cm3/mole.

Record CHEMP.8

Free format for 4 parameters, or Format (4E10.4)

SOLAM, SOLBM, SOLCM, SOLDM

The chemical solubility is calculated from the polynomial SOLUBILITY = SOLAM + SOLBM*T + SOLCM* T2T^2T2+ SOLDM* T3T^3T3. If data for the solubility as a function of temperature are available, then SOLAM, SOLBM, SOLCM, and SOLDM should be calculated from a polynomial fit of the data. If such data are not available (the usual case), the solubility will be assumed to be constant, and SOLAM should be set equal to the known solubility, with SOLBM, SOLCM, and SOLDM set equal to 0.

SOLAM constant for chemical solubility in water, mole fraction.

SOLBM constant for chemical solubility in water, mole fraction/KKK.

SOLCM constant for chemical solubility in water, mole fraction/K2K^2K2.

SOLDM constant for chemical solubility in water, mole fraction/K3K^3K3.

Record CHEMP.9

Free format for 3 parameters, or Format (3E10.4)

OCKM, FOXM, ALAMM

FOCM default value for fraction of organic carbon in soil, used for all domains for which no specific value is provided in record ROCKS.1.1.

ALAMM decay constant for biodegradation of VOC, s-1. Biodegradation is assumed to take place only in the aqueous phase, and to follow a first order decay law, MVOC(t) = MVOC,0 * exp (-λ t). The decay constant λ = ALAMM is expressed in terms of the half life T1/2 of the VOC as follows: λ = (ln 2) / T1/2. Default is ALAMM = 0.

Repeat records CHEMP.2 through CHEMP.9 for a total of NumHyC (up to eighteen) different organic chemicals.

Used in: TMVOC

EXAMPLE:

CHEMP

2 //Number of organic chemicals

BENZENE // Name of the first chemical

562.2, 48.2, 0.271, 0.212, 0.0 //CHEMP.3 for BENZENE

353.2, -6.98273, 1.33213, -2.62863, -3.33399 //CHEMP.4 for BENZENE

78.114, -.3392E+02, 0.4739E+00, -.3017E-03, 0.7130E-07 //CHEMP.5 for BENZENE

885., 289.00, 0.770E-05, 273.10, 1.52 //CHEMP.6 for BENZENE

0.4612E+01, 0.1489E+03, -.2544E-01, 0.2222E-04, 259.0 //CHEMP.7 for BENZENE

0.411E-03, 0.000E+00, 0.000E+00, 0.000E+00 //CHEMP.8 for BENZENE

0.891E-01, 0.001, 0.0 //CHEMP.9 for BENZENE

n-DECANE // Name of the second chemical

617.700, 21.200, 0.249, 0.489, 0.00 //CHEMP.4 for n-DECANE

447.300, -8.56523, 1.97756, -5.81971, -0.29982 //CHEMP.4 for n-DECANE

142.286, -7.913E+0, 9.609E-1, -5.288E-4, 1.131E-7 //CHEMP.5 for n-DECANE

730.000, 293.000, 1.000E-5, 293.000, 1.600 //CHEMP.6 for n-DECANE

0.0, 0.0, 0.5900, 293.000, 603.000 //CHEMP.7 for n-DECANE

3.799e-7 //CHEMP.8 for n-DECANE, constant solubility

0.000 //CHEMP.9 for n-DECANE. The first parameter is 0.0, others set to default values.

PreviousCBMDANextCOFT

Last updated 1 year ago

Two options are available for calculating the NAPL liquid viscosity. The liquid viscosity constants VLOAM - VLODM for the desired NAPL may be assigned data given in Table 9-8 of Reid et al. (1987), and the viscosity will be calculated using a polynomial fit to actual viscosity data (Eq. 4.3.11 in). Alternatively, VLOAM and VLOBM may be set equal to 0, and VLOCM and VLODM are assigned equal to a reference viscosity and a reference temperature, respectively. In this case, the viscosity is calculated from a more general (and less accurate) empirical correlation (Van Velzen et al., 1972).

OCKM chemical organic carbon partition coefficient Koc (see Eq. 4.4.2 in ), m3/kg.

8️⃣
TMVOC user manual
TMVOC user manual