Christine Doughty's Profile - TOUGH User Community

Skip to main content

Latest Posts

Re: EOS7C WARNING. H2O MASS FRACTION IN THE GAS WAS LESS THAN ZERO IN 1 GRID BLOCKS
In the past, I have had better success using EOS7C when I specify a very small amount of all components as initial conditions. I use something like 1e-10 for mass fractions,…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH3
- Reported - view
Re: Relative permeability curves in tough2
If none of the existing relative permeability options in TOUGH2 provide Krl and Krg that you like, then you will have to create your own. You could take the Krl function from one IRP option and the…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH2
- Reported - view
Re: Gridblocks to represent sink and source
There are a few problems with your GENER file: 1. A blank line means you are all finished reading GENER, so move the one you have now near the top of the file to the end. 2.…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH+
- Reported - view
Re: Gridblocks to represent sink and source
Dear Mahmud, At first I thought that the entries in the GENER block in your input file were just in the wrong columns. Here is what I thought it should look like: GENER A2R00Inj01 144 76 1 COM1 2.…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- 1
- TOUGH+
- Reported - view
Re: Gridblocks to represent sink and source
I think you want one well to be constant flow rate (an injection well) and the other well to be constant pressure (a production well), is that right? One possibility is to specify one mass source in…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH+
- Reported - view
Re: SAVE file not consistent with OUTPUT
Dear Haiyan, I am surprised to see slight differences in P and T between the SAVE and OUTPUT files. Can you verify that the SAVE file is actually for the same time as the output file (kcyc=17, t=0.…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH3
- Reported - view
Re: meshmaker error, increments problem
Oh, that's right - you mentioned TOUGH+ in your first post. Since we were using the TOUGH3 forum category I forgot that. I do not have much TOUGH+ experience.…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH3
- Reported - view
Re: meshmaker error, increments problem
You are on the right track with your GENER block, but there are a couple of things to check. First, if you want to inject heat (units Watts, like with an underground heater),…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH3
- Reported - view
Re: meshmaker error, increments problem
There are a bunch of lines at the start of your input file that TOUGH does not understand. That is what the message "input conversion error" is telling you.…
Attachment
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH3
- Reported - view
Re: Confining beds option geothermal doublet
Hi Virginie, If you post your input file and a citation for the Gringarten solution you are trying to match, I will take a look and see if I can suggest any changes that might give you a better…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH2
- Reported - view
Re: No convergence in NSIN at - non-isothermal simulation T2WELL/ECO2N
The code actually recovers from the error message you show, and continues running. The outputs at 5 and 10 seconds look pretty normal. But at about 18 seconds,…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH2
- Reported - view
Re: applying heat load in BHE and varying input T?
Hi Hannah, Your heat source term looks correct (when I display it with a uniform-width font, so I can check the column alignment). Check the temperature in element A1000,…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH2
- Reported - view
Re: Format of INCON and SAVE file
@Jianhua Jiang Unfortunately, the bug is such that TOUGH does think it is using the correct permeability, so if you use the OUTPU block to request that permeability be shown for each block,…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH3
- Reported - view
Re: Format of INCON and SAVE file
@Jianhua Jiang If you have the source code, check that the bug in MULTI.f90 has been fixed (see my previous post and the TOUGH Bugs and Fixes (https://tough.lbl.…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- TOUGH3
- Reported - view
Re: Relative permeability curves in tough2
If you are using TOUGH3 (or iTOUGH I think), the option IRP=12, ICP=12 invokes hysteretic relative permeability and capillary pressure, which includes the option for Krg_max<1 (with RP(5)).…
- Staff Scientist
- Christine_Doughty
- 2 yrs ago
- 1
- TOUGH2
- Reported - view
Re: Using TOUGH ECO2N for Core flooding
I agree with all Kenny's suggestions, and in addition you could try the following: 1. Different relative permeability function - liquid relative permeability krl controls how much water can flow…
- Staff Scientist
- Christine_Doughty
- 3 yrs ago
- 1
- TOUGH2
- Reported - view
Re: EOS7c and NCG injection
I would need to see more of your output to know for sure, but I have experienced this sort of thing also. In the original EOS7c, it was a strict requirement that some CH4 always be present,…
- Staff Scientist
- Christine_Doughty
- 3 yrs ago
- TOUGH3
- Reported - view
Re: How to install TOUGH3 on Window OS?
I have compiled and run TOUGH3 in an Ubuntu window under Windows 10. There should be a script file in your TOUGH3 release package called compile_T3_UbuntuTerminal.sh This script calls configure.…
- Staff Scientist
- Christine_Doughty
- 3 yrs ago
- Questions & Answers
- Reported - view
Re: Changes of hydraulic conductivity driven by effective stress
TOUGH has been coupled to the commercial software FLAC and also another code ROCMECH, to do coupled fluid flow and transport and geomechanics. If you look at the list of TOUGH publications at https:…
- Staff Scientist
- Christine_Doughty
- 3 yrs ago
- Questions & Answers
- Reported - view
Re: MeshMaker: Element depth
Dear Ludovic, Sorry for the late reply, I guess I have not been checking the Forum regularly enough. I think you are correct about MeshMaker doing the right thing,…
- Staff Scientist
- Christine_Doughty
- 3 yrs ago
- Pre & Post Processors
- Reported - view
Re: Calculating Total Heat Generated in an Area
@Christine Doughty Oops, I said the wrong thing about mop(5)=9 output. First is the gas phase, then the liquid phase.
- Staff Scientist
- Christine_Doughty
- 3 yrs ago
- TOUGH2
- Reported - view
Re: Calculating Total Heat Generated in an Area
@Christine Doughty One more thing I forgot. The fluid specific enthalpies use a reference temperature of 0 C, so for consistency that should also be used with rock internal energy.…
- Staff Scientist
- Christine_Doughty
- 3 yrs ago
- TOUGH2
- Reported - view
Re: Calculating Total Heat Generated in an Area
Dear Alborz, Sorry for the late reply - I forgot to check the Forum for a few days. My first thought when answering your original query was to try to use heat flux,…
- Staff Scientist
- Christine_Doughty
- 3 yrs ago
- TOUGH2
- Reported - view
Re: Calculating Total Heat Generated in an Area
This is kind of an ad hoc approach, but if pressure changes are small, so little heat is being created by mechanical work, then you might be able to approximate the heat content of each grid block…
- Staff Scientist
- Christine_Doughty
- 3 yrs ago
- TOUGH2
- Reported - view
Re: Can TOUGH3 model temperature below zero degree?
TOUGH3 cannot handle temperatures below zero degrees. In fact, using default parameters, it cannot handle temperatures below 1 C. If you set mop2(11)=1 (see Table 7 in User's Manual),…
- Staff Scientist
- Christine_Doughty
- 3 yrs ago
- TOUGH3
- Reported - view

1 2 3 4 5 6

Back to top

Home

View all topics

TOUGH Website