Christine Doughty's Profile - TOUGH User Community

Skip to main content

Latest Posts

Re: T2Well (Error message)
Dear Refaat, According to the message, your grid has 243220 connections but the code is set to have maximum of 57000 connections. In the T2 file, you need to find the parameter MNCON at set it to…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: TOUGH3; long running time (3D model)
Hi Refaat, I tried your input file, using TOUGH3 with ECO2N and got the same result as you. I tried decreasing the permeability in the monitoring-well material to 10^-10 and it runs the 20 time…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Mobility concept in TOUGH (Parameter MOP(11) )
Hi Diego, Mobility generally refers to the product of absolute permeability and relative permeability. If you have different absolute permeabilities for sandstone and mudstone,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: Phase transitions with tough3-eos7r
Hi Mikey, Using your input file and tough3-eos7r running in an Ubuntu Window on a Windows 10 laptop, I get the attached output. The error messages at step 25 look like yours,…
Attachment
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: TOUGH3; Running time is too long
Dear Refaat, When I run your problem on a Linux machine, I think I am getting results like yours - a small time step of about 10 s that does not grow. When this happens,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: TOUGH3; Running time is too long
Yingqi, I also noticed this problem (on Linux, but not on Windows) - the code wouldn't even start to run and gave me an error message about bad dimensions.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: TOUGH3; Running time is too long
Refaat. Could you explain what you mean by "instability of the numerical results"? I do not have specific suggestions for solver options. See the TOUGH3 manual for MOP(21) and the SOLVR block for…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: TOUGH3; Running time is too long
Dear Refaat, I do not see anything obviously wrong with your input file. When I try to run TOUGH3/eco2n on a Windows 10 machine, using an Ubuntu window,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Temperature limit
TMVOC has an upper limit of T=350 C and P=1e8 Pa. It checks at the start of each time step whether T and P are within bounds, and it stops if not. You should get a message "CANNOT FIND PARAMETERS AT…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TMVOC
- Reported - view
Re: Why does the partial pressure of carbon dioxide change so much when SG is greater than zero？
@Liheng Wang You need to run your simulation longer to get to steady state. At steady state, Pco2 in the middle element is the same as Pco2 in the high-pressure boundary,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: Why does the partial pressure of carbon dioxide change so much when SG is greater than zero？
@Liheng Wang It looks like everything I need is in flow.inp. Thanks.
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: Why does the partial pressure of carbon dioxide change so much when SG is greater than zero？
Very interesting problem, and nice illustrations. I have no experience with EOS2, but I think what is happening is that for the T=99 case when you have single-phase liquid conditions for all times,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: Steady-state simulation of a low-type geothermal system
Dear Hakim, You should use the Command Line option in Petrasim, where you specify the name of the input file you want to use (in your case, the T1_modified file).…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: How to improve the stability of the numerical results without changing the grid discretization ?
I can think of one easy thing for you to try: increase rock compressibility COM (second line in ROCKS entry for each material). Your current value of 3.71E-10 is pretty small.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Format of INCON and SAVE file
Yep, I think this is one of the most confusing aspects of preparing TOUGH3 input. Please read through the TOUGH3 User’s Guide carefully. There are several ways to get anisotropy and heterogeneity in…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH3
- Reported - view
Re: Format of INCON and SAVE file
TOUGH3 implements both the TOUGH2 method and the iTOUGH2 method for element-by-element permeability. If the material SEED is present in the ROCKS block, then the TOUGH2 method is used.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Cross flow between stacked layers
I don't know of a direct way to do this in TOUGH3, but you could use COFT on a series of connections between the stacked layers and add up the mass flow rates with a utility program or Python script…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Technical question regarding mass conservation in TOUGH3
Dear Refaat, Very nice plot! I have run the radial CO2 injection problem many times, but I never thought to plot the results that way. I just reran it and I see the same effect as you - liquid mass…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH3
- Reported - view
Re: Specifying an injection well with a constant volumetric rate
I checked the TOUGH3 User's Guide and source code, and I did not find a constant volumetric rate option for GENER. I think some other versions of TOUGH do have this option - maybe iTOUGH2?
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH3
- Reported - view
Re: General question about Capillary Pressure
In general, fluid flow is driven by pressure gradients, gravity, and capillarity. I think it is most important to include capillary pressure when pressure-driven flow is small.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: How to deactivate the mutual solubility of CO2 - Brine mixture in TOUGH3/ECO2N
I have done a similar study to investigate how different processes affect the temperature profile. This link https://drive.google.com/file/d/1Ryrih0MNjvinYUukg1GidvBxlz77Ys2o/view?…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 2
- TOUGH3
- Reported - view
Re: How to deactivate the mutual solubility of CO2 - Brine mixture in TOUGH3/ECO2N
I do not know of a way to deactivate the mutual solubility between CO2 and brine in TOUGH3/ECO2N. I have used the same approach as you have to minimize dissolution.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: How to identify the cause of an incomplete simulation run
In the ICP=7 formulation for capilliary pressure, Pcmax < 1.0 acts a flag that changes how the Pc versus Sl curve is calculated for small values of Sl.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: How to identify the cause of an incomplete simulation run
@Refaat G Hashish You asked for the code to take 9999 time steps (PARAM block, first line, columns 5-8) and the code took 9999 time steps and stopped. It is doing just what you asked.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH3
- Reported - view
Re: How to identify the cause of an incomplete simulation run
I don't see any obvious problems with your input file. Could you post your output file?
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view

1 2 3 4 5 6

Back to top

Home

View all topics

TOUGH Website