-
Compiling TOUGHREACT V3.0-OMP using Cygwin and gfortran
Hey everyone, I have so far been able to compile all EOS modules of TOUGHREACT V3.0-OMP using gfortran and the Cygwin Windows utility. This was done using the following options:…
-
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.…
-
Pre- & Post- processors - What are their Pros and Cons?
Greetings I need to develop reactive-transport models with TOUGHREACT and run simulations as part of a project on "subsurface barrier formation as a CO2 leakage mitigation and remediation…
-
Tuning diffusion coefficients for water isotopes
Hi @Eric Sonnenthal I have a few questions about modeling diffusion for water isotopologues in an evaporating soil. Basically, the setup (following Singleton et al.…
-
Compiling and Running TOUGHREACT v3.32 using Ubuntu/Linux
Hey everyone, I've been able to get the newest version of TOUGHREACT working with my workstation, which has Linux/Ubuntu (v 17.04) as its OS. Setting things up didn't go without hitches, though,…
-
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?…
-
Pressure drop using GENER
Dear Tough community, I'm trying to run a 1D simulation with initial and boundary waters. I'm using GENER to drive the flow. The program runs, but I get a drop in pressure for each step in flowdata.…
-
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.…
-
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.…
-
TOUGHREACT compilation using gfortran MinGW-error message
Dear Community, I would like to compile TOUGHREACT (version: TR-OMP_Source_v3.3.2-OMP) using gfortran with MinGW. I typed make eos1 -f makefile_gf, and had the attached error message.…
-
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,…
-
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.…
-
set different diffusion coefficient for different layer
As far as I know, only one solute diffusion coefficient can be set for all aqueous species,which is not practical. Is it possible to set different diffusion coefficient of one solute for different…
-
No convergence during CO2 injection into heterogeneous mineralogy
Hello, I am very stuck on a reactive transport modeling project, and was wondering if someone would be willing and able to help me. I have initialized a homogenous (100% quartz) 2-D aquifer system…
-
Setting up initial water chemistry
Hello, I am trying to set up a TOUGHREACT model of reactive transport during upwelling of CO2 and salt, starting by emulating the chemistry presented by Zhang et al., 2016 (attached).…
-
Lithofacies random fields in TOUGHREACT?
Hi, I am trying to apply random fields both in flow and rock chemistry. I know I can do one for the flow field using permeability modifier (PMX) option equipped in TOUGH2 (flow.…
-
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…
-
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.…
-
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,…
-
Error in reading INCON - Problem 5, eco2n_1D Radial
Hello, In problem 5, eco2n_1D Radial, I am using the input files provided (flow, solute, chemical, CO2TAB, .dat, and executable) to run a q=0 kg/s flow in order to find the equilibrium conditions…
-
Using TOUGHREACT with zero porosity boundary
Hey everyone, I have recently tried to use a zero porosity boundary as I run a TOUGHREACT 2.0. In TOUGH2 manual, I have found out that this is allowable.…
-
Atmospheric carbonation of concrete with EOS4
Hi all, We are trying to run a case simulating atmospheric carbonation of concrete with Toughreact V3/EOS4. Running the case without reactivity to check the transport of a tracer works fine (see…
-
Confused by the choice of boundaries based on the experiment
Dear Sir or Madam, Sorry to disturb you. I want to simulate the 1D contaminant transport using TOUGHTRACT v3.0. My experimental model is the simple constant injection with constant concentration C0…
-
Initial Geochemical Equilibrium between Minerals and Brine
Hello there, I'm simulating CO2 injection into sandstone aquifer of synthetic brine, NaCl. I've read in several papers in which TOUGHREACT w/ ECO2N is used and something is mentioned like "batch…