-
TOUGH+HYDRATES convergence issue in phase change
I want to calculate the depth interval at which gas hydrates forms, given P and T initial conditions, fixed on top. The idea is to run the model until it reaches the equilibrium,…
-
EOS3-sink/source rates are piecewise constant, how to set up?
As described in the title, I want to set up sink/source rates which are piecewise constant,as shown in figure attached, I tried to set up the input file, but i'm confused about it?…
-
Heat generation in EOS3
Hello, I am a relatively new TOUGH2 user. I have developed a model for heat and fluid flow in a basin with a lake. I am attempting to have atmospheric boundary conditions are the top of the model…
-
Volume and Mass Balance Error in flow.out
Hello, everyone. I am using EWASG model to simulate chemical reaction (Ba+2 + SO4-2 = BaSO4) in porous media. Although I define the initial water and assign it to the matrix,…
-
TOUGHREACT_Debugging
Hello Everyone: I am pretty new to TOUGHREACT. I am trying to build a model to simulate Ba++ and So4-- reaction in porous media. The attachments are the input files I modified based on P8_EOS1 from…
-
SAVE file not consistent with OUTPUT
Hello, I am running TOUGH3 EOS7C for a gas reservoir modeling. The SAVE file will be used as INCON but I realized that some SAVE files are not consistent with the OUTPUT file.…
-
GENER block injection fluid
Hello, There is an option for injection component on the TOUGH+HYDRATE GENER data block defined by COM1, COM2, COM3, etc. However it doesn't specify how we can inject a particular fluid into an…
-
applying heat load in BHE and varying input T?
Hi all, I am using T2Well with TOUGH2 to model a borehole heat exchanger and was wondering if it was possible to apply a heat load to the BHE? I have tried something of the form below:…
-
T2WELL/ECO2N problem with mass rate
Dear TOUGH2 community users, I am new to T2Well and I am trying to simulate pure gas-CO2 injection with several radial wellbore cells and only one big reservoir cell,…
-
3D Visualization Software
Hi, guys! Some questions have been bothering me for a long time, I need help! OK, I have built a 3D model. But when visualizing with tecplot,the result is not what I want.…
-
Time-Dependent Material Properties
Hello, Is there a way to implement time-dependent material properties (block ROCKS) in TOUGH3? I am interested in assigning time varying permeabilities (prescribed as a function of time).…
-
Incorporating dip angle in meshmaker
Hello, I am trying to model a gently dipping aquifer, and am setting up the grid using the online meshmaker utility. The input field called "Angle between Y-axis and the horizontal" is meant to help…
-
Can TOUGH3 model temperature below zero degree?
Hi All, I am new to TOUGH and plan to use it in modeling a borehole thermal energy storage system. However, the ground temperature surrounding the system may drop to below zero degree.…
-
FOFT keyword no printing
hi. I have been trying to select the elements of my interest to print out in the forward.out or another time series file like in the manual examples by using FOFT in my forward input file but…
-
TOUGH+HYDRATE pre/post scripts
If there is any active user on T+H, is anyone interested in discussing pre and post processing options. After struggling for some time, I have managed to put some scripts in python to process both…
-
T2Well questions
Hi there, Starting to play around with T2Well and having good fun. ;) Two things (at least) that I didn't understood: 1- ROCKS data block: the domain name for wellbore cells starting with 'x'.…
-
A question about the additional mechanism under the Minerals entry in the chemistry.inp file
Dear All, Hello everyone! I just started learning toughreact software. At the moment, I am a bit confused about the description of the additional mechanism under the Minerals entry in the chemistry.…
-
EOS7 - Brine Definition & Tracing flow of brines
Hello everyone I am looking for more clarity on the definition of brine in EOS7. In the TOUGH2 user's guide V2, it states that, "The brine is modeled as NaCl solution,…
-
Generation Time Table
Dears, I have this problem while constructing the generation time table using (GENER). I have a table of 5 time-points (rows) which is illustrated below. When I run the code, it crashes.…
-
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.…
-
injecting VOC mass through virtual cells
Hi everyone, I am trying to simulate a 1D diffusion scenario where the domain is a column comprised of 200 grid cells (domain is 1m*1m*20m with 200 cells in z-direction) and every grid cell has been…
-
Two time steps with restart gives different results than one without restart
Hello, Do you know why a restart during 6 hours (for example) after one time step (6 hours for example) is different from one time step of 12 hours? In my simulation with eos4,…
-
What is the use of 'aquatic complexes'
I am wondering how aquatic complexes may change the simulation results. below are my findings: 1. aquatic complexes, if declared in 'chemical.inp', are assumed to be at equilibrium with basic…
-
Compilation TOUGHREACT on Windows using gfortran and Cygwin
I am having some problems with compiling the original source code of TOUGHREACT. I have tried using all of the compiler flags and the ones suggested on this forum.…
-
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.…