-
How to handle time step control for proper chemistry computation?
Dear all, I am struggling with some Toughreact basics and maybe some of you can explain me shortly how to handle the time step control in a proper way to avoid mistakes and miscalculations.…
-
restart with surface complexation
Hello, I have been trying to run a restart with surface complexation. First, I get an error like that: "Restart time in INCHEM file is not the same as starting time for flow calculations. timetot= 3.…
-
zero flux for ions at boundary condition
Hello, When simulating a drying process of a porous media it is possible in TOUGH2 to apply a capillary pressure at the boundary condition. However, when doing such a simulation with a porous media…
-
'nkdd' field in TOUGHREACT input file
Hi- I participated 2014 TOUGH(REACT) short course and new to here. I am using PetraSim only to generate basic input files and modify them to run on our university's cluster computer (PetraSim…
-
Changing boundary condition for chemistry in EOS3, TOUGHREACT
Dear all, I would appreciate it very much any help with the following: In the problem I model I use constant concentration chemical boundary condition,…
-
Restart problem
Hello, I am running a series of simulations which were running well until the last restart. With the most recent restart, when I wish to restart the program immediately stops and gives this error…
-
Convert TOUGHREACT output to Tecplot readable output
Has someone generated a utility to convert TOUGHREACT output files into files that are readable by TECPLOT? I found the EXT utility for converting TOUGH2 output into TECPLOT-formatted files,…
-
Boundary water
Hi all, I'm modeling CO2 underground injection and we are injection a mixture of brine and supercritical CO2. Under the system pressure and temperature and the injection rates,…
-
version of shared libraries
I just tested the Linux executable of TOUGHREACTv3 on our Linux cluster and I got the following error message before it crushed: /home/ycda/bin/treactv3omp_eos1_linux_intel: /vend/intel/intel_cc_10.…
-
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…