-
flow.out - residuals and equations
Hello, I am troubleshooting model. I have set MOP(1) = 1 and trying to figure out why my time-stepping continues to get reduced. I have automatic time-stepping turned on.…
-
Error when use MESHMAKER to generate XYZ mesh in TOUGHREACT ECO2N sample problem
Hello dear, I am trying to learn the ECO2N module in TOUGHREACT, so I began at the sample Problem 5 – 1-D Radial Model for CO2 Sequestration in a Deep Saline Formation.…
-
EOS1 non-convergence
Hello！ I am a beginner in TOUGHREACT. I ran it based on the input file provided by the author of a paper. However, when I ran toughreactV2.0, I found that the model did not converge and the program…
-
Reservoir Temperature Change
Hi all, In the current simulation I'm running, the reservoir temperature change is far less significant than I would expect (50 kg/s for 30 yrs with only localized grid block temp change).…
-
Carbonate cement distribution in 1D model
Hi, tough users I created a 1D model (EOS1) to simulate the carbonate cement. The initial temperature (15-26 C) and pressure (0.2-5Mpa, equal to 500 m depth) distribution are included in INCON file.…
-
TOUGHREACT and MODFLOW grid
Is TOUGHREACT AND MODFLOW grids interchangeable? Let's say I have a 3D MODFLOW grid for a project and now I want to use the same grid, potential same flow solution and transfer it to TOUGHREACT to…
-
EOS2 non-convergence
I'm new to TOUGHREACT and have built a simple mesh and created a temperature and pressure gradient with INCON. Now, I am trying to run flow.inp through TOUGHREACT (EOS2) before I move to the solute.…
-
How to set up well bore
Hi TOUGH community, I am trying to incorporate a well bore into my simulations in TOUGHREACT, and I am confused on how to set up my mesh. Please, kindly Assist.
-
TOUGHREACT Workshop Files
Hi, How can we get the files which are mentioned in the TOUGHREACT workshop handout here: https://tough.lbl.gov/assets/files/02/support/TOUGHREACT_Sample_Problems.…
-
how to set up initial chemical conditions?
Hello, I've just started using Toughreact, and I am wondering how to set up the chemical initial condition. Should I run a simulation to let TR calculate the brine chemical composition at equilibrium,…
-
CONSTANT pH PROBLEM
Hello everyone. After running a simulation (injection of CO2 into a saline aquifer), I observed that in my result, the pH remained constant at all times for my chemical initialization, injection,…
-
Batch Model with Brine Injection and Secondary Minerals Precipitation
I am attempting to set up a batch reaction model. Initially, the formation water equilibrates with calcite, gypsum, albite, illite, and k-feldspar. Then,…
-
Speciation with brine using Pitzer
I am trying to run a simple speciation problem with a brine composition in contact with some minerals (like Albite, Barite, Illite, K-feldspar etc.). The idea is to use the results of the speciation…
-
chemistry did not converge at node
Hi，tough community users, I try to carry out the 1D fluid-rock modeling for meteoric water dissolution during 10e4 years. In the runlog file, I occasionally found this error,…
-
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…