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.

8replies Oldest first
  • Oldest first
  • Newest first
  • Active threads
  • Popular
  • 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

    Like
  • Kindly find the attached input file

    Like
  • Yashesh, it also appears you're trying to run this in parallel.  Are you using TOUGH3, by chance?

    Like
    • Mikey Hannon Yes I am using TOUGH 3

      Like
  • 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

    Like
  • 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

    Like
    • Yingqi Zhang Thanks a lot it works now. I really appreciate your help

      Like
  • Thank you, everyone, I really appreciate the help. 

    Like
Like Follow
  • 9 mths agoLast active
  • 8Replies
  • 56Views
  • 3 Following