Christine Doughty's Profile - TOUGH User Community

Skip to main content

Latest Posts

Re: Incomplete simulation run
Please check that the radial spacing of your grid is what you really want. You start with a well radius of 0.3 m, and according to the MESHMAKER input appended to your input file,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
TOUGH3 on PC: CPU timing routine

I recently used Ubuntu terminal to install TOUGH3 on a Windows 10 Laptop. In order to compile, I had to change a call to the CPU timing routine. In the esd-toughlib source code, file parallel_stubs.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1reply
- 142
- TOUGH3
Re: TOUGH+RealGasBrine for CO2 injection
@Yashvardhan Verma I don't know if this will work for your problem, but when I wanted to inject supercritical CO2 into water and have no dissolution, I specified the water to be saturated with CO2…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH+
- Reported - view
Re: How to output all primary and secondary variables using the OUTPU data block in TOUGH3?
Dear Refaat, It looks like you have some tabs in your input file. In general, tabs give TOUGH3 trouble and should be removed. The variable names should begin in column 1.…
Attachment
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH3
- Reported - view
Re: Steady-state simulation of a low-type geothermal system
Dear Hakim, I might be confused, but it looks like you are using the two-waters option for EOS1, since your MULTI block is 2 3 2 6. According to the User's Guide,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH2
- Reported - view
Re: Steady-state simulation of a low-type geothermal system
Dear Hakim, I ran a sample problem injecting into a closed volume with EOS1 and the pressure was able to go way above 800 bars without any obvious problems,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH2
- Reported - view
Re: Steady-state simulation of a low-type geothermal system
Dear Hakim, The file you attached T1.dat is the correct format for a TOUGH input file. However, since I haven't used Petrasim for a really long time, I cannot be sure this is the input file that…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH2
- Reported - view
Re: Steady-state simulation of a low-type geothermal system
Dear Hakim, I looked at your Excel file, and as far as I can tell, it all looks reasonable. However, to have a chance of helping you, I need to look at the actual TOUGH input file.…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH2
- Reported - view
Re: Steady-state simulation of a low-type geothermal system
Dear Saibi, I downloaded your zip file and extracted what should be the input file, but when I tried to read it I found it is not an ASCII (i.e. plain text) file,…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- TOUGH2
- Reported - view
Re: Tough2-FLAC 3D
I forwarded your query to Jonny Rutqvist (the person I know who knows the most about TOUGH/FLAC), and here is his reply: I do not have a routine that convert TOUGH2 coordinates (mid-element…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- 1
- Questions & Answers
- Reported - view
Re: Issue with using EXT for extracting data from TOUGH2 MESH and OUTPUT files
Hmm, I am not sure Ext will work for such a big mesh. If anyone else has experience with this, please chime in. The compilation error messages look to me like your computer memory is not big enough…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: Issue with using EXT for extracting data from TOUGH2 MESH and OUTPUT files
I can't tell from the error message what is wrong. Could you post your MESH and OUTPUT files? If they are too big, try creating a sample problem with a smaller MESH and output at just one time in…
- Staff Scientist
- Christine_Doughty
- 4 yrs ago
- TOUGH2
- Reported - view
Re: How to print K(LIQ) and K(GAS)
You could also try setting KDATA=3 (column 4 in the first line of the PARAM block). I can't remember about TOUGH2-ECO2N for sure, but for some TOUGH versions this adds an extra section at each of…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- TOUGH2
- Reported - view
Re: Unexpected gas phase evolving
During a simulation with no gas phase initially present, if pressure drops somewhere a gas phase may form. The process could be water boiling, or a dissolved species coming out of solution.…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- TOUGH+
- Reported - view
Re: No molecular diffusion and hydrodynamic dispersion in T2VOC
I am not really an expert in T2VOC, so others may want to answer also, but if you want molecular diffusion in all phases, you should use TMVOC instead of T2VOC.…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- T2VOC
- Reported - view
Re: Executable or Source code?
Dear Diego, It's not a silly question, and I hope I can help you understand the difference. If you buy the executable only, you will be able to run the code with all the options available in the…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- TOUGH2
- Reported - view
Re: Set An Initial Temperature Greater Than 100 In ECO2N Module
Dear Henry, There are other limitations in ECO2N for T>103.04 C besides just CO2TAB. If you want T=120 C, you need to use ECO2N V2. Check to see if your current version of TOUGH has ECO2N V2.…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- TOUGH2
- Reported - view
Re: Failure in linear equation soluation
Dear Alborz, In my experience, a variety of problems could result in this error message. If you set MOP(1)=1, you will get a one-line message every time step that gives you information about how the…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- TOUGH3
- Reported - view
Re: Hysteretic RELP and CAP
Hmm, I have had trouble with convergence in EOS7c from time to time also. If you have not already tried it, I suggest trying IRP=12 and ICP=12 and turning off the hysteresis by setting…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- 1
- TOUGH2
- Reported - view
Re: Hysteretic RELP and CAP
Hi Ramin, What version of TOUGH are you using? My TOUGH3 does not have ICP=9 as an option. In looking at the source code, I see that if you pick an ICP that does not exist,…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- 1
- TOUGH2
- Reported - view
Re: Partitioning Tracer-EOS7C
@Ramin Dear Ramin, Yes, that is what I meant: the Henry's coefficient for the tracer component does not change with temperature. Whether it should change with temperature depends on what tracer you…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- 1
- TOUGH2
- Reported - view
Re: Partitioning Tracer-EOS7C
Dear Ramin, 1. See the EOS7C User's Manual, Section 2.2, in particular, the text leading up to Equation (13), which describes how the Henry's coefficient is calculated for the NCG component…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- 1
- TOUGH2
- Reported - view
Re: Five-spot Geothermal production/injection
@Alfredo Battistelli Hi Alfredo, This is the first time I've heard of PATTY, and it sounds very useful. Do you have the source code and/or a copy of LBL-20739 (which I searched for here but cannot…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- TOUGH2
- Reported - view
Re: Non-convergence when upper boundary saturated in CO2
Hi Nathalie and Mikey, I tried running the problem with TOUGH2/ECO2N V2 on my Windows 7 32-bit laptop. It failed after 12-15 steps when RE1 was 1.e-5 or 1.e-6. With RE1=1.…
- Staff Scientist
- Christine_Doughty
- 5 yrs ago
- TOUGH2
- Reported - view
Re: anisotropic permeability for TOUGH3-EOS7
Dear Cheng-Kuo, For now, I want to wait until Yoojin can respond to think about Cases 2 and 3. For case 1, you said it ran okay. From what I can tell, you have mop2(20)=0 so the code should take the…
- Staff Scientist
- Christine_Doughty
- 6 yrs ago
- TOUGH2
- Reported - view

1 2 3 4 5 6

Back to top

Home

View all topics

TOUGH Website