-
Convergence problem chemistry
Hello, everyone, I was running a water-rock interaction model with Toughreact. The runlog.out file prompted "Warning: chemistry did not converge at node A1114 (routine NEWTONEQ),…
-
Visualize data using MATLAB - Github link inside
Hello TOUGH/TOUGHREACT community, I would like to share some code I developed on MATLAB to visualize pressure profile, gas saturation, pH, etc. This code only works for 2D data. https://github.…
-
Compiling and Running TOUGHREACT v3.32 using Ubuntu/Linux
Hey everyone, I've been able to get the newest version of TOUGHREACT working with my workstation, which has Linux/Ubuntu (v 17.04) as its OS. Setting things up didn't go without hitches, though,…
-
No change in minerals vol. fraction
Dear All, I am a beginner in TOUGHREACT simulation and I have issues with getting results from simulation. This simulation performs simulation and there is no error indicated in the runlog file,…
-
Lithofacies random fields in TOUGHREACT?
Hi, I am trying to apply random fields both in flow and rock chemistry. I know I can do one for the flow field using permeability modifier (PMX) option equipped in TOUGH2 (flow.…
-
What does negative H+ concentration mean in TOUGHREACT?
I am running some simuations, the results seem to make sense (including pH). But why does the h+ concentrations are negative? Is there a special calculation to convert h+ to pH?
-
Confusion about CTOT and CGUESS
Hello, I am quite confused about how to provide initial chemistry to TOUGHREACT. If I want to tell TOUGHREACT that my pH is 6 and calcium concentration is 2mol/kg H2O, how should I proceed?…
-
Stuck on Performing Simulation, and chdump.out stops writing resutls
It shows multi TIME=0 CHEM.ITRERATIONS is 1 for all of them So does this mean it doesn't converge on CHEM? Thanks
-
New to Tough simulator, question about Mesh
Hi, I am new to Tough simulator, and I have built a model in Gmsh, it would be great if anyone tells me what material I should look at to use my model. Thanks!
-
Can TOUGHREACT automatically calculate injection rate into each element?
Hello, I was wondering if TOUGHREACT can automatically calculate the injection rate into each element without the user interfering? For example, I would like to injection 90 kg/s into a layer that is…
-
TOUGHREACT ECO2N
I want to use TOUGHREACT V3.0 ECO2N simulate the mineral changes caused by the co injection of different concentrations of SO2 and CO2, mainly observing gypsum precipitation.…
-
Stopping in main due to flow nonconvergence
Hello everyone, I am using ECO2N to simulate CO2 injection into a formation with relatively low permeability. I am having trouble with an error popping up.…
-
Error in reading INCON - Problem 5, eco2n_1D Radial
Hello, In problem 5, eco2n_1D Radial, I am using the input files provided (flow, solute, chemical, CO2TAB, .dat, and executable) to run a q=0 kg/s flow in order to find the equilibrium conditions…
-
An issue (experience) with very large mesh
Hi all Just to share my experience. I have a very large mesh (6E5 grid cells). Surprisingly the name of some elements start with "END" (like END 1, .. END10 etc).…
-
TOUGHREACT compilation using gfortran MinGW-error message
Dear Community, I would like to compile TOUGHREACT (version: TR-OMP_Source_v3.3.2-OMP) using gfortran with MinGW. I typed make eos1 -f makefile_gf, and had the attached error message.…
-
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.…
-
Sequestered mass of CO2 not calculated/printed out by ECO2N
Hi there, I'm having a look at the example P5_eco2n (TOUGHREACT v332), on CO2 geological sequestration. To my surprise, the sequestered CO2 (output variable 'SmCO2') appears to be zero regardless of…
-
Looking for training sessions of Toughreact using Petrasim
Hi Everyone, I am new user of TR and I have research project that investigates CCS in Kazakhstan. Since the project can cover the training expenses, I would like to ask if anyone can explain how to…
-
set different diffusion coefficient for different layer
As far as I know, only one solute diffusion coefficient can be set for all aqueous species,which is not practical. Is it possible to set different diffusion coefficient of one solute for different…
-
Large runlog.out file
Hi I am having problems with large runlog.out files causing my disk size to get full quickly. I have made sure that the parameters responsible for this in REACT block are not made greater than zero…
-
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,…
-
Nonconvergence of TOUGREACT v332 (ECO2N)
I am performing 2-D radial simulation using TOUGHREACT (ECO2N). When I carried out the simulation using v332, the message "Stopping in main due to flow nonconvergence" appeared after several time…
-
Failure restarting Toughreact simulation.
Hi Everyone, I want to run a model with the geochemical data of a previous model. I follow the manual instructions changing the names of the SAVE and savechem files to INCON and inchem.…
-
Initial Geochemical Equilibrium between Minerals and Brine
Hello there, I'm simulating CO2 injection into sandstone aquifer of synthetic brine, NaCl. I've read in several papers in which TOUGHREACT w/ ECO2N is used and something is mentioned like "batch…
-
TOUGHREACT ECO2N
Hi, I'm quite new in using TOUGHREACT. Can anyone please clarify one question I've been having for quite a while now? The system I'm trying to model is just simple mixing between brine defined by the…