from CoolProp import CoolProp as CP
from PDSim.misc.datatypes import Collector
import numpy as np
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
from scipy.odr import *
import textwrap
fluid_REF = 'Propane'
Tcrit_REF = CP.Props(fluid_REF, 'Tcrit')
omega_REF = CP.Props(fluid_REF, "accentric")
molemass_REF = CP.Props(fluid_REF, 'molemass')
rhocrit_REF = CP.Props(fluid_REF, 'rhocrit')
Zcrit_REF = CP.DerivTerms('Z', Tcrit_REF, rhocrit_REF, fluid_REF)
fluid = 'DimethylEther'
molemass = CP.Props(fluid, 'molemass')
Ttriple = CP.Props(fluid, 'Ttriple')
Tcrit = CP.Props(fluid, 'Tcrit')
omega = CP.Props(fluid, "accentric")
rhocrit = CP.Props(fluid, 'rhocrit')
pcrit = CP.Props(fluid, 'pcrit')
Zcrit = CP.DerivTerms('Z', Tcrit, rhocrit, fluid)
N = 12
RHO, TTT, RHO0, TTT0 = Collector(), Collector(), Collector(), Collector()
rhomax = CP.Props('D', 'T', Ttriple, 'Q', 0, fluid)
#Build a database of "experimental" data
for T in np.linspace(Ttriple, Tcrit + 50, 80):
"phir", "dphir_dTau", "d2phir_dTau2", "dphir_dDelta", "d2phir_dDelta2",
"d2phir_dDelta_dTau", "d3phir_dDelta2_dTau", "phi0", "dphi0_dTau",
"d2phi0_dTau2", "dphi0_dDelta", "d2phi0_dDelta2",
"IsothermalCompressibility"
]
multiply = ["VB", "dBdT", "VC", "dCdT", "dphir_dDelta", "dpdT"]
keys = multiply[:]
fluid = "n-Pentane"
T = numpy.array([30, 100, 150, 210]) + 273.15
D = numpy.array([710, 20, 20, 210])
for key in keys:
liquid = CP.DerivTerms(key, T[0], D[0], fluid)
twophase = CP.DerivTerms(key, T[1], D[1], fluid)
gaseous = CP.DerivTerms(key, T[2], D[2], fluid)
supercrit = CP.DerivTerms(key, T[3], D[3], fluid)
print '{:<25}: {:>10.5f}; {:>10.5f}; {:>10.5f}; {:>10.5f}'.format(
key, liquid, twophase, gaseous, supercrit)
for key in keys:
liquid = CP.DerivTermsU(key, T[0], D[0], fluid, 'SI')
twophase = CP.DerivTermsU(key, T[1], D[1], fluid, 'SI')
gaseous = CP.DerivTermsU(key, T[2], D[2], fluid, 'SI')
supercrit = CP.DerivTermsU(key, T[3], D[3], fluid, 'SI')
print '{:<25}: {:>10.5f}; {:>10.5f}; {:>10.5f}; {:>10.5f}'.format(
key, liquid, twophase, gaseous, supercrit)
print
