Christine Doughty's Profile - TOUGH User Community

Skip to main content

Latest Posts

Re: Not to achineve the end time and no error
The code seems to be pretty unhappy, reducing time step frequently. To know more about what is happening, I recommend setting MOP(1)=1 (to see iteration information), MOP(5)=3 (to see phase changes),…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH2
- Reported - view
Re: How to consider the evaporation effect of water affected by temperature in tough2
Dear Sue, I am hesitant to send you the source code for the subroutine that calculates saturation pressure as a function of temperature in TOUGH3, because it is not designed to be used in…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: How to output viscosity in the output results
Use the OUTPU option in TOUGH3, described in the TOUGH3 User's Guide starting on page 83. See Table 8 for a description of all variables that can be output for each grid block.
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: How to consider the evaporation effect of water affected by temperature in tough2
@sue In the TOUGH3 source code, open file thrm.f90 and look for subroutine sat and the subroutines it calls.
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: How to consider the evaporation effect of water affected by temperature in tough2
All versions of TOUGH can simulate evaporation. The instructions you show give three options for how to specify a evaporation at the ground surface. More details on this topic may be found in the…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: Infiltration into dry soil
Dear Masha, When I look at your output file, I do not see any iteration failure. You asked the code to run 10 steps and it did. However, there appear to be problems in your input file.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH3
- Reported - view
Re: Initialization (T2Well)
I think that should work. You should remove the lateral connections between the monitoring well and the adjacent surface grid blocks that have been given a huge volume,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: How to consider the evaporation effect of water affected by temperature in tough2
Dear Sue, The partitioning of water into liquid and gaseous phases is determined by steam tables, which are built into TOUGH. For a given temperature and pressure condition in a grid block,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: T2Well (Error message)
@Refaat G Hashish When you bought T2Well did you buy just the executable or the executable plus the source code? If you just bought the executable, the instructions I gave you for changing the file…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Coupled reservoir-wellbore model initialization (using T2Well)
I am generally familiar with how T2Well works, but I have not actually used it. I recommend that you look at the user's guide, and examine the governing equations for wellbore flow to see whether…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Coupled reservoir-wellbore model initialization (using T2Well)
I think you are correct that T=1 is too small for the EOS you are using. SATB is just the name of the subroutine where the code fails, but the important thing here is that for some reason the…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Coupled reservoir-wellbore model initialization (using T2Well)
Try decreasing the wellbore permeability by a factor of 10. If that does not help decrease it more. You want wellbore permeability to be very high compared to reservoir permeability,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Coupled reservoir-wellbore model initialization (using T2Well)
I recommend making the permeability of the well material equal to the reservoir permeability for the initialization run. With no sources or sinks, it does not need to have high permeability.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: T2Well (Error message)
I don't think so. As I recall, providing output in excel format is a new feature with TOUGH3. For TOUGH2, there is a free utility program called Ext (available for download,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: T2Well (Error message)
T2Well has TOUGH2 built into it, so it can model coupled reservoir-wellbore problems. We do not currently have a "T3Well" where TOUGH3 is coupled to a wellbore simulator.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
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
- 5 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
- 5 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
- 5 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
- 5 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
- 5 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
- 5 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
- 5 yrs ago
- TOUGH2
- Reported - view

1 2 3 4 5 6

Back to top

Home

View all topics

TOUGH Website