Christine Doughty's Profile - TOUGH User Community

Latest Posts

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 agoMon, January 24, 2022 at 5:22 PM UTC
- TOUGH3
- Reported - view
Re: Trapping CO2 in ECO2N
Please see the thread "Residual and solubility trapping" in the TOUGH2 forum.
- Staff Scientist
- Christine_Doughty
- 3 yrs agoTue, November 23, 2021 at 3:55 PM UTC
- Questions & Answers
- Reported - view
Re: Compressibility On/Off in ECO2M
Pore compressibility describes how porosity changes with pressure change. You can turn it off by setting COM in the ROCKS block to zero. Then porosity will be constant (unless you have precipitation…
- Staff Scientist
- Christine_Doughty
- 3 yrs agoMon, November 22, 2021 at 6:50 PM UTC
- 1
- TOUGH3
- Reported - view
Re: Residual and solubility trapping
I have not used ECO2M in a while, so if anyone else has better information, please chime in. But here is what I suggest: For Solubility Trapping (mass of CO2 dissolved in the aqueous phase) there…
- Staff Scientist
- Christine_Doughty
- 3 yrs agoMon, November 22, 2021 at 6:06 PM UTC
- 2
- TOUGH2
- Reported - view
Re: Salt precipitation dependency on rate (and grid size)
Dear Godarzi, In my experience, precipitation of salt accompanying CO2 injection is one of the most grid-sensitive problems there is. With a coarse enough grid you may not see precipitation at all,…
- Staff Scientist
- Christine_Doughty
- 3 yrs agoWed, August 18, 2021 at 9:16 PM UTC
- TOUGH2
- Reported - view
Re: Output results of specific grid versus time
Dear Refaat, You should include the FOFT block in your TOUGH3 input. You specify a list of elements (grid cells) for which state variables at each time step are recorded in csv files.…
- Staff Scientist
- Christine_Doughty
- 3 yrs agoThu, July 15, 2021 at 5:03 PM UTC
- TOUGH3
- Reported - view
Re: Columns Headings in COFT (T2Well)
I don't know T2Well in particular, but in general for identifying COFT entries I make sure TOUGH is printing out at least one main output (search for the string OUTPUT DATA),…
- Staff Scientist
- Christine_Doughty
- 3 yrs agoWed, June 16, 2021 at 5:21 PM UTC
- TOUGH3
- Reported - view
Re: Specify the time-dependent BHP
I don't have too much experience using the DELV option, but I think the only way to have a decreasing bottomhole pressure over time would be to restart the problem at a series of times,…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoTue, May 18, 2021 at 12:34 AM UTC
- TOUGH2
- Reported - view
Re: TOUGH Workshop 2021
Yes, we are tentatively planning a remote Beginner's Short Course for this October or November, and possibly other short courses also. Contact Helen Prieto (HGPrieto@lbl.…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoMon, April 12, 2021 at 9:55 PM UTC
- Questions & Answers
- Reported - view
Re: The order of ELEME and CONNE will influence computation results
I suspect there are other TOUGH users who can give you a much more technical reply, but the general guidance I have been given about constructing grids is to be systematic about ordering elements…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoThu, April 1, 2021 at 3:34 AM UTC
- TOUGH+
- Reported - view
Re: Alternative for running immiscible brine-CO2 system
I think what you are asking for is a system where you have brine and supercritical CO2, and no brine evaporates into the CO2 phase and no CO2 dissolves into the aqueous phase.…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoThu, April 1, 2021 at 3:22 AM UTC
- TOUGH2
- Reported - view
Re: How to account for evaporation Sw initial < Swir
Hi Carla, If you represent infiltration as a mass source of water, then Sw will increase at your infiltration locations and eventually Sw > Swir and liquid will be mobile.…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoMon, March 15, 2021 at 4:40 PM UTC
- TOUGH2
- Reported - view
Re: How does one track internal progress while the solver is iterating?
I am able to edit the output file while the code is running in an Ubuntu window. I use the text editor vi. I get a warning message when the file is updated,…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoFri, February 19, 2021 at 12:26 AM UTC
- TOUGH3
- Reported - view
Re: Non-isothermal solution crashes
Yannick, According to the output file you attached, the code completes one step and then crashes while printing output. Near the top of the output file are some error messages regarding the input…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoWed, February 17, 2021 at 12:32 AM UTC
- TOUGH3
- Reported - view
Re: CANNOT FIND PARAMETERS AT ELEMENT *00001*
I just thought of one more thing. Your first element name is 00000. This makes me kind of nervous, since TOUGH can read zeros as blanks and vice versa.…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoWed, February 10, 2021 at 9:17 PM UTC
- TOUGH3
- Reported - view
Re: CANNOT FIND PARAMETERS AT ELEMENT *00001*
I think there is an inconsistency in your input file. I haven't used ECO2M for a while, and never in TOUGH3, but I found the following in the User's Guide.…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoWed, February 10, 2021 at 9:11 PM UTC
- TOUGH3
- Reported - view
Re: CANNOT FIND PARAMETERS AT ELEMENT *00001*
I notice in your output file that files MESHA and MESHB exist. This should only occur when you have not changed your ELEME and CONNE blocks (or MESH file) from the previous time you ran TOUGH3.…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoWed, February 10, 2021 at 7:00 PM UTC
- TOUGH3
- Reported - view
Re: injecting VOC mass through virtual cells
Dear Alireza, Your input file does not match what your message says. You do not have any mass sources- these would be specified in the GENER block and yours is empty.…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoMon, January 25, 2021 at 9:09 PM UTC
- TMVOC
- Reported - view
Re: T2Well-ECO2N
TOGA does handle three-phase flow, and when I asked Lehua about the possibility of using it for CO2 in liquid, gaseous, and supercritical phases (like ECO2M), he said that formally it should work,…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoFri, January 8, 2021 at 10:12 PM UTC
- TOUGH2
- Reported - view
Re: Not to achineve the end time and no error
Whenever a code is not running well, it is helpful to set these MOPs to learn more about what is happening. Look at your output file and look for the string " ST ".…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoFri, January 8, 2021 at 6:10 PM UTC
- 1
- TOUGH2
- Reported - view
Re: T2Well-ECO2N
Dear Godarzi, To help you decide whether you need ECO2M versus ECO2N, below is the phase diagram of CO2. If your problem will not cross the saturation line (shown in black) or get near the critical…
- Staff Scientist
- Christine_Doughty
- 4 yrs agoFri, January 8, 2021 at 3:45 AM UTC
- TOUGH2
- Reported - view
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 agoFri, January 8, 2021 at 3:04 AM UTC
- 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 agoMon, December 21, 2020 at 4:09 PM UTC
- 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 agoMon, December 21, 2020 at 3:32 PM UTC
- 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 agoTue, December 15, 2020 at 5:06 PM UTC
- TOUGH2
- Reported - view

1 2 3 4 5 6

Christine Doughty

Re: Can TOUGH3 model temperature below zero degree?

Re: Trapping CO2 in ECO2N

Re: Compressibility On/Off in ECO2M

Re: Residual and solubility trapping

Re: Salt precipitation dependency on rate (and grid size)

Re: Output results of specific grid versus time

Re: Columns Headings in COFT (T2Well)

Re: Specify the time-dependent BHP

Re: TOUGH Workshop 2021

Re: The order of ELEME and CONNE will influence computation results

Re: Alternative for running immiscible brine-CO2 system

Re: How to account for evaporation Sw initial < Swir

Re: How does one track internal progress while the solver is iterating?

Re: Non-isothermal solution crashes

Re: CANNOT FIND PARAMETERS AT ELEMENT 00001

Re: CANNOT FIND PARAMETERS AT ELEMENT 00001

Re: CANNOT FIND PARAMETERS AT ELEMENT 00001

Re: injecting VOC mass through virtual cells

Re: T2Well-ECO2N

Re: Not to achineve the end time and no error

Re: T2Well-ECO2N

Re: Not to achineve the end time and no error

Re: How to consider the evaporation effect of water affected by temperature in tough2

Re: How to output viscosity in the output results

Re: How to consider the evaporation effect of water affected by temperature in tough2

Home