-
Restarting TOUGHREACT Simulations
I'm trying to restart TOUGHREACT simulations by renaming SAVE to INCON and savechem to inchem. In between simulations, I'm adjusting the initial waters associated with certain elements by changing…
-
Error when trying to output the time series
So, I have 40000+ elments, and I want to output all of them's data, so I put all the elements in the solute.inp, for example But then I have the error Elements for which to output time series:…
-
Can TOUGHREACT automatically calculate injection rate into each element?
Hello, I was wondering if TOUGHREACT can automatically calculate the injection rate into each element without the user interfering? For example, I would like to injection 90 kg/s into a layer that is…
-
CONVP block in solute.inp
Hello everyone: I have a question about implementing CONVP (chemical convergence bounds for time stepping) in solute.inp. I firstly make MOPR(10)=2 to make simulator read CONVP block.…
-
Restart time in inchem file is not the same as starting time in INCON - 1D Radial CO2 injection - Problem 5
Hello, I removed the GENER block from the flow.inp file in order to run an equilibration simulation for 10,000 years. I then wanted to use the 10,000th year data to be used as my initial conditions…
-
Understanding SOLVR - From problem 5, 1-D Radial Injection of CO2 flow.inp file
Hello, I have been looking for what the 3rd, 4th, and 5th lines of SOLVR in the 1D Radial CO2 injection problem are but I cannot locate them in the manual.…
-
Associating the variables to the numbers in a flow.inp file
Hello everyone, I am trying to understand what is what in the ELEME and CONNE sections of the flow.inp file for "1-D radial flow problem for CO2 injection into a saline aquifer" (problem 5 ECO2N).…
-
Problem 5 –1-D Radial Model for CO2 Sequestration in a Deep Saline Formation (ECO2N
Hello everyone, I have been working with the example "Problem 5 –1-D Radial Model for CO2 Sequestration in a Deep Saline Formation (ECO2N)". This problem has a constant injection of CO2 for 10 years.…
-
Issue modeling a soil with zone of unsaturated vapor
Hi, I am using EOS3 to model a evaporating soil and have encountered some behavior with soil air relative humidity that I think is unrealistic. My 1D setup has a static block at the bottom to act as…
-
Running TOUGHREACT LBNL compiled on an AMD CPU
Hi, I am using an AMD Threadripper 24 core CPU and this cannot use kmp_affinity of OpenMP so you get an warning message #72 and #71. kmp_affinity only supported by Intel processors and kmp_affinity…
-
How can I run the example problems in UBUNTU 20.04.2?
Hello, I am new to using ToughReact (V3.32), and I have no idea how to run the executable already read the installation file but, I have not been able to solve the problem. Use ubuntu 20.04.…
-
TOUGHREACT (v1.2) not considering fluid density due to salinity?
Hello, I am using EOS9 to model a plume of leachate entering the top of an aquifer to evaluate the resulting changes to water & mineral chemistry. I have noticed that when injected,…
-
Problem running TOUGHREACT in Petrasim
Hi, I am new to TOUGHREACT. I set up a simulation with material properties and other parameters. Created MESH. However, TOUGHREACT Halts Simulation prematurely with an error message "MESHMAKER HAS…
-
Initialization failed in TFREACT Econ
Dear everyone, Currently I am learning to build a simple model using EOS 7 (3,4,3,6), incorporating salt and CO2 injection. I have been stucking in the initialization step for quite a few days.…
-
Reconciling Flow and Chemical Inputs
Hello, I have a couple questions that are very basic and therefore hopefully rather straightforward for someone to answer. I am just getting started using TOUGHREACT,…
-
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.…
-
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.…
-
Empty FOFT file
Hello, I was wondering if someone could help me with the following problem. I am running ToughReact1.2(eos3) and would like to use the FOFT keyword to obtain time-dependent data for selected…
-
Failure of flow.inp in TOUGHREACT v3.0 generated by PetraSim
Dear Sir or Madam, Sorry to disturb you. When I run the flow.inp generated by the PetraSim (TOUGHREACT v1.2) in TOUGHREACT v3.0, and I found it could not work.…
-
Effect of the diffusion coefficient(D0) on the TOUGHREACT
Dear Sir or Madam, Sorry to disturb you. I found a problem when I run the example model in the 8.1 Problem 1 – Aqueous Transport with Adsorption (Linear Kd) and Decay (EOS9) in TOUGHREACT User Guide…
-
Helium diffusion into CO2 plume
Hi, I am running toughreact V3 example problem 5 with diffusion switched on and helium added to the fluid and CO2 phase. During and after CO2 injection i expect the helium to diffuse into the CO2…
-
MOPR(1) = 4 (constant velocity)
Hi, I know most of simulation uses the first option of MOPR(1) as 0: reactive transport or 2: only flow with no chemical reaction but when do we use 4:…
-
Convergence problem in initialization of water composition
Hello, everyone, When I run my TOUGHREACT model, the runlog.out file displayed "ERROR (convergence problem in initialization of water composition) Please adjust convergence criteria regarding…
-
How to simulate the surface complex in PetraSim
Hello, everyone, Sorry to disturb you. I wanted to simulate the surface complex between aqueous Europium+3 and -S-sOH and -S-wOH (the defined group) in clay.…
-
How to define the tortuosity in TOUGHREACT
Hi. I am confused by the concept of the tortuosity. I found it was defined as is the ratio of the average length of the flow paths through the medium to the straight-line length across the medium…