TOUGHREACT Category

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…

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…

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).…

Error in Problem 6 2-‐D Radial Model for CO2 Sequestration in Deep Saline Formation from TOUGHREACT V4

Hi everyone, I have been trying to run TOUGHREACT V4 Deep Salt Formation 2-D Radial Model Problem 6 with the instructions in the manual. however, I have only obtained results for the first 10 years…

TOUGHREACT error while loading shared libraries: libifport.so.5:, but used gfortran complied

I was trying to run some test files using treactv413omp_eco2n_linux, but I ran into the error: error while loading shared libraries: libifport.so.5: cannot open shared object file:…

Adsorbed amounts included in concentrations?

Could you please confirm whether the total concentrations that are output from a TR simulation include the adsorbed amount when using Kd linear adsorption?

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.…

Problem 5 – 1-‐D Radial Model for CO2 Sequestration in a Deep Saline Formation (ECO2N)

Hello.. How do you get the results to graph the figure Figure 1.5—3 Cumulative CO2 sequestration by carbonate mineral precipitation for different times (a) CO2 sequestrated?

Output flux of aqueous species

I would like to output the flux of all aqueous species across specified connections. Is this possible? I am using TR3.32. Thanks, Heather

LINEQ

I am trying to model chemistry reactions in porous media. To build the model from the base, I run the case without injection/production or any chemistry reaction. From the flow.out, the MAX.…

CO2 injection in a depleted oil field

Hello, what module do you recommend for the simulation of CO2 injection in a depleted oil field? Some examples would also be helpful.

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…

Toughreact V4.13 COFT output

I have problems with the COFT output. I actually doen't know how to read it. I'm using Toughreact 4.13 Executeable, and the output of the COFT file give only some of the Indexes,…

inchem and INCON in TR- difference with TMVOC

Guys Just realized that if you remove +++ at the end of INCON, TR does not read inchem file. This is different with e.g TMVOC Kind regards Kaveh

Reaction rate >0 at equilibrium

I am running a simple test problem representing dissolution of a hypothetical "glass" mineral in a single grid cell. The solid phase consists entirely of glass, the pore fluid is water.…

Recompiling with higher number of ELEME/CONNE

Hi I remember I read somewhere which parameters to be modified to recompile the code (omp-v3.3) to include higher number of elements and connections. However I cant find it.…

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.…

Radioactive chain decay

Hi Is it possible to simulate radioactive chain decay (i.e parents produce daughter elements) in TOUGHREACT? It is more complicated than the simple decay equation,…

Mass balance issue on mineral phase during CO2 injection in radial system

Hi everyone, we use a typical case of CO2 injection (see details in Andre et al. 10.1016/j.enconman.2007.01.006) at 200 b in a radial system, so it was using the 2007 version of toughreact.…

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…

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.…

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.…

CO2 solubility ECO2n carbonated water

Hi everyone, I'm working on CO2 mineralization in the subsurface where carbonated water is injected into reactive minerals, similar to the CarbFix project.…