Files
CoolProp/dev/Tickets/SF-89.py
2014-05-14 12:46:24 +02:00

34 lines
1.1 KiB
Python

'''
Created on 20 Sep 2013
@author: jowr
'''
# New example with R407F mixture
from pyrp.refpropClasses import RefpropSI
import CoolProp.CoolProp as cp
p = 30000
T = 273.15
ref = False
if ref:
xkg=[0.473194694453358,0.205109095413331,0.321696210133311]
names="R32|R125|R134a"
RP = RefpropSI()
RP.SETUPFLEX(xkg=xkg, FluidNames=names)
T_A,p_A,D_A,Dl_A,Dv_A,q_A,e_A,h_A,s_A,cv_A,cp_A,w_A = RP.PQFLSH(p, 0)
T_B,p_B,D_B,Dl_B,Dv_B,q_B,e_B,h_B,s_B,cv_B,cp_B,w_B = RP.PQFLSH(p, 1)
T_C,p_C,D_C,Dl_C,Dv_C,q_C,e_C,h_C,s_C,cv_C,cp_C,w_C = RP.TQFLSH(T, 0)
hlb = h_A/1000.
hrb = h_B/1000.
h200 = h_C/1000.
print "Refprop: ", hlb, hrb, h200
else:
R407F='REFPROP-MIX:R32[0.473194694453358]&R125[0.205109095413331]&R134a[0.321696210133311]'
#R407F='REFPROP-MIX:R32[0.651669604033581]&R125[0.122438378639971]&R134a[0.225892017326446]'
hlb=cp.Props('H','P',30,'Q',0,R407F) # 30 kPa saturated liquid
hrb=cp.Props('H','P',30,'Q',1,R407F) # 30 kPa saturated vapour
h200=cp.Props('H','T',273.15,'Q',0,R407F) # saturated liquid at 0C IIR
print "CoolProp: ", hlb, hrb, h200