TOUGH3 Category

Heterogeneous field simulation

I was simulating water injection into oil reservoir with a heterogeneous permeability/porosity field using TOUGH3-EOS8. The heterogeneous permeability/porosity fields were realized by modifying the…

Error in installation

Can somebody help in rectifying the error encountered as shown in attached figure. Thank you in advance

PETSC configure command failed-when compile TOUGH3

I got the "PETSC configure command failed" error at compiling TOUGH3. I need some help. Thanks. cd /Users/wuhao/TOUGH3/TOUGH3v1.0/TOUGH3-Code/esd-tough3/build/tpl-build/tpls/petsc/petsc-3.7.…

GENER BLOCK : water pumping EOS7r two phase problem

Hello everyone, I want to pump water from a certain elements using EOS7R. I'm confused between using: 1. "COM1" as a type and a negative value of the rate (< 0) to mean that I extract water from this…

Is there a way to specify the component being produced from a source well?

My producing well is producing fluid through the ESP and venting gas through the casing annulus. Both flows are controlled independently. I am trying to history match the test and would like to…

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

Mass balance of ECO2N

I was simulating constant CO2 injection into a 2D aquifer (salt not considered). For such 2-phase 2-component system, there should be 2 mass balance equations for CO2 and H2O, respectively,…

TOUGH3 vs TOUGH2-MP Printout of Mass fractions in single phase

Our understanding is that adding all component mass fractions of a phase add up to one. In TOUGH3 (or at least using eos3), it would appear that this is not the case.…

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

Simple model test with ECO2N

I was testing the ECO2N module compiled by myself (Linux version). I built a simple model with 2 grids. I set pressure difference of 2MPa and gas saturation difference of 0.…

Hysteresis and non-wetting phase entrapment saturation

Hello everyone, First of all, thank you all for your help. I have a question related to gas permeability. In a wetting process where \( k_g(1 - S_{\delta_{gr}}) = k_g(S_{l^*}) = 0 \),…

cannot find -laztec

Hi folks, I am doing a clean TOUGH3 installation on a Linux system and want to use the PETSc solver. I have no interest in using the AZTEC solvers. But it seems I cannot finish building TOUGH3…

T2WELL deals with time scale inconsistencies

Hello, everyone. I would like to ask how T2WELL deals with the problem of inconsistency of wellbore and reservoir models in time scale?

T2WELL and TOGA coupling

Hello, everyone. Recently our group would like to do about CO2 injection into reservoir for CO2 flooding and sequestration based on wellbore-reservoir coupling model.…

Does TMVOCBio in TOUGH3 have any option for 3-phase hysteresis

Hi all, What I recall is that TMVOC did not have Hysteresis options for three-phase capillary pressure and relative permeabilities (see the link). Seems it is still the case in TMVOCbio (in TOUGH3).…

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…

how to debug Tough3 in vscode

Hi all, When configuring, I changed the -- build type to debug, which can generate an executable file called tough3-eco2n in the tough3-install folder. Then,…

Metis and parallel computing in Tough

Hi All, I found many descriptions about Metis in the User's guide of Tough2-MP. However, there is little description of Metis in Tough3. So, is Metis used in Tough3 for parallel computing? Moreover,…

MESHMAKER CORNERS

Good Morning All, I believe that a version of MESHMAKER (Meshmaker 2?) exists that writes a CORNERS attribute defining the nodes of an unstructured grid.…

How to visualize TOUGH MESH file in pyvista?

Hello TOUGH3 Forum, I got the above files, and want to visualize it with some python packages, e.g., pyvista or pytecplot. Now I can only read "MESH" file with the following code through toughio:…

call TOUGH in jupyter notebook and PyTOUGH

Hi All, I am using PyTOUGH to automate TOUGH simulation. However, when I try to call TOUGH in Jupyter notebook with this command subprocess.call(["path_to_my_tough3_executable(tough3-eos1)",…

Run ends early

Hi, I was working on a larger block mesh using a doublet configuration for an aquifer, over 10 years. When I made the mesh finer, the simulation timed out and quit after 16 hours.…

zero pore pressure boundary condition

Can TOUGH3 define pore pressure Dirichlet boundary condition as zero values? In order to validate the coupled THM code, I have to compare the numerical solution to analytical solution.…

TOUGH3-FLAC3D time step setting

Hi, I am trying to couple TOUGH3 and FLAC3Dv7.0 through python scripting engine built in FLAC3D software. I am confused about the time step control after reading some published papers (e.g.…

MAX.RES. is small but DELTEX keeps decreasing

Hi All, I am running TOUGH3 EOS7C. The simulation seems stuck at extremely small time step. Please see below for the failure message. You can see that the residual is already very small but the time…