0

TOUGH + Phase equilibrium Relationship...where am i wrong?

TOUGH+ has two phase equilibrium relationship

1. Kamath

2. Moridis

I tried plotting  Moridis Generalized equation taken for the manual and i find that its plot is different from the pic in the manual

Why is it so, am i missing something here?

I am trying to determine hydrate stability zone so its important to me

3 replies

null
    • Castleberry
    • 1 yr ago
    • Reported - view

    Here is the update, due to copying and pasting there were errors in the code

    Finally i got the correct relationship as shown in the picture, Do you guys think this is correct representation of these two relationships 

      • Castleberry
      • 1 yr ago
      • Reported - view

       Update:

      At T= 13.67celcius

      the moridis equation yields Peq as 10.8988MPa

      but when i use T+H to run simulation, it gives error for element 541 with T=13.67 Celsius stating hydrate cannot exist,, equilibrium pressure is 13.28369MPa

      I know that initialization is erroneous since hydrate cannot remain stable a this condition, but he problem is when i use the equation in the manual it yield 10.98MPa at 13.677 but when i run it says higher pressure is needed to make it stable. Could there be error in the manual equation?

       

      ERROR:

      Gas phase evolves at element "  541   " ==> Px  =  1.21189961E+07  Peq =  1.32883698E+07
       !!!!!!!!!!   Routine Phases2_AqH_Eql: ERRONEOUS DATA INITIALIZATION   ==>  STOP EXECUTION <<<<<<<<<<

    • kenny
    • 1 yr ago
    • Reported - view

    Have you compared your equilibrium pressure plot with the the figure presented in the T+H user manual? I guess something wrong with your calculation for the equilibrium pressure.

Content aside

  • 1 yr agoLast active
  • 3Replies
  • 50Views
  • 2 Following