verificationTable.py

   1 import XSteam #unit XSteam
   2 import refprop #unit refprop
   3 import fluid #unit SI
   4
   5
   6
   7 #SI
   8 x = 4e6
   9 y = 793.15
  10
  11 #------------------------------------------------------------------------------
  12
  13 fluid = fluid.FluidREFPROP(None,'WATER')
  14 fF = fluid.h_pT(x,y)
  15 '''fF = fluid.h_ps(x,y)
  16 fF = fluid.h_px(x,y)
  17 fF = fluid.h_Tx(x,y)
  18 fF = fluid.s_pT(x,y)
  19 fF = fluid.s_ph(x,y)
  20 fF = fluid.x_ph(x,y)
  21 fF = fluid.x_ps(x,y)'''
  22 resultF = fF
  23 print fF
  24
  25 #------------------------------------------------------------------------------
  26 wm  = refprop.function('wm',0,0,'WATER')
  27 b = wm/1000.0
  28 '''
  29 f = refprop.function('tsat_p',100,0,'WATER') #GOOD
  30 f = refprop.function('tsat_s',5.5*1e3*b,0,'WATER')#GOOD
  31 f = refprop.function('t_ph',4000,45.912*1e3*b,'WATER') #GOOD T:[10:500]+273.15 K p:[1:40]bar
  32 f = refprop.function('t_ps',100,8.45*1e3*b,'WATER')#GOOD
  33 f = refprop.function('t_hs',457*1e3*b,1.5*1e3*b,'WATER')  #difference legere with xsteam, mais probablement dans des zones speciales
  34 f = refprop.function('psat_t',373.15,0,'WATER') #GOOD
  35 f = refprop.function('psat_s',4.8*1e3*b,0,'WATER')#Good meme si legere difference.
  36 f = refprop.function('p_hs',3800*1e3*b,7.5*1e3*b,'WATER')#Ne fonctionne que dans certaine zone !
  37 '''
  38
  39 f = refprop.function('h_pt',x/1e3,y,'WATER')/b
  40 '''f = refprop.function('h_ps',x/1e3,y*b,'WATER')/b
  41 f = refprop.function('h_px',x/1e3,y,'WATER')/b
  42 f = refprop.function('h_Tx',x,y,'WATER')/b
  43 f = refprop.function('s_pT',x/1e3,y,'WATER')/b
  44 f = refprop.function('s_ph',x/1e3,y*b,'WATER')/b
  45 f = refprop.function('x_ph',x/1e3,y*b,'WATER')
  46 f = refprop.function('x_ps',x/1e3,y*b,'WATER')'''
  47 result = f
  48
  49
  50 #------------------------------------------------------------------------------
  51 fX = XSteam.function('h_pT',x/1e5,y-273.15)
  52 #fX = XSteam.function('h_ps',x/1e5,y/1e3)
  53 #fX = XSteam.function('h_px',x/1e5,y)
  54 #fX = XSteam.function('h_Tx',x-273.15,y)
  55 #fX = XSteam.function('s_pT',x/1e5,y-273.15)
  56 #fX = XSteam.function('s_ph',x/1e5,y/1e3)
  57 #fX = XSteam.function('x_ph',x/1e5,y/1e3)
  58 #fX = XSteam.function('x_ps',x/1e5,y/1e3)
  59 resultX = fX
  60
  61
  62
  63
  64 #print result
  65 print 'refprop :',result, 'xsteam :',resultX,'fluid.py :',resultF
  66
  67
  68 '''
  69 h = j/mole
  70
  71 wm = g/mole
  72                 b = kg/mole = wm/a
  73 a = 1000 g/kg
  74
  75 j/kg => h/b = h*a/wm
  76 '''