TOUGHREACT Category

Problems precipitating phases at high ionic strengh using trpitz_v1.21

Hello everybody, I am working on a geothermal application and I try to model halite precipitation during the cooling of a brine in a borehole. Even when the brine is oversaturated by far with…

Timestep problems while including O2(g) in calculations (V3-OMP / ECO2N)

Hello everybody, I have found the next problem for you to resolve :-). I am using Toughreact V3.0-OMP, but the problem would appear with V2, too, I assume.…

Prescribed boundary conditions for chemistry but constant flow rate BCs for flow

Hi all, I'm using TOUGHREACT with ECO2N module. In the simulations, the flow is injected at a constant flow rate and due to dissolution/precipitation, the pressure would change at the boundary cell.…

Role of ico2h2o and mass balance

I have run a number of simulations using Toughreact Versions 1.2 and 3OMP and the ECO2N EOS. In trying to understand the mass balance I have tried runs with 0 entered for ico2h20 and 1 or 2 as well.…

Questions on aqueous kinetics, Sample Problem 10.

Hello everybody, I have some questions concerning the representation of aqueous kinetics reactions and especially the sample problem 10 from the Toughreact V3-OMP description. 1) In table 1.10-2 (p.…

Recommendations on using TOUGHREACT to perfome reactive transport modeling of hydrogen in a depleted gas reservoir

Dear All, I am dealing with the problem of reactive transport modeling of hydrogen in a depleted gas reservoir. I was wondering if anyone has already done reactive transport modelling of H2 with…

ECO2N - possible to solve for temperature in sample problem?

Hello All, I have been working with the (*rcc3*) example from the ECO2N user manual: Radial Flow from a CO2 Injection Well. This problem has a constant temperature of 45C.…

Link error - ToughreactV3.0

Hi! I have the following error message when compiling TOUGHREACTV3.0 using Microsoft Visual Studio 2013: fatal error LNK1248: image size (F3F2E000) exceeds maximum allowable size (80000000) Does…

Please help with 2D version of 1D radial flow ECO2N example

Hello Everyone, I am very new to TOUGHREACT, and in all honesty, I know I am in over my head, but I refuse to give up. I apologize in advance, although I find my problem very difficult for me it…

Input conversion error, unit 33

Hi there, I've modified a sample input file flow.inp (in P5-ECO2N) to use our dimensions and parameter values. But after running the toughreact executable, this is the error message I've got.…

initial volume fraction of the mineral

Hello everybody, I have just one trouble, could "initial volume fraction of the mineral" be superieur 1 ? Thank a lot.

TR-V3-OMP: Composition of gas zone / VOLG

Hello everybody, I am using the new Toughreact Version 3.0 OMP with ECO2N EOS. I took the Sample problem no. 5 (1D radial model for CO2 Sequestration) and tried to add some SO2 as tracer gas in the…

How to assign pore pressure gradient

Hi, I have a question in assigning the pore pressure gradient in TOUGHREACT? I tried to assign the pore pressure to each zone by the txt file, but failed. Is there anyone has the experience? Thanks.…

TOUGHREACT ECO2N, temperatures higher than 110 ˚C

I am using Petrasim 5.3.1203 for a scenario to inject water and 5% CO2 in a deep geothermal reservoir to examine the mineral interactions and their impact on the fluid flow.…

Question about Restart in TOUGHREACT

Hello Everybody, I am struggling to get running a restart and I would really much appreciate your insights. I'm using eos1 in toughreact_V2. The first run is a 2d model in which I have 8 minerals,…

Failed mass balance for CO2 if "Two Fluid Phases conditions" employed

Hello everyone, I am trying to run the geochemical equilibration of brine with a sandstone rock composed of several minerals (mostly quartz 70%, around 4% carbonates and some feldspars and clay…

How to set-up a model that reacts only on the surface

I am trying to set-up a model that simulates a leaching experiment. The material is very impermeable (flow is not supposed to go through) and based on existed data,…

Difference between TOUGHREACT 1.2 and 2

Hi all, I am observing quite different (quantitative) results between 0D and 1D simulations run with the 1.2 version of TOUGHREACT/ECO2N and the 2.0 version.…

Different output results for different max time step (DELTMX) values

~~I am running a Toughreact model. I run this model with the time step of 1200 seconds (DELTEN=1200 s) using different max time step values (DELTMX= Infinite, 15 years,…

Rapidly reducing temperature in cells and MINC problem

Hello, I'm trying to recreate the simulation of the Coso EGS as published by Kovac et al on the Stanford Geothermal workshop in 2006. I have two questions that I would be grateful to have some…

Question about constant rate

Hi, everyone, I am trying to add the constant mass injection case in my model. Based on the manual, I uses the "mass" to represent the injector by using the dummy element.…

what if a primary species is not defined in the database?

The modeling I am working on involves some initial species that are not defined as primary species in the database. Although I have figured out how to define a new mineral in the database,…

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?

Restart　timestep

When restarting a run after long time calculation (e.g.1000year) and continuing the run until 10000year, I can't start the restart run properly. In the SAVE file,…

Why it takes unbelievable long time to run kinetics?

It seems that TOUGHREACT starts to choke when I add kinetics to the chemistry of the my model. The model runs just fine for equilibrium. But once I add kinetics,…