TOUGH3 Category

Non-isothermal solution crashes

I am running a CO2 single injection well into a dipping saline aquifer. The model seems to be running fine in isothermal conditions but crashes when non-isothermal is turned on.…

CANNOT FIND PARAMETERS AT ELEMENT *00001*

Hello, The input file I am using does not seem to be reading the INCON block. I am getting the following error: +++++++++ CANNOT FIND PARAMETERS AT ELEMENT *00001* XX(M) = 0.000000E+00 0.…

How to identify the cause of an incomplete simulation run

Dears, I have a problem to identify the cause of the incomplete simulation run for the attached input file. The simulation ends before the specified total simulation time without error messages!…

Questions about TOUGH3 setup for 3D random permeability fields

Our lab wants to use TOUGH3 for the CO2 storage reservoir simulations of 3D random permeability fields. We are not sure about how to set up the 3D anisotropic permeability correctly in TOUGH3;…

How to output viscosity in the output results

Hello every one, I want to know how to output water and CO2 viscosity at every grid block in the output results file. Thanks in advance

Infiltration into dry soil

Hello every one, I am new to work with TOUGH3 software. I tried to simulate 1D infiltration into a 1m high column of sandy clay loam with an initial head pressure of -8 m.…

Initialization (T2Well)

Hello everyone, I have the following question regarding initialization procedure for simulation using T2Well. I want to simulate the changes at wellhead of a monitoring well that extends to the…

T2Well Output

Dears, Can the results obtained from T2Well be output in Excel file format (like TOUGH3)? Thanks in advance

TOUGH3; long running time (3D model)

Hi everyone, I am trying to build 3D model for CO2 injection in underground saline aquifer that is upper bounded by impermeable strata which extends to the surface.…

gravity-capillary equilibrium in EOS3

hello, I am using EOS3 (TOUGH3-eos3 windows executable) to simulate a gravity capillary equilibrium along a 1-D column but I keep running to convergence failure issue.…

How to build windows executable file(.exe) through TOUGH3 code?

I tried to build an executable file by mingw32 on the UBUNTU platform, but I failed.

Problem with visualizing material type in ParaView

Hi every one, I have run initialization simulation to obtain initial pressure and temperature conditions that follow hydrostatic pressure gradient and geothermal temperature gradient, respectively.…

Visualizing grid blocks

Dears, I have the attached MESH file generated using MESHMAKER. I would like to know if there is a simple Python code that enables me to read and visualize the grid blocks in Pyvista.…

Specifying initial conditions for layered reservoir

Dears, I am trying to model CO2 injection in stratified saline aquifer reservoir. I have to provide initial conditions for my model considering hydrostatic pressure gradient and geothermal…

Numerical instabilities introduced after printout

I'm running a simulation with intermittent printouts using TOUGH3. I noticed that, after some of those printouts, numerical instabilities are seemingly introduced that stall the remainder of the run.…

Compiling against more up to date PETSC version

Hi All, Just a question regarding compiling TOUGH3 against a PETSC versions higher than the 3.7 branch. Has anyone tried this? Is it possible? I know at some stage a number of the PETSC function…

Phase transitions with tough3-eos7r

Hey everyone, I wanted to share an issue I've been having lately. I'm attempting to run a simulation with a vaporization front using tough3-eos7r. I've successfully run the problem with eos7,…

TOUGH3; Running time is too long

Dears, I am trying to run a simulation using 2D model for CO2 injection in stratified reservoir. I have 650 grid cells in the radial direction and 79 grid cells in the z-direction as illustrated…

TOUGH-3 Parallel compile on Windows 10 / Ubuntu terminal

Hi All, Has anyone managed to compile TOUGH3 in PARALLEL on Windows 10 with Ubuntu Terminal? I have followed Curt's instructions and managed to compile in serial (i.e.…

TOUGH2Matlab

Dears, I would like to know if TOUGH2Matlab can be used to post-process the results obtained using TOUGH3. I have tried to use TOUGH2Matlab to read the MESH and OUTPUT files obtained using TOUGH3,…

Printing outputs to a single file when running in parallel

Hey everyone, I'm running through a complicated simulation using tough3-eos7 with lots of phase shifts due to boiling. I noticed that, when I run the problem in parallel ("mpiexec -n <# of…

How to improve the stability of the numerical results without changing the grid discretization ?

Dears, I have the following simulation case for CO2 injection in stratified aquifer using 2D grid model. I have slight fluctuations in the pressure results.…

Question on DELV control

Hi again ! I have a question about the use of the Deliverability well control : What is the parameter HG in the GENER keyword used for ? What happens when it's defaulted ? (HG thickness of layer, m;…

INCON format in the *203t ECO2n samples

In the samples for Tough3 ECO2n, there are several cases whose name ends in 203t where the INCON keyword looks strange to me. For each element, there are 6 lines for the values of the primary…

Specifying material types for stratified heterogeneous reservoir

Dears, I am trying to model CO2 injection in heterogeneous aquifer that is composed of 23 layers with different porosity, permeability, and capillary pressure data.…