-
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.…
-
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.…
-
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?
-
Thermodynamic Database for CO₂ Storage in Icelandic Basalts
Hello everyone, I've been running sample problem 5 on CO₂ storage in sandstone using TOUGHREACT V3.32. I'm interested in running this problem considering a basaltic formation, however,…
-
TOUGHREACT/EOS5 serpentinization
Hello All, Is TOUGHREACT / EOS5 able to simulate the process of serpentization that can lead to the production of natural hydrogen? I know that EOS5 is intended for transporting hydrogen in the…
-
Simulation of the non-isothermal analysis of co2 injection
Hi, everyone i am now tyring to simulate the co2 injection under the non-isothermal situation byTOUGHREACT V1.2 ECO2N). For simplicity, i cosed the reactive transport in this simulation and focus on…
-
How to perform simulations of water-rock interactions to obtain near-equilibrium water chemistry.
Hello everyone, I am currently studying TOUGHREACT and using Petrasim for preprocessing and postprocessing, with the aim of researching CO2 geological storage issues.…
-
Output of permeability and initial value of porosity and permeability
Hi,everyone i am using TOUGHREACT v1.2/ECO2N, and i have some issues (1) i find that the file 'solid.out" only shows one permeability, what orientation it is ?…
-
HELP!!! Regarding the issue of being unable to simulate flow given multiple rock types~THX!
I am currently using TOUGHREACT V3 software and ECO2N module for flow simulation. When all the first domain is given in ELEME, the simulation is successful;…
-
The TOUGHREACT V4.13's calculation time is long
Hi everyone, I have been trying to run TOUGHREACT V4 Deep Salt Formation 2-D Radial Model Problem 6 with the instructions in the manual. However, I encountered the following two troubles:…
-
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.…
-
how to use RESTART to inherit the minerals fractions and change P T?
Hi, My model is 1D with 10*10*400 meter. Firstly, I renamed the save and savechemical to INCON and inchemical in the first simulation, and copied these two files to the second simulation folder.…
-
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…
-
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…
-
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…