-
Difference between 'e' and 'd' in power expression
Hi, it might be a minor problem though (because I am a newbie to scripting), Sometimes the power expressions are written in 'd', not 'e' (e.g. 1.0d-9, not 1.0e-9 especially in DIFUN filed).…
-
Documentation on the Numerical Aspects in Forming and Solving the Equations
Hello, There are many papers on Toughreact but I am searching for a documentation/paper on forming the equations i.e. accumulation, convection, reaction terms and solving these equations.…
-
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.…
-
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…
-
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.…
-
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…
-
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.…
-
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…
-
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,…
-
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,…
-
TOUGHREACT Workshop Files
Hi, How can we get the files which are mentioned in the TOUGHREACT workshop handout here: https://tough.lbl.gov/assets/files/02/support/TOUGHREACT_Sample_Problems.…
-
CONSTANT pH PROBLEM
Hello everyone. After running a simulation (injection of CO2 into a saline aquifer), I observed that in my result, the pH remained constant at all times for my chemical initialization, injection,…
-
Batch Model with Brine Injection and Secondary Minerals Precipitation
I am attempting to set up a batch reaction model. Initially, the formation water equilibrates with calcite, gypsum, albite, illite, and k-feldspar. Then,…
-
chemistry did not converge at node
Hi,tough community users, I try to carry out the 1D fluid-rock modeling for meteoric water dissolution during 10e4 years. In the runlog file, I occasionally found this error,…
-
how to set initial temperature distribution
Hi,tough community users, I want to prepare the initial or stable pressure and temperature conditions for EGS. INDOM was introduced into my flow file. The UPPE1 reservoir is the top mesh having 1.…
-
About the fluid velocity in the governing equation
Hi, everyone, I'm currently studying the governing equations underlying TOUGHREACT, particularly in the context of the ECO2N module. The entire solution process is as follows: (1)For the fluid flow…
-
total concentration of sulfate increasing over time
Hello, I ran a chemical initialization simulation for 10 million years to achieve steady state. In order to figure out how much dissolution and precipitation I have in my system,…
-
low injection rate cannot work in 2D model
Hi, everyone. I am so sorry to interrupt you again for a same question. I want to simulate the dissolution and precipitation of sandstone in 2D model, which were induced by meteoric water injection.…
-
V1.2 TOUGHREACT ECO2N IN PetraSim:Problems with initialization for flow and chemical field
Hi, I am new to TOUGHREACT and using v 1.2 ECO2N on PETRASIM to model CO2 Sequestration. And I have created a 2-D isothermal rectangle model and designed a temperature gradient of 3.…
-
TOUGHREACT - Dissolution and precipitation
Hi all, I am a new user of Toughreact. I am testing a real case dataset, whereas I want to quantify the amount of CO2 that will be mineralized in a mafic formation.…
-
Therakin.dat question
Hello, I would like to know what some of the numbers in therakin are. 'dolomite-2' 184.401 64.365 4 -2.0000 'h+' 1.0000 'ca+2' 1.0000 'mg+2' 2.0000 'hco3-' 'dolomite-2' 6.9420 5.7530 4.2300 2.6830 0.…
-
About the toughreactEco2n v1.2 not have the “savechem”files
Hello everyone, I was using TOUGHREACT V1.2 to run a initial progress of HC field(no thermal field). There were some problems when adding the chemical field,…