-
Postprocessing ToughReact (or Tough2)
Hi. Any preferred tool for post processing, including visualization and analysis ability for ToughReact? Seems that Tecplot is used by some people. But the vendor states that the output files (or…
-
Difference between surface area for sorption and initial mineral reactive surface area
Hi, everyone, I am confused by the concepts of S_AREA (surface area for sorption) and AMIN (initial mineral reactive surface area) If I choose IFLAG=0 and IMFLG=0 and only one mineral for the…
-
Setting Boundary Conditions in Toughreact
Hello All I know how to impose Dirichlet and Neumann conditions for a single column volume like A in my sketch (one dimensional case). In the image, the blue cell represents a injection block.…
-
Convergence problem with different mesh size
Dear Sir or Madam, Sorry to disturb you. I was running a simple 1D transport model with surface complexation, but it never ends. However, I found it worked when I increased the mesh size from ∆x =1.…
-
How we can define a reaction with two minerals (such as reaction 1) in Toughreact’s thermodynamic database?
I am trying to simulate the precipitation of uraninite using reactive mass transport modeling. In my reactive transport model, Uraninite precipitates through the following reaction:…
-
Initializing simulation with sorbed species
First, this is my first foray into building up my own TOUGHREACT simulations, so I apologize if this is a question I should have figured out by reading the manual more closely.…
-
how to stop overwriting output file of concentration everytime
Hey everyone, I set to monitor chemical concentrations for two specific elements in solute.inp. Like this: Writing control variables [7i5]: 1 2 2 4 0 1 1 ! NWTI,NWNOD,NWCOM,NWMIN,IWCOMT,iconflag(=1:…
-
How to implement a linear thermal gradient?
Greetings: I'm trying to impose a linear thermal gradient over time, i have seen elsewhere that one way to kept the temperature constant is by specifying high densities (1E15 kg/m3) in all the…
-
Cation exchange and surface complex input files
Hello everyone, I'm simulating the cation exchange and surface complex between Eu+3 and clay, and I found only the 8.2 problem 2 and 8.10 problem 10 in TOUGHREACT v2.…
-
Does the ExtraCell function run well
I want to simulate a small flooding, the core is tiny and I want to add an ExtraCell, where can I find the manual about the meanin of the parameters?
-
The reactive portion of the model seems not to run
Hi! I'm a new user of this interesting forum. I'm trying to simulate a diagenetic process by using PetraSim 5.4 (TOUGHREACT and EOS3). Therefore, for each new simulation,…
-
salt, chemistry and permeability reduction in ECO2N
I am trying to simulate salt precipitation and corresponding permeability reduction upon CO2 injection using the ECO2N EOS in TOUGHREACT. Initially without chemistry.…
-
Problem on the number of grid blocks for time evolution > 15
Dear Sir or Madam, I am sorry to disturb you. I wanted to print out 60 grid blocks for time evolution printout. And I just listed all the grid numbers in several lines (15 in every line) in the…
-
no runlog.out file for error
Hi everyone, I am a beginner using toughreact. Currently, I am building up a 3 layered single phase model. After I built up everything, the model stop running after 1.8 seconds and shows me error.…
-
The simulation of the general three types of boundaries
Dear Sir or Madam, I just want to simulate the general three types of boundaries (without any specific cases) in the inlet end: Dirichlet, Neumann, and Robin.…
-
Simulate the cation exchange and surface complex between heavy metal and clay
Hello, everyone. I am trying to simulate the cation exchange and surface complex in clay. However, my simulation results are very strange. Could someone please help me figure out how to how to…
-
Error compiling the utility kswitch_v2.3.f
I'm attempting to compile the kswitch_v2.3.f source code that came with TOUGHREACT v3.32. I'm using the gfortran compiler, and if this is relevant, I'm running this on Ubuntu v 17.10 as my OS.…
-
Modelling the incongruent dissolution
Hi all, I'm trying to simulate CSH(Ca-Si-H2O) dissolution, which is usually considered as an incongruent solubility behavior. I'm thinking how to simulate the reaction of CSH with high Ca/Si ratio…
-
Convergence problem for the surface complexation model
Dear Sir or Madam, Sorry to disturb you. I am Lang. I created a 1D transport model with surface complexation. I found it could not achieve convergence,…
-
Time stepping with 'infinite volume' boundaries
Hello, I was wondering if you all could help me understand the time stepping of a soil model I’m building. I would like to model annual soil carbonate formation and dissolution dynamics in a soil.…
-
CEC: cation exchange capacity (meq/100 g of solid).
In ToughReact, does the cation exchange capacity have the time vector?Does it represents meq/100 g of solid/s?
-
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…
-
Cells with time-dependent concentration?
Hi, I found that time-dependant boundary conditions in flow are possible in TOUGH2 by using FOFT and GOFT blocks, but is it also possible to have a time-dependent concentration at certain element?…
-
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.…
-
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).…