-
initialization of water composition
Hi, Recently, I tried to use toughreact (EOS1) to simulate fluid-rock interaction including dissolution and precipitation. But it always reminded me “convergence problem in initialization of water…
-
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…
-
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.…
-
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.…
-
How do I conduct a batch modeling of NaCl-water with primary minerals
Hello everyone, I am trying to conduct a batch modelling of NaCl-water with primary minerals in a rock sample over 10-year period to determine the brine composition and primary species in the rock.…