-
Is TOUGHREACT only suitable for NaCl salt?
Hello everyone, I am conducting an experiment using different salts and distilled water to prepare brine. I was wondering if I can simulate the experiment using TOUGHREACT as it appears Toughreact…
-
Use input files for injection rates at many timesteps?
Hello, I am trying to include many CO2 and brine injection rates over a range of timesteps in TOUGHREACT. I was wondering if I would be able to call a file with all of the data I want in GENER,…
-
PHYSICAL MEANING OF INPUT RECORDS
Hello, I am using Toughreact version 4.13 to solve CO2 storage problems. However, I still don't understand the physical meaning of PRIMARY SPECIES and AQUEOUS COMPLEXES,…
-
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…
-
ELEME & CONNE blocks
Hello everyone, I am trying to understand how to calculate grid blocks volumes and interface areas of heat exchange in the ELEME block of Toughreact flow.inp file.…
-
How to confirm mole fractions of injected CO2 and SO2
Hello, I'm new to Toughreact codes. I'm simulating 1-D flow through a rock sample of dimension 50mmx100mm. CO2 is injected with SO2 impurity (10% of CO2 volume fraction) for 24 hours (1 day) and the…
-
cguess and ctot in chemical input
Dear all, I would like to ask what does cguess and ctot mean in the chemical input in ToughReact? I have read the user guide, I have also searched in papers and still haven't found an answer what…
-
How to calculate what portion of a constituent adsorbs.
Hello there, I ran a 1-cell chemical initialization simulation that includes surface complexation for 1000 years. For some constituents such as Cr, the total sorbed moles are more than their initial…
-
changes in chemical.inp
Dear all, I am a totally new Toughreact user so please forgive me if my question is a bit stupid. I want to modify the chemical.inp from one of the Toughreact examples (1d-CO2).…
-
Stopping in main due to flow nonconvergence
Hello everyone, I am using ECO2N to simulate CO2 injection into a formation with relatively low permeability. I am having trouble with an error popping up.…
-
OUTPT section
I use cubes to construct an octahedron prism. The cross-section shape is shown in the attachment. The largest distance in three directions is 41, 41, 291. How do I write the OUTPT section? Thanks!
-
Permeability-Porosity relation
I use the TOUGHREACT EOS1 to model chemical reactions in porous media. In the Permeability-Porosity relation, I choose the Verma-Pruess correlation. The equation is k/ki=(phi-phi_c/phi_i-phi_c)^n.…
-
Multiple questions regarding data entry for a surface complexation problem
Hi there, I am working on a coal combustion product (CCP) leachate project. The goal is to understand the natural attenuation and reversibility for some remediation drivers at CCP sites,…
-
How to monitor the temperature change of the sink with time?
Hi gays, I'm trying to simulate plume geothermal systems using toughreact. I want to output the temperature change of the sink at each time step, but I don't know how to set it.…
-
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.…