-
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…
-
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...…
-
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:…
-
PHYSICAL MEANING OF INPUT RECORDS
Hello, I am using Toughreact version 4.13 to solve CO2 storage problems. However, I still don't understand the physical meaning of PRIMARY SPECIES and AQUEOUS COMPLEXES,…
-
Multiple questions regarding data entry for a surface complexation problem
Hi there, I am working on a coal combustion product (CCP) leachate project. The goal is to understand the natural attenuation and reversibility for some remediation drivers at CCP sites,…
-
Question on modeling water isotopes – after Singleton et al. (2004)
Hi Eric, I have been playing around with your water isotopes model of Singleton et al. (2004) and have a question about how fractionation factors are implemented in TOUGHREACT.…
-
Convergence problem chemistry
Hello, everyone, I was running a water-rock interaction model with Toughreact. The runlog.out file prompted "Warning: chemistry did not converge at node A1114 (routine NEWTONEQ),…
-
Species appear two times in time.out
Hi, I have a problem in reading result file. Some primary species in <time.out> appear two times; one column with the header of "t_xxx" and the other start with just "xxx" without "t_" ahead.…
-
Is there a limitation in the number of grid blocks to be printed out?
Hi- I ran two simulations with different number of grid blocks to be printed out in 'time.out' The first case with 10 grid blocks had no problem at all but the other one with 100 blocks encountered…
-
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…
-
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,…
-
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…
-
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,…
-
Volume and Mass Balance Error in flow.out
Hello, everyone. I am using EWASG model to simulate chemical reaction (Ba+2 + SO4-2 = BaSO4) in porous media. Although I define the initial water and assign it to the matrix,…
-
TOUGHREACT_Debugging
Hello Everyone: I am pretty new to TOUGHREACT. I am trying to build a model to simulate Ba++ and So4-- reaction in porous media. The attachments are the input files I modified based on P8_EOS1 from…
-
A question about the additional mechanism under the Minerals entry in the chemistry.inp file
Dear All, Hello everyone! I just started learning toughreact software. At the moment, I am a bit confused about the description of the additional mechanism under the Minerals entry in the chemistry.…
-
Two time steps with restart gives different results than one without restart
Hello, Do you know why a restart during 6 hours (for example) after one time step (6 hours for example) is different from one time step of 12 hours? In my simulation with eos4,…
-
What is the use of 'aquatic complexes'
I am wondering how aquatic complexes may change the simulation results. below are my findings: 1. aquatic complexes, if declared in 'chemical.inp', are assumed to be at equilibrium with basic…
-
Compilation TOUGHREACT on Windows using gfortran and Cygwin
I am having some problems with compiling the original source code of TOUGHREACT. I have tried using all of the compiler flags and the ones suggested on this forum.…
-
Failure restarting Toughreact simulation.
Hi Everyone, I want to run a model with the geochemical data of a previous model. I follow the manual instructions changing the names of the SAVE and savechem files to INCON and inchem.…
-
Failure to set up a constant boundary concentration by input a large volume in the first cell
Dear Sir or Madam, Sorry to disturb you. I set up a simple 1-D contaminant transport model using TOUGHREACT v3.0 following the 8.1 Problem 1 – Aqueous Transport with Adsorption (Linear Kd) and Decay…
-
Geochemical Equilibrium between 2 water compositions
Hello: I understand how to obtain the speciation of one water composition in one grid block for a given thermodynamic conditions. But I want to know if it's possible to obtain the geochemical…
-
TOUGHREACT ECO2N
Hi, I'm quite new in using TOUGHREACT. Can anyone please clarify one question I've been having for quite a while now? The system I'm trying to model is just simple mixing between brine defined by the…
-
Cooling cell
Hello and thanks for maintaining this forum, very useful indeed. I am running a variation of the 5 spot problem in Toughreact. Essentially this is a polygonal mesh and there is a 10 meter high…
-
Mass Balance for trace gases in TOUGHREACT V3-OMP
Dear colleagues, I stumbled about a little problem concerning the mass balance for trace gases using the new TOUGHREACT V3-OMP version. I am simulating the injection of CO2 containing trace gases…