TOUGHREACT Category

Error message: Components missing:h2o

Hello All! I am trying to run a batch equilibrium simulation with several mineral phases including goethite and others having h2o in the reactions equation.…

how to set initial temperature distribution

Hi，tough community users, I want to prepare the initial or stable pressure and temperature conditions for EGS. INDOM was introduced into my flow file. The UPPE1 reservoir is the top mesh having 1.…

About the fluid velocity in the governing equation

Hi, everyone, I'm currently studying the governing equations underlying TOUGHREACT, particularly in the context of the ECO2N module. The entire solution process is as follows： （1）For the fluid flow…

During CO2 sequestration at basalt formation, time step size rapidly decrease

hello, everyone. I am graduated student, during master degree. I'm constructing simulation model about basalt CO2 mineralization. I am utilizing TOUGHREACT 1.2.…

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

What is the popular way to visualize 3D results from TOUGHREACT

Hi, I have unstructured grid created by GMSH, what is the popular way to visualize 3D results from TOUGHREACT I tried TOUGHIO, but it won't read my flow.out properly

total concentration of sulfate increasing over time

Hello, I ran a chemical initialization simulation for 10 million years to achieve steady state. In order to figure out how much dissolution and precipitation I have in my system,…

low injection rate cannot work in 2D model

Hi, everyone. I am so sorry to interrupt you again for a same question. I want to simulate the dissolution and precipitation of sandstone in 2D model, which were induced by meteoric water injection.…

How to show the permeability and porosity changes？

I'm using TOUGHREACT1.2/ECO2N with Petrasim,and i found that my result shows that the permeability and porosity is the same as my initial value,but I have opened the KCPL to consider feedback effects…

How to set reasonable injection rate?

Hi, everyone. I want to simulate the dissolution and precipitation induced by meteoric water injection. Firstly, I acquired the annual rain precipitation about 100 mm/a.…

How could charge diffierent boundary water at same model?

I set two types of boundary water in same model. And I want to charge first boundary water into reservoirs during 100-50 Ma, then charge the second one by same boundary elements sicne 50 Ma to…

TOUGHREACT model boundary grid Settings

Hello, I want to set up flux-free boundary conditions in TOUGHREACT, how do I set up the related mesh?

Tough code

Hello, I am a doctoral student at Xi'an Jiaotong University. Due to the regional restrictions, it is unable to buy the ToughReact code, can you share a ToughReact source code? Very grateful

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

V1.2 TOUGHREACT ECO2N IN PetraSim:Problems with initialization for flow and chemical field

Hi, I am new to TOUGHREACT and using v 1.2 ECO2N on PETRASIM to model CO2 Sequestration. And I have created a 2-D isothermal rectangle model and designed a temperature gradient of 3.…

TOUGHREACT - Dissolution and precipitation

Hi all, I am a new user of Toughreact. I am testing a real case dataset, whereas I want to quantify the amount of CO2 that will be mineralized in a mafic formation.…

Therakin.dat question

Hello, I would like to know what some of the numbers in therakin are. 'dolomite-2' 184.401 64.365 4 -2.0000 'h+' 1.0000 'ca+2' 1.0000 'mg+2' 2.0000 'hco3-' 'dolomite-2' 6.9420 5.7530 4.2300 2.6830 0.…

time.out won't output the data of the last timestep

So my last timestep is 8.64e3 sec, but the time.out the last data in time.out is Is there any setting? I don't want to go the .tec file to find the cell by XYZ, thanks!

About the regular and radial mesh

Hi everyone, I was using TOUGHREACT V1.2 ECO2N to run a simple injection simulation，and When I set up multiple injection points(such as 5 points and each is spaced 20m apart, that is,…

FLOW VECTOR CANCELLING EACH OTHER FOR SOURCE/SINK TERMS

Dear, I notice that for all the source and sink elements flow vector with opposite direction are cancelling each other. To avoid this I would like to output the absolute total flow,…

Toughreact and eos5

Hello everyone, I am new at using Toughreact and would like to know if there is any version of Toughreact that works with EOS5 (Water, Hydrogen). If not,…

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

Restarting TOUGHREACT Simulations

I'm trying to restart TOUGHREACT simulations by renaming SAVE to INCON and savechem to inchem. In between simulations, I'm adjusting the initial waters associated with certain elements by changing…

Error when trying to output the time series

So, I have 40000+ elments, and I want to output all of them's data, so I put all the elements in the solute.inp, for example But then I have the error Elements for which to output time series:…

About the toughreactEco2n v1.2 not have the “savechem”files

Hello everyone, I was using TOUGHREACT V1.2 to run a initial progress of HC field(no thermal field). There were some problems when adding the chemical field,…