TEMPERATURE = 0.000000E+00 OUT OF RANGE IN SAT
Hi All,
I am a newbie to the TOUGH simulator and trying to create a small model for brine injection using EOS7 module.
While I am trying to run the model, I get the following error message
TOUGH STATUS: Preprocessing step completed.
TOUGH STATUS: Total preprocessing time = 0.000000000000000E+000
parallel_abort(): EOS7: DATA INITIALIZATION ERROR
abort: Fortran Abort Called
Also, my output file shows the following:
MOLECULAR DIFFUSIVITY OF WATER, BRINE, AND AIR THROUGH THE GASEOUS AND AQUEOUS PHASES, (FDDIAG(PHASE,COMP)) [M**2/S]:
PHASE 1 = GAS; PHASE 2 = AQUEOUS
PHASE COMP PHASE COMP PHASE COMP PHASE COMP PHASE COMP PHASE COMP
-1- -1- -1- -2- -1- -3- -2- -1- -2- -2- -2- -3-
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
TEMPERATURE = 0.000000E+00 OUT OF RANGE IN SAT
+++++++++ CANNOT FIND PARAMETERS AT ELEMENT *A11 1* XX(M) = 0.000000E+00 0.000000E+00 0.000000E+00
!!!!!!!!!!! ERRONEOUS DATA INITIALIZATION !!!!!!!!!! STOP EXECUTION----------
From the output file, I understand that I need to change the value for Temperature but do not know where I can make that change (I believe it must be in PARAM or GENER)
I understand this is a small problem for the experienced people here but your suggestion will help me in learning the simulator and I am thankful for the suggestion in advance.
8 replies
-
Hello Yashesh,
From the error message it seems your initialization of the problem is not physical. Can you attach your input file? Please include the ELEM, CONNE and INCON blocks.
Thank you
Yingqi
-
Kindly find the attached input file
-
Yashesh, it also appears you're trying to run this in parallel. Are you using TOUGH3, by chance?
-
The primary variables for eos7 are (in sequence):
P - PRESSURE T - TEMPERATURE X - AIR MASS FRACTION S+10. - (GAS PHASE SATURATION + 10.)
In your param block you specified 0 for temperature for all elements that are not in INCON block. Please also check if your INCON block specifications are correct.
Yingqi
-
Sorry Yashesh, I should re-explain your input. What I said was not exactly correct.
Based on what you had in Multi: (2,2,2,6) , you have three primary VARIABLES which are (P, XB, T)
In the param block you set T=0. INCON block you did not specify T, so T=0 for those elements as well. I replaced the third variable by 25.0, removed INCON. It runs fine (you forgot the decimal point in the total simulation time).
How I figured out was I run your input file and saw the output:
OPTIONS SELECTED ARE: (NK,NEQ,NPH,NB) = (2,2,2, 6)
NK = 2 - NUMBER OF FLUID COMPONENTS
NEQ = 2 - NUMBER OF EQUATIONS PER GRID BLOCK
NPH = 2 - NUMBER OF PHASES THAT CAN BE PRESENT
NB = 6 - NUMBER OF SECONDARY PARAMETERS (OTHER THAN COMPONENT MASS FRACTIONS)For NB = 6, diffusion is "off", for NB = 8, diffusion is "on"
AVAILABLE OPTIONS for (NK,NEQ,NPH,NB):
(3,3,2,6 or 8) - WATER, BRINE, AIR; ISOTHERMAL; VARIABLES (P, XB, X OR S+10, T)
(3,4,2,6 or 8) - WATER, BRINE, AIR; NON-ISOTHERMAL; VARIABLES (P, XB, X OR S+10, T)
(2,2,2,6 or 8) - WATER, BRINE, NO AIR; ISOTHERMAL; VARIABLES (P, XB, T)
(2,3,2,6 or 8) - WATER, BRINE, NO AIR; NON-ISOTHERMAL; VARIABLES (P, XB, T)Sorry for the confusing previously
Yingqi
-
Thank you, everyone, I really appreciate the help.