TOUGHREACT Category

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

modify output flowdata.tect

Dear all, I 'm using EOS1 with TOUGHREACT V4.13omp and I would like to modify the source code to output in the flowdata.tec file the density of the different phases with their specific internal…

runtime in Problem 6 2-‐D Radial Model for CO2 Sequestration in Deep Saline Formation from TOUGHREACT V4

Hello everyone! I am working on Problem 6 2-‐D Radial Model for CO2 Sequestration in Deep Saline Formation from TOUGHREACT V4.0. I already got the results for the first 10 years...…

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

Describe the GENER block in context of the sample problem p6_eco2n

Could someone kindly provide an explanation of the GENER block, along with a discussion of its defined parameters, focusing on sample problem "p6_eco2n_2D_radial"?

Don't understand what causes this error

The first picture is the error from .exe, the second error is in the VERS, and only flow.out has been created flow.out:…

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…

SMCO2 is not increased

Hello I am using Toughreact v3 for simulation for CO2 injection and our simulation runs well, however, CO2 is not mineralized after the simulation has finished.…

Last bock in CO2TAB file

Hello everyone: I have difficulties to understand the last block in CO2TAB file and search for help from here. I know the previous blocks are pressure, temperature, viscosity,…

TOUGHREACT/ECO2N:CANNOT FIND PARAMETERS AT ELEMENT * 1*

Hi Everyone, I'm new at using ECO2N of TOUGHREACT. And the following problem is encountered when I used PetraSim to simulate. It's a simple model. And I model a strata from -2260m to -2900m below…

CONVP block in solute.inp

Hello everyone: I have a question about implementing CONVP (chemical convergence bounds for time stepping) in solute.inp. I firstly make MOPR(10)=2 to make simulator read CONVP block.…