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