TOUGH3 Category

TOUGH3 installation issue

I realized many people have difficulty in compilation of TOUGH3. After polishing the installation cmake scripts and reorganizing the installation package, compilation of TOUGH3 should be more easy,…

How to debug tough3 correctly

I want to debug this code. When configuring tough3-eco2n, the command used is"/configure.sh --build-type=DEBUG --eos=eco2n --no-x11 --petsc-dir=/home/zhigao/petsc-3.17.…

Coupled reservoir-wellbore model initialization (using T2Well)

Hello everyone, I have started using T2Well to simulate coupled reservoir-wellbore problems recently. While initializing the numerical model, I find the following problem;…

EOS5 but Helium in stead of Hydrogen

Hello, I am trying to simulate a year long helium injection with varying rates using EOS5. I'm wondering if I would need to change something in the source code for this to run correctly,…

Two phase flow (Diffusion vs advection) with mesh sizes

Dear all, dear Finsterle, I have a question regarding the mesh size and its impact on diffusion in two-phase flow. I'm curious to know if the mesh size affects the diffusion from element0 to element1.…

EOS7MG-A new TOUGH3 module for multiple gas mixture simulation

Recently, we have developed a new TOUGH3 EOS modules named ESO7MG. It is designed for modeling water, brine, and multiple gas species and their mixtures including hydrogen in subsurface multiphase…

An error running the "CO2 leakage along a fault" example in toughio tutorial

Hi, When I was trying to run the INFILE generated by the "generate_model_parameters_input_file" script, I had this error message: "forrtl: severe (408): fort: (2):…

AZTEC solver - rapidly increasing temperatures going out of range

I am running an inj/prod well doublet in a geothermal reservoir with free gas at the top of the reservoir so using EOS7C (with 4-12 cores depending). The AZTEC solver crashes rapidly b/c the…

GOFT files only printing for one component (COM1)

Running EWASG and GOFT files are only printing for COM1 and not for COM3 (see my INFILE attached). I have two generation elements per node (the injection /production wells move two components,…

ECO2M: Automatic switch in injection constraint

Hello, I'm running ECO2M to simulate a CO2 injection. If it's possible to inject the CO2 under a pressure/BHP constraint, would it be then possible to simulate the injection initially under…

Display temperature and pressure change in 2D plot

Hi everyone, I am trying to visualize pressure and temperature results for 3D model that has initial pressure and temperature values described by hydrostatic pressure and geothermal temperature…

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…

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.…

Questions Regarding Local Refinement and Mesh Size Limitations (number of elements)

1. Local Mesh Refinement I am importing my mesh from ECLIPSE and not using the MeshMaker utility. I would like to apply local refinement, similar to the approach shown in the snippet below.…

How to implement anisotropic permeability tensors and unstructured mesh in TOUGH3?

Hi everyone, I'm working on a model involving complex geological structures with anisotropic flow behavior, and I'm planning to use TOUGH3 because of its support for full permeability tensors and…

Invitation for TOUGH4 testing

TOUGH4 is the next generation of TOUGH codes, featuring numerous advancements and enhanced modeling capabilities. More information can be obtained from the following webpage: https://lbl-2.gitbook.…

*** PETSC Solver failed to converge *** in parallel mode (TOUOGH3_v1.12- EOS1)

Hi, I have a problem with code parallelization. The code runs smoothly on my laptop, but as soon as I move to the HPC and use multiple processors, it says *** PETSC Solver failed to converge ***.…

Infinite Acting Aquifer

I have a TOUGH3 model (EOS7C) that is showing signs that the model extent is not wide enough in the East and South directions (pressures increasing over time).…

TOUGH3 Installation Fails Due to BLAS Package URL Change

I am writing to seek your assistance with an issue I encountered during the installation of TOUGH3. The installation process fails due to a recent change in the URL of the BLAS package,…

Mass balance assessment with TOUGH3

I am looking to validate the observed mass balance agains the theoretical (initial + source/sink - migration form domain into boundary cells) in a TOUGH3 with Neumann boundary conditions.…

How should the enthalpy be correctly specified when simulating the process of air injection into an aquifer using TOUGH3/EOS3?

Hi I am currently attempting to conduct a simulation for CAESA, specifically modeling the process of air injection into an aquifer. Due to persistent convergence issues,…

possibly bug in viscosity calculation in TOUGH3_EWASG

Dear developers, First of all, I would like to say thank you for developing this excellent codes. Here I am reporting a tiny bug I found in viscosity calculation in TOUGH3 EWASG,…

TOUGH3 output and ext2

Hello all, I am using the famous ext2.exe file to extract data from TOUGH3-eos3 output file but it is not working (although it works perfectly fine with my TOUGH2 outputs).…

Use of AMRES in TOUGH3

It appears that the `AMRES` parameter is defunct in TOUGH3 (see Main_Comp.f90). However, I am having issues of exceesive residuals which break PETSc. Are there any plans to implement a fix for…

T2Well/ECO2N mesh software recommendation

Hi everyone! I'm new to using T2Well/ECO2N and was wondering which meshing software you recommend for creating an anticline structure with different layers? I've been exploring MeshMaker, PetraSim,…