self.Tt = Tt
self.omega = omega
self.HEOS = HEOS
# Store the propoerties in a dict of CP_fluid instances
coolprop_fluids = {}
if has_CoolProp:
for CASRN in coolprop_dict:
HEOS = AbstractState("HEOS", CASRN)
coolprop_fluids[CASRN] = CP_fluid(
Tmin=HEOS.Tmin(),
Tmax=HEOS.Tmax(),
Pmax=HEOS.pmax(),
has_melting_line=HEOS.has_melting_line(),
Tc=HEOS.T_critical(),
Pc=HEOS.p_critical(),
Tt=HEOS.Ttriple(),
omega=HEOS.acentric_factor(),
HEOS=HEOS)
def CoolProp_T_dependent_property(T, CASRN, prop, phase):
r'''Calculates a property of a chemical in either the liquid or gas phase
as a function of temperature only. This means that the property is
either at 1 atm or along the saturation curve.
Parameters
----------
T : float
Temperature of the fluid [K]
self.Tt = Tt
self.omega = omega
self.HEOS = HEOS
# Store the propoerties in a dict of CP_fluid instances
coolprop_fluids = {}
if has_CoolProp:
for CASRN in coolprop_dict:
HEOS = AbstractState("HEOS", CASRN)
coolprop_fluids[CASRN] = CP_fluid(
Tmin=HEOS.Tmin(),
Tmax=HEOS.Tmax(),
Pmax=HEOS.pmax(),
has_melting_line=HEOS.has_melting_line(),
Tc=HEOS.T_critical(),
Pc=HEOS.p_critical(),
Tt=HEOS.Ttriple(),
omega=HEOS.acentric_factor(),
HEOS=HEOS)
class MultiCheb1D(object):
'''Simple class to store set of coefficients for multiple chebshev
approximations and perform calculations from them.
'''
def __init__(self, points, coeffs):
self.points = points
self.coeffs = coeffs
self.N = len(points) - 1
class Fluid:
"""Represents a fluid."""
def __init__(self, name: str, eos: str = 'HEOS'):
"""Initiate a Fluid instance.
:param ID: id of fluid
"""
self.name = name
self.eos = eos.upper()
self.state = AbstractState(self.eos, self.name)
if self.state.name() not in ["NitrousOxide"]:
raise NotImplementedError(
f'{self.state.name()} is not supported currently')
self.Tmax = 309.5 # self.state.Tmax()
self.Tmin = 182.23 # self.state.Tmin()
self.pmax = self.state.pmax()
self.pmin = 0
def copy(self):
fluid = Fluid(self.name, self.eos)
fluid.set_state(P=self.state.p(), T=self.state.T())
return fluid
def set_state(self,
D: float = None,
P: float = None,
T: float = None,
Q: float = None,
H: float = None,
S: float = None,
U: float = None):
args_present = {
'D': D is not None,
'P': P is not None,
'T': T is not None,
'Q': Q is not None,
'H': H is not None,
'S': S is not None,
'U': U is not None
}
if sum(args_present.values()) != 2:
raise ValueError(
f'Must have exactly 2 arguments: {sum(args_present.values())} provided'
)
args = []
if args_present['D']:
if args_present['P']:
args = CP.DmassP_INPUTS, D, P
elif args_present['T']:
args = CP.DmassT_INPUTS, D, T
elif args_present['Q']:
args = CP.DmassQ_INPUTS, D, Q
elif args_present['H']:
args = CP.DmassHmass_INPUTS, D, H
elif args_present['S']:
args = CP.DmassSmass_INPUTS, D, S
elif args_present['U']:
args = CP.DmassUmass_INPUTS, D, U
elif args_present['P']:
if args_present['T']:
args = CP.PT_INPUTS, P, T
elif args_present['Q']:
args = CP.PQ_INPUTS, P, Q
elif args_present['H']:
args = CP.HmassP_INPUTS, H, P
elif args_present['S']:
args = CP.PSmass_INPUTS, P, S
elif args_present['U']:
args = CP.PUmass_INPUTS, P, U
elif args_present['T']:
if args_present['Q']:
args = CP.QT_INPUTS, Q, T
elif args_present['H']:
args = CP.HmassT_INPUTS, H, T
elif args_present['S']:
args = CP.SmassT_INPUTS, S, T
elif args_present['U']:
args = CP.TUmass_INPUTS, T, U
elif args_present['Q']:
if args_present['H']:
args = CP.HmassQ_INPUTS, H, Q
elif args_present['S']:
args = CP.QSmass_INPUTS, Q, S
elif args_present['U']:
raise ValueError('Invalid combination: Q and U')
elif args_present['H']:
if args_present['S']:
args = CP.HmassSmass_INPUTS, H, S
elif args_present['U']:
raise ValueError('Invalid combination: H and U')
elif args_present['S'] and args_present['U']:
args = CP.SmassUmass_INPUTS, S, U
else:
raise ValueError('Invalid combination')
try:
self.state.update(*args)
except ValueError:
pass
def Vf(self, S: float):
x = self.chi()
return 1 / (1 /
(1 + S * ((1 - x) / x) * self.rhog([self.state.T()])[0] /
self.rhol([self.state.T()])[0]))
def chi(self):
x = self.state.Q()
if x >= 0:
return x
if self.state.rhomass() > self.rhol([self.state.T()])[0]:
return 0
else:
return 1
def viscosityl(self, T):
b1 = 1.6089
b2 = 2.0439
b3 = 5.24
b4 = 0.0293423
theta = (self.state.T_critical() - b3) / (T - b3)
return b4 * np.exp(b1 * (theta - 1)**(1 / 3) + b2 *
(theta - 1)**(4 / 3))
def viscosityg(self, T):
b1 = 3.3281
b2 = -1.18237
b3 = -0.055155
Tr = T / self.state.T_critical()
return 0.001 * np.exp(b1 + b2 * (1 / Tr - 1)**(1 / 3) + b3 *
(1 / Tr - 1)**(4 / 3)) if (
self.Tmin < T < self.Tmax) else None
def viscosity(self):
T = self.state.T()
mug = self.viscosityg(T)
mul = self.viscosityl(T)
x = self.chi()
rho = self.state.rhomass()
rhol = self.rhol([T])[0]
rhog = self.rhog([T])[0]
return mug * x * rho / rhog + mul * (1 - x) * rho / rhol
def dvg_dP_saturation(self):
dP = 0.002
rho1 = PropsSI('D', 'P', self.state.p() + dP / 2, 'Q', 1, self.name)
rho0 = PropsSI('D', 'P', self.state.p() - dP / 2, 'Q', 1, self.name)
dvg = 1 / rho1 - 1 / rho0
return dvg / dP
def rhol(self, T: Iterable):
return [
PropsSI('D', 'T', t, 'Q', 0, self.eos + "::" + self.name)
for t in T
]
def rhog(self, T: Iterable):
return [
PropsSI('D', 'T', t, 'Q', 1, self.eos + "::" + self.name)
for t in T
]
def psat(self, T: Iterable):
return [
PropsSI('P', 'Q', 0.5, 'T', t, self.eos + "::" + self.name) if
(self.Tmin < t < self.Tmax) else None for t in T
]
def tsat(self, P: Iterable):
return [
PropsSI('T', 'P', P, 'Q', 0.5, self.eos + "::" + self.name) if
(self.pmin < p < self.pmax) else None for p in P
]
def hl(self, T: Iterable):
return [
PropsSI('H', 'T', t, 'Q', 0, self.eos + "::" + self.name) if
(self.Tmin < t < self.Tmax) else None for t in T
]
def hg(self, T: Iterable):
return [
PropsSI('H', 'T', t, 'Q', 1, self.eos + "::" + self.name) if
(self.Tmin < t < self.Tmax) else None for t in T
]
def sl(self, T: Iterable):
return [
PropsSI('S', 'T', t, 'Q', 0, self.eos + "::" + self.name) if
(self.Tmin < t < self.Tmax) else None for t in T
]
def sg(self, T: Iterable):
return [
PropsSI('S', 'T', t, 'Q', 1, self.eos + "::" + self.name) if
(self.Tmin < t < self.Tmax) else None for t in T
]
def __repr__(self):
return f'Fluid({self.name}: p={self.state.p() / 100000:.1f}bar, t={self.state.T():.1f}K)'
def __init__(self, symbol="N2", T=530.0, P=1000.0, child=1):
'''Init generic Fluid'''
self.symbol = symbol.upper()
# http://www.coolprop.org/_static/doxygen/html/class_cool_prop_1_1_abstract_state.html
AS = AbstractState("HEOS", self.symbol)
self.AS = AS
self.name = AS.name()
self.T = T
self.P = P
self.child = child
Tsi = TSI_fromEng(T)
Psi = PSI_fromEng(P)
AS.update(CP.PT_INPUTS, Psi, Tsi)
self.WtMol = AS.molar_mass() * 1000.0
self.Tc = Teng_fromSI(AS.T_critical())
self.Pc = Peng_fromSI(AS.p_critical())
self.Dc = Deng_fromSI(AS.rhomass_critical())
self.Ttriple = Teng_fromSI(AS.Ttriple())
try:
self.Tfreeze = Teng_fromSI(AS.T_freeze())
except:
self.Tfreeze = self.Ttriple
self.Tmin = Teng_fromSI(AS.Tmin())
self.Tmax = Teng_fromSI(AS.Tmax())
#self.Pmin = Peng_fromSI( AS.pmin() ) # missing from AbstractState
self.Pmax = Peng_fromSI(AS.pmax())
dcIdeal = self.Pc * self.WtMol / self.Tc / 10.729
self.Zc = dcIdeal / self.Dc
try:
TnbpSI = PropsSI("T", "P", 101325, "Q", Q_LIQUID, self.symbol)
self.good_nbp = True
self.Tnbp = Teng_fromSI(TnbpSI)
except:
#print('WARNING... "%s" failed Normal Boiling Point Calculation.'%self.symbol)
Ttriple = PropsSI(self.symbol, 'Ttriple')
#print(' Using Triple Point = %g degK as Tref'%Ttriple)
self.good_nbp = False
self.Tnbp = 'N/A'
if self.good_nbp:
if self.Tnbp < 536.67:
self.Tref = self.Tnbp # if NBP is low, use NBP as ref
else:
self.Tref = 536.67 # 536.67R = (SATP, Standard Ambient T P) = 77F, 25C
else:
self.Tref = Teng_fromSI(Ttriple) + 0.1
self.Pref = 14.7
#print( 'About to call setTP #1 with Tref, Pref=',self.Tref,self.Pref )
self.setTP(self.Tref, self.Pref)
#print( 'Back from call setTP')
self.Href = self.H
#print( 'About to call setTP #2')
self.setTP(T, P)
#print( 'Back from call setTP')
if child == 1:
self.dup = EC_Fluid(symbol=self.symbol,
T=self.T,
P=self.P,
child=0)
self.calcdFreezePt = 0
#fluid = "R407C"
#fluid = "R32[0.381109419953993]&R125[0.179558888662016]&R134A[0.439331691383991]"
#fluid = "R32[0.40]&R125[0.15]&R134A[0.45]"
fluids = ["R32", "R125", "R134a"]
moles = [0.381109419953993, 0.179558888662016, 0.439331691383991]
for backend in ["HEOS", "REFPROP"]:
state = AbstractState(backend, '&'.join(fluids))
state.set_mole_fractions(moles)
print(state.get_mole_fractions())
print(state.get_mass_fractions())
try:
print("Normal routines")
print(state.T_critical())
print(state.p_critical())
print(True)
except:
try:
print("Routine for all points")
states = state.all_critical_points()
for crit_state in states:
print(crit_state.T)
print(crit_state.p)
print(crit_state.stable)
except:
print("All failed")
print("\n-----------------\n")
def get_critical_state(
state: CoolProp.AbstractState) -> CoolProp.AbstractState:
"""Create a new state instance and update it with the critical point data
Parameters
----------
state : CoolProp.AbstractState
Returns
-------
CoolProp.AbstractState
"""
crit_state = CP.PyCriticalState()
crit_state.T = np.nan
crit_state.p = np.nan
crit_state.rhomolar = np.nan
crit_state.stable = False
try:
crit_state.T = state.T_critical()
crit_state.p = state.p_critical()
crit_state.rhomolar = state.rhomolar_critical()
crit_state.stable = True
except:
try:
for crit_state_tmp in state.all_critical_points():
if crit_state_tmp.stable and (crit_state_tmp.T > crit_state.T
or
not np.isfinite(crit_state.T)):
crit_state.T = crit_state_tmp.T
crit_state.p = crit_state_tmp.p
crit_state.rhomolar = crit_state_tmp.rhomolar
crit_state.stable = crit_state_tmp.stable
except:
raise ValueError("Could not calculate the critical point data.")
new_state = CoolProp.AbstractState(state.backend_name(),
'&'.join(state.fluid_names()))
masses = state.get_mass_fractions()
if len(masses) > 1:
new_state.set_mass_fractions(masses)
# Uses mass fraction to work with incompressible fluids
# try: new_state.build_phase_envelope("dummy")
# except: pass
# TODO: Remove this hack ASAP
# Avoid problems with https://github.com/CoolProp/CoolProp/issues/1962
BE = state.backend_name()
FL = '&'.join(state.fluid_names())
if BE == "HelmholtzEOSBackend" and FL == "Ammonia":
crit_state.T *= 1.001
msg = ""
if np.isfinite(crit_state.rhomolar) and np.isfinite(crit_state.T):
try:
new_state.specify_phase(CoolProp.iphase_critical_point)
new_state.update(CoolProp.DmolarT_INPUTS, crit_state.rhomolar,
crit_state.T)
return new_state
except Exception as e:
msg += str(e) + " - "
pass
try:
new_state.update(CoolProp.DmolarT_INPUTS, crit_state.rhomolar,
crit_state.T)
return new_state
except Exception as e:
msg += str(e) + " - "
pass
if np.isfinite(crit_state.p) and np.isfinite(crit_state.T):
try:
new_state.specify_phase(CoolProp.iphase_critical_point)
new_state.update(CoolProp.PT_INPUTS, crit_state.p, crit_state.T)
return new_state
except Exception as e:
msg += str(e) + " - "
pass
try:
new_state.update(CoolProp.PT_INPUTS, crit_state.p, crit_state.T)
return new_state
except Exception as e:
msg += str(e) + " - "
pass
raise ValueError("Could not calculate the critical point data. " + msg)
