mirror of
https://github.com/CoolProp/CoolProp.git
synced 2026-01-15 00:48:18 -05:00
35 lines
1.2 KiB
Python
35 lines
1.2 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: %s %s %s" % (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: %s %s %s" % (hlb, hrb, h200))
|