[go: up one dir, main page]
Menu ▾ ▴

[8231f4]: / dev / Tickets / 114.py  Maximize  Restore  History

Download this file

77 lines (61 with data), 2.7 kB 

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
import matplotlib

matplotlib.use('Qt4Agg')



import numpy 

import CoolProp

import CoolProp.CoolProp as CP

import matplotlib.pyplot as plt

from CoolProp.CoolProp import FluidsList



print "Testing TTSE for revision "+CoolProp.__gitrevision__



def compareProperty(fluid="",p=0,what=""):

    global c_diff, c_unit, c_exce

    

    if p==0:

        p = 0.75*CP.Props(fluid,"pcrit")

  

    Delta_T = 50

    T_bub = CP.Props("T","P",p,"Q",0,fluid)

    T_dew = CP.Props("T","P",p,"Q",1,fluid)

    h_bub = CP.Props("H","P",p,"Q",0,fluid)

    h_dew = CP.Props("H","P",p,"Q",1,fluid)

    

    T_1 = T_bub-Delta_T

    if T_1 < CP.Props(fluid,"Tmin"):

        T_1 = CP.Props(fluid,"Tmin")+0.5*(T_bub-CP.Props(fluid,"Tmin"))

    T_2 = T_dew+Delta_T

    

    h_1 = CP.Props("H","P",p,"T",T_1,fluid)

    h_2 = CP.Props("H","P",p,"T",T_2,fluid)

  

    h_liq = numpy.linspace(h_1,h_bub,num=100)

    h_vap = numpy.linspace(h_dew,h_2,num=100)

    

    T_liq = CP.Props("T","P",p,"H",h_liq,fluid)-T_bub

    T_vap = CP.Props("T","P",p,"H",h_vap,fluid)-T_dew

    

    X_liq_STDV = CP.Props(what,"P",p,"H",h_liq,fluid)

    X_vap_STDV = CP.Props(what,"P",p,"H",h_vap,fluid)

    

    CP.enable_TTSE_LUT(fluid)

    X_liq_TTSE = CP.Props(what,"P",p,"H",h_liq,fluid)

    X_vap_TTSE = CP.Props(what,"P",p,"H",h_vap,fluid)

    

    if numpy.max([X_liq_STDV/X_liq_TTSE,X_vap_STDV/X_vap_TTSE])>1.25 or numpy.min([X_liq_STDV/X_liq_TTSE,X_vap_STDV/X_vap_TTSE])<0.75:

        c_diff += 1

        print ""

        print "There were problems with "+what+" for "+fluid

        print "Relative difference liquid: "+str(numpy.mean((X_liq_STDV-X_liq_TTSE)/X_liq_STDV))

        print "Relative difference vapour: "+str(numpy.mean((X_vap_STDV-X_vap_TTSE)/X_vap_STDV))

        print "Average factor liquid: "+str(numpy.mean(X_liq_STDV/X_liq_TTSE))

        print "Average factor vapour: "+str(numpy.mean(X_vap_STDV/X_vap_TTSE))

        #plt.plot(numpy.append(T_liq,T_vap),numpy.append(X_liq_STDV,X_vap_STDV),label=what+", standard")

        #plt.plot(numpy.append(T_liq,T_vap),numpy.append(X_liq_TTSE,X_vap_TTSE),label=what+", TTSE")

        ##plt.show(block=True)

        #plt.savefig("/home/jowr/tmp/viscosity/"+fluid+".png")

        #plt.clf()

    



#toTest  = ["L","V"]

toTest  = "V"

#fluids = ["n-Pentane","R134a"]

fluids = FluidsList()

c_diff = 0

c_unit = 0

c_exce = 0



for fluid in fluids:

    try: 

        compareProperty(fluid=fluid,what=toTest)

    except ValueError:

        c_exce += 1

        print "An exception occurred for "+toTest+" with "+fluid



print "Finished testing TTSE:"

print "Errors occurred in "+str(c_exce)+" out of "+str(len(fluids))+" fluids"

print "and differences occurred in "+str(c_diff)+" fluids."