-
Building regions in PyTough
How can I build regions using pytough? I thought regions may be created by assigning rock properties to block centers, but it seems like that's not the way to build regions.…
-
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…
-
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.…
-
How to choose the module to simulate the Terzaghi's one-dimensional consolidation?
Hi everybody, as stated in the title，i want to simulate Terzaghi's one-dimensional consolidation, but i'm not sure which module of TOUGH2 i should choose, can any body tell me? thanks very much.…
-
Gas density in EOS2 is unrealistic(5000kg/m3)
Hello all: I set up a model to simulate CO2 displace water in core by EOS2 module. The inlet is CO2 mass inlet boundary condition, and the outlet is constant pressure boundary condition(10MPa).…
-
Temperatures reaching out of range in EOS2 despite no GENER block
Hi all, I am trying to run a model operating at temperatures close to the maximum modelled by TOUGH (360C) due to the system it is based on. Since this is only 13 C cooler than what is achievable by…
-
slightly increased liquid viscosity
Hi, dear all I have one question about eos1. To construct a sensitivity matrix 'co' in subroutine 'multi' , we obtain secondary parameters using slightly increased primary variables,…
-
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…
-
To have "Component mass in place" in a single cell
Hello everyone , In the output file of TOUGH2, for each plot time value, there is the -component mass in place- table, in which I think, are reported the mass of each component referred to the total…
-
How to calculate the mechanical dispersion in TOUGHREACT
Hello. I am inquiring about how to know the mechanical dispersion in the PetraSim/TOUGH. When I input the diffusion coefficient and tortuosity/porosity,…
-
The CO2 plume migraion is very strange.
The CO2 plume migraion is very strange. Can anybody tell me why is this case? Thanks very much. Sawyer.
-
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 determine the influence of boundary size?
Hi all, when i use tough2(ECO2N), i'm not sure whether the boundary size i chosen is appropriate, can anybody tell me how to determine the boundary size? thans very much! sawyer
-
How to set up source water volume as boundary condition
I have one question regarding how to set up the water source volume in petrasim/tough2. For example, the reservoir is bottom water type. There is water driven at the bottom of the reservoir,…
-
Mineral volume fraction greater than 1 in outsolid.dat
Hello, I have a problem with TOUGHREACT PITZER. I'm running a model with calcite, anhydrite and dolomite whose initial volume fractions are 0.99, 0 and 0.01 respectively.…
-
Coupling between flow and processes
Hi there, I am wondering how the fluid flow and heat transport are coupled in Tough2. Are fluid and heat flow equations solved iteratively ? I would be thankful if could help me in this regard.…
-
what is the meaning of the running time is negative?
Hello everyone, can anybody tell me why the running time is negative? As shown in the figure? Thanks very much! sawyer
-
Debugging
Hello I am running a problem and the simulation was halted at the middle of simulation. Here is the text that I found in the output file. How can I figure out and fix the problem?…
-
How to handle time step control for proper chemistry computation?
Dear all, I am struggling with some Toughreact basics and maybe some of you can explain me shortly how to handle the time step control in a proper way to avoid mistakes and miscalculations.…
-
TOUGH 2 ECO2M temperature condition
Hi, It is described in the ECO2M user guide that TOUGH2 simulates CO2 migration under temperature condition up to 110 °C. I am wondering if the temperature condition is about 115 °C-125 °C,…
-
production mass volume calculation
Hi there, For a constant boundary pressure, and constant production rate gas reservoir, is there a way to tell how much gas and how much water produced from outputs? Thanks Jing
-
VOL. GENER TYPE
Hi everyone, Wandering in the t2f.f file, I found that there is a GENER type not referenced in the TMVOC manual : "VOL." I assume it allows specifying an injection or pumping volumetric rate.…
-
Problem with RF1 input file
Hello everyone, I'm starting to use TOUGH2 for the first time. I tried to run the sample-input file RF1, but the output file gives me an error and I don't know how to proceed.…
-
'nkdd' field in TOUGHREACT input file
Hi- I participated 2014 TOUGH(REACT) short course and new to here. I am using PetraSim only to generate basic input files and modify them to run on our university's cluster computer (PetraSim…
-
Dirichlet boundary
Hi all I am trying to make a boundary with pressure variations. I have two questions: 1- Should I just follow the equation 27 in TOUGH2 manual (ideal gas law)?…