def calculo(self):
ind1 = self.Comp1.currentIndex()
ind2 = self.Comp2.currentIndex()
if ind1 != ind2:
zi = arange(0.025, 1., 0.025)
id1 = self.indices[ind1]
id2 = self.indices[ind2]
x = [0]
y = [0]
for z in zi:
try:
fraccion = [0.] * len(self.indices)
fraccion[ind1] = z
fraccion[ind2] = 1 - z
mez = Mezcla(tipo=3,
fraccionMolar=fraccion,
caudalMasico=1.)
tb = mez.componente[0].Tb
corr = Corriente(T=tb, P=101325., mezcla=mez)
T = corr.eos._Dew_T()
corr = Corriente(T=T, P=101325., mezcla=mez)
while corr.Liquido.fraccion[0] == corr.Gas.fraccion[
0] and corr.T < corr.mezcla.componente[1].Tb:
corr = Corriente(T=corr.T - 0.1, P=101325., mezcla=mez)
x.append(corr.Liquido.fraccion[0])
y.append(corr.Gas.fraccion[0])
except:
pass
x.append(1)
y.append(1)
self.rellenar(x, y)
a7, b7 = 0.0705233, -0.044448
a8, b8 = 0.504087, 1.32245
a9, b9 = 0.0307452, 0.179433
a10, b10 = 0.0732828, 0.463492
a11, b11 = 0.006450, -0.022143
Bo = (a1 + b1 * cmp.f_acent) * cmp.Vc
Ao = (a2 + b2 * cmp.f_acent) * R_atml * cmp.Tc * cmp.Vc
Co = (a3 + b3 * cmp.f_acent) * R_atml * cmp.Tc**3 * cmp.Vc
gamma = (a4 + b4 * cmp.f_acent) * (cmp.Vc)**2
b = (a5 + b5 * cmp.f_acent) * (cmp.Vc)**2
a = (a6 + b6 * cmp.f_acent) * R_atml * cmp.Tc * (cmp.Vc)**2
alfa = (a7 + b7 * cmp.f_acent) * (cmp.Vc)**3
c = (a8 + b8 * cmp.f_acent) * R_atml * cmp.Tc**3 * (cmp.Vc)**2
Do = (a9 + b9 * cmp.f_acent) * R_atml * cmp.Tc**4 * cmp.Vc
d = (a10 + b10 * cmp.f_acent) * R_atml * cmp.Tc**2 * (cmp.Vc)**2
Eo = (a11 + b11 * cmp.f_acent *
exp(-3.8 * cmp.f_acent)) * R_atml * cmp.Tc**5 * cmp.Vc
return Ao, Bo, Co, Do, Eo, a, b, c, d, alfa, gamma
_all = [BWRS]
if __name__ == "__main__":
from lib.mezcla import Mezcla
mezcla = Mezcla(2,
ids=[10, 38, 22, 61],
caudalUnitarioMolar=[0.3, 0.5, 0.05, 0.15])
eq = BWRS(340, 101325, mezcla)
print(eq.x)
__doi__ = {
"autor": "Twu, C.H., Coon, J.E., Cunningham, J.R.",
"title": "A New Generalized Alpha Function for a Cubic Equation of "
"State Part 1. Peng-Robinson equation",
"ref": "Fluid Phase Equilibria 105 (1995) 49-59",
"doi": "10.1016/0378-3812(94)02601-v"
},
OmegaA = 0.457235528921
OmegaB = 0.0777960739039
def _alfa(self, cmp, T):
Tr = T / cmp.Tc
alpha0 = Tr**(-0.171813) * exp(0.125283 * (1 - Tr**1.77634)) # Eq 11
alpha1 = Tr**(-0.607352) * exp(0.511614 * (1 - Tr**2.20517)) # Eq 12
# Eq 9
alpha = alpha0 + cmp.f_acent * (alpha1 - alpha0)
return 0, alpha
if __name__ == "__main__":
from lib.mezcla import Mezcla
mix = Mezcla(5, ids=[4], caudalMolar=1, fraccionMolar=[1])
eq = PRTwu(300, 9.9742e5, mix)
print('%0.0f %0.1f' % (eq.Vg.ccmol, eq.Vl.ccmol))
eq = PRTwu(300, 42.477e5, mix)
print('%0.1f' % (eq.Vl.ccmol))
# f_acent = (0.01100, 0.02200, -0.00200)
# x = (0.50000, 0.30000, 0.20000)
from lib.mezcla import Mezcla
from lib.compuestos import Componente
ch4 = Componente(2)
ch4.Tc, ch4.Pc, ch4.f_acent = 190.564, 4599200, 0.011
o2 = Componente(47)
o2.Tc, o2.Pc, o2.f_acent = 154.581, 5042800, 0.022
ar = Componente(98)
ar.Tc, ar.Pc, ar.f_acent = 150.687, 4863000, -0.002
mix = Mezcla(5,
customCmp=[ch4, o2, ar],
caudalMolar=1,
fraccionMolar=[0.5, 0.3, 0.2])
print(2222)
eq = PRMathiasCopeman(800, 34933409.8798343, mix)
print(eq._phir(800, 5000, eq.yi))
from lib.mezcla import Mezcla
from lib.compuestos import Componente
ch4 = Componente(2)
ch4.Tc, ch4.Pc, ch4.f_acent = 190.564, 4599200, 0.011
o2 = Componente(47)
o2.Tc, o2.Pc, o2.f_acent = 154.581, 5042800, 0.022
ar = Componente(98)
ar.Tc, ar.Pc, ar.f_acent = 150.687, 4863000, -0.002
mix = Mezcla(5,
customCmp=[ch4, o2, ar],
# mix = Mezcla(5, ids=[46, 2], caudalMolar=1, fraccionMolar=[0.2152, 0.7848])
# eq = SRK(144.26, 2.0684e6, mix)
# print(eq.rhoL.kmolm3)
# # Ejemplo 6.6, Wallas, pag 340
# mezcla = Mezcla(2, ids=[1, 2, 40, 41], caudalUnitarioMolar=[0.3177, 0.5894, 0.0715, 0.0214])
# P = unidades.Pressure(500, "psi")
# T = unidades.Temperature(120, "F")
# eq = SRK(T, P, mezcla)
# print(T)
# print(eq.x)
# print([x*(1-eq.x)*5597 for x in eq.xi])
# print([y*eq.x*5597 for y in eq.yi])
# print(eq.Ki)
# Example 6.6 wallas
P = unidades.Pressure(20, "atm")
mix = Mezcla(5, ids=[23, 5], caudalMolar=1, fraccionMolar=[0.607, 0.393])
eq1 = SRK(300, P, mix)
eq2 = SRK(400, P, mix)
print(eq1._Dew_T(P))
print(eq2._Dew_T(P))
# eq = SRK(500, P, mezcla)
# print(eq._Dew_T(P))
kw["a"] = self.tita
kw["dat"] = damt
kw["datt"] = damtt
kw["dattt"] = damttt
print(tau, delta, kw)
p = CubicHelmholtz(tau, delta, **kw)
p["tau"] = tau
p["delta"] = delta
print("fir: ", p["fir"])
print("fird: ", p["fird"] * delta)
print("firt: ", p["firt"] * tau)
print("firdd: ", p["firdd"] * delta**2)
print("firdt: ", p["firdt"] * delta * tau)
print("firtt: ", p["firtt"] * tau**2)
print("firddd: ", p["firddd"] * delta**3)
print("firddt: ", p["firddt"] * delta**2 * tau)
print("firdtt: ", p["firdtt"] * delta * tau**2)
print("firttt: ", p["firttt"] * tau**3)
print("P", (1 + delta * p["fird"]) * R * self.T * self.rho - self.P)
self.fir = p
if __name__ == "__main__":
from lib.mezcla import Mezcla
mix = Mezcla(tipo=5,
caudalMolar=1,
ids=[2, 47, 98],
fraccionMolar=[0.5, 0.3, 0.2])
eq = PRMathiasCopeman(800, 34937532, mix)
self.Bi = [bi*self.P/R/self.T for bi in self.bi]
self.Ai = [ai*self.P/(R*self.T)**2 for ai in self.ai]
def _lib(self, cmp, T):
a0 = 0.42747*R**2*cmp.Tc**2/cmp.Pc
alfa = (self.T/cmp.Tc)**-0.5
b = 0.08664*R*cmp.Tc/cmp.Pc
return a0*alfa, b
def _GEOS(self, xi):
am, bm = self._mixture("RK", xi, [self.ai, self.bi])
delta = bm
epsilon = 0
return am, bm, delta, epsilon
if __name__ == "__main__":
from lib.mezcla import Mezcla
from lib import unidades
# mix = Mezcla(2, ids=[10, 38, 22, 61],
# caudalUnitarioMolar=[0.3, 0.5, 0.05, 0.15])
# eq = RK(340, 101325, mix)
mezcla = Mezcla(2, ids=[1, 2, 40, 41],
caudalUnitarioMolar=[0.31767, 0.58942, 0.07147, 0.02144])
P = unidades.Pressure(485, "psi")
T = unidades.Temperature(100, "F")
eq = RK(T, P, mezcla, flory=1)
def readStatefromJSON(self, state):
self._thermo = state["thermo"]
self._bool = state["bool"]
self.tipoTermodinamica = state["thermoType"]
self.T = unidades.Temperature(state["T"])
self.P = unidades.Pressure(state["P"])
self.x = unidades.Dimensionless(state["x"])
self.M = unidades.Dimensionless(state["M"])
self.Tc = unidades.Temperature(state["Tc"])
self.Pc = unidades.Pressure(state["Pc"])
self.h = unidades.Power(state["h"])
self.s = unidades.Entropy(state["s"])
self.rho = unidades.Density(state["rho"])
self.Q = unidades.VolFlow(state["Q"])
self.SG = unidades.Dimensionless(state["SG"])
self.tipoFlujo = state["fluxType"]
self.mezcla = Mezcla()
self.mezcla.readStatefromJSON(state["mezcla"])
self.ids = self.mezcla.ids
self.componente = self.mezcla.componente
self.fraccion = self.mezcla.fraccion
self.caudalmasico = self.mezcla.caudalmasico
self.caudalmolar = self.mezcla.caudalmolar
self.fraccion_masica = self.mezcla.fraccion_masica
self.caudalunitariomasico = self.mezcla.caudalunitariomasico
self.caudalunitariomolar = self.mezcla.caudalunitariomolar
# Define the advanced thermo component if is necessary
if self._thermo == "freesteam":
self.cmp = freeSteam.Freesteam()
elif self._thermo == "iapws":
self.cmp = iapws97.IAPWS97()
elif self._thermo == "refprop":
self.cmp = refProp.RefProp(ids=self.ids)
elif self._thermo == "coolprop":
self.cmp = coolProp.CoolProp(ids=self.ids)
elif self._thermo == "meos":
eq = state["meos_eq"]
self.cmp = mEoS.__all__[mEoS.id_mEoS.index(self.ids[0])](eq=eq)
# Load advanced properties from global phase
if self._thermo != "eos":
self.cmp._readGlobalState(self, state)
# Load state for phases
if self._thermo in ["iapws", "freesteam"]:
self.Liquido = ThermoWater()
self.Gas = ThermoWater()
elif self._thermo in ["coolprop", "meos"]:
self.Liquido = ThermoAdvanced()
self.Gas = ThermoAdvanced()
elif self._thermo == "refprop":
self.Liquido = ThermoRefProp()
self.Gas = ThermoRefProp()
else:
self.Liquido = Mezcla()
self.Gas = Mezcla()
self.Liquido.readStatefromJSON(state["liquid"])
self.Gas.readStatefromJSON(state["gas"])
if state["liquid"]:
self.Liquido.sigma = unidades.Tension(state["liquid"]["sigma"])
class Corriente(config.Entity):
""" Class to model a stream object
Parameters:
-T: temperature, Kelvin
-P: Pressure, Pa
-x: quality
-caudalMasico: mass flow in kg/s (solid component excluded)
-caudalMolar: molar flow in kmol/s (solid component excluded)
-caudalVolumetrico: volumetric flow in m3/s
-caudalUnitarioMolar: Array with molar flows for each component
-caudalUnitarioMasico: Array with mass flows for each component
-fraccionMolar: Array with molar fractions for each component
-fraccionMasica: Array with mass fractions for each component
-mezcla: instance of class Mezcla to define all mixture variables
-solido: Instance of class Solid to define all solid variables
-caudalSolido: mass flow, in kg/h for solids component, array for each
solid component
-diametroMedio: mean diameter for particules, in micrometer
-distribucion_fraccion: Array with mass fraccion of solid particle
distribution
-distribucion_diametro: Array with particle diameter of solid particle
distribution, in micrometer
-notas: Description text for stream
Aditional parameter for define a corriente with different properties than
the project config, components for debug, thermodynamic method for add per
stream configuration:
-ids: Array with index of components
-solids: Array with index of solids components
-K: Name of method for K values calculation, ej: "SRK", "Lee-Kesler"
-alfa: Name of method for alpha calculation, available only for some K
methods
-mix: Mixing rule for calculate mix properties
-H: Name of method for enthalpy calculation, ej: "SRK", "Lee-Kesler"
-Cp_ideal: Bool to choose method for ideal cp:
0 - Ideal parameters
1 - DIPPR parameters
-MEoS: Use meos equation if is available
-iapws: Use iapws97 standard for water
-GERG: Use GERG-2008 equation if is available
-freesteam: Use freesteam external library for water
-coolProp: Use coolProp external library if is available
-refprop: Use refProp external library if is available
"""
kwargs = {"T": 0.0,
"P": 0.0,
"x": None,
"caudalMasico": 0.0,
"caudalVolumetrico": 0.0,
"caudalMolar": 0.0,
"caudalUnitarioMolar": [],
"caudalUnitarioMasico": [],
"fraccionMolar": [],
"fraccionMasica": [],
"mezcla": None,
"solido": None,
"caudalSolido": [],
"diametroMedio": 0.0,
"distribucion_fraccion": [],
"distribucion_diametro": [],
"ids": [],
"solids": None,
"K": "",
"alfa": "",
"mix": "",
"H": "",
"Cp_ideal": None,
"MEoS": None,
"iapws": None,
"GERG": None,
"freesteam": None,
"coolProp": None,
"refprop": None}
status = 0
msg = QApplication.translate("pychemqt", "Unknown variables")
kwargs_forbidden = ["entrada", "mezcla", "solido"]
solido = None
def __init__(self, **kwargs):
self.kwargs = Corriente.kwargs.copy()
self.__call__(**kwargs)
def __call__(self, **kwargs):
if kwargs.get("mezcla", None):
kwargs.update(kwargs["mezcla"].kwargs)
if kwargs.get("solido", None):
kwargs.update(kwargs["solido"].kwargs)
if kwargs.get("caudalUnitarioMasico", []):
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMolar"] = []
self.kwargs["fraccionMasica"] = []
self.kwargs["caudalMasico"] = None
self.kwargs["caudalVolumetrico"] = None
self.kwargs["caudalMolar"] = None
elif kwargs.get("caudalUnitarioMolar", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["fraccionMolar"] = []
self.kwargs["fraccionMasica"] = []
self.kwargs["caudalMasico"] = None
self.kwargs["caudalVolumetrico"] = None
self.kwargs["caudalMolar"] = None
elif kwargs.get("caudalMasico", None) and \
kwargs.get("fraccionMolar", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMasica"] = []
self.kwargs["caudalVolumetrico"] = None
self.kwargs["caudalMolar"] = None
elif kwargs.get("caudalMasico", None) and \
kwargs.get("fraccionMasica", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMolar"] = []
self.kwargs["caudalVolumetrico"] = None
self.kwargs["caudalMolar"] = None
elif kwargs.get("caudalMasico", None):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["caudalVolumetrico"] = None
self.kwargs["caudalMolar"] = None
elif kwargs.get("caudalMolar", None) and \
kwargs.get("fraccionMolar", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMasica"] = []
self.kwargs["caudalMasico"] = None
self.kwargs["caudalVolumetrico"] = None
elif kwargs.get("caudalMolar", None) and \
kwargs.get("fraccionMasica", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMolar"] = []
self.kwargs["caudalMasico"] = None
self.kwargs["caudalVolumetrico"] = None
elif kwargs.get("caudalMolar", None):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["caudalMasico"] = None
self.kwargs["caudalVolumetrico"] = None
elif kwargs.get("caudalVolumetrico", None) and \
kwargs.get("fraccionMolar", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMasica"] = []
self.kwargs["caudalMasico"] = None
elif kwargs.get("caudalVolumetrico") and \
kwargs.get("fraccionMasica", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMolar"] = []
self.kwargs["caudalMasico"] = None
elif kwargs.get("caudalVolumetrico", None):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["caudalMolar"] = []
self.kwargs["caudalMasico"] = None
elif kwargs.get("fraccionMasica", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMolar"] = []
elif kwargs.get("fraccionMolar", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMasica"] = []
elif kwargs.get("x", None) and self.kwargs["T"] and self.kwargs["P"]:
self.kwargs["T"] = 0.0
elif kwargs.get("T", 0.0) and self.kwargs["x"] and self.kwargs["P"]:
self.kwargs["P"] = 0.0
elif kwargs.get("P", 0.0) and self.kwargs["T"] and self.kwargs["x"]:
self.kwargs["x"] = None
self.kwargs.update(kwargs)
for key, value in list(self.kwargs.items()):
if value:
self._bool = True
break
logging.info('Calculate STREAM')
kw_new = {}
for key, value in list(kwargs.items()):
if self.__class__.kwargs[key] != value:
kw_new[key] = value
logging.debug('kwarg; %s' % kw_new)
if self.calculable:
statusmsg = (
QApplication.translate("pychemqt", "Underspecified"),
QApplication.translate("pychemqt", "Solved"),
QApplication.translate("pychemqt", "Ignored"),
QApplication.translate("pychemqt", "Warning"),
QApplication.translate("pychemqt", "Calculating..."),
QApplication.translate("pychemqt", "Error"))
status = statusmsg[self.status]
logging.debug('%s %s' % (status, self.msg))
QApplication.processEvents()
self.status = 1
self.calculo()
self.msg = ""
elif self.tipoFlujo:
if self.kwargs["mezcla"]:
self.mezcla = self.kwargs["mezcla"]
else:
self.mezcla = Mezcla(self.tipoFlujo, **self.kwargs)
elif self.tipoSolido:
if self.kwargs["solido"]:
self.solido = self.kwargs["solido"]
else:
self.solido = Solid(**self.kwargs)
if self.solido:
self.solido.RhoS(self.kwargs["T"])
@property
def calculable(self):
# Thermo definition
self.tipoTermodinamica = ""
if self.kwargs["T"] and self.kwargs["P"]:
self.tipoTermodinamica = "TP"
elif self.kwargs["T"] and self.kwargs["x"]:
self.tipoTermodinamica = "Tx"
elif self.kwargs["P"] and self.kwargs["x"]:
self.tipoTermodinamica = "Px"
# Mix definition
self.tipoFlujo = 0
if self.kwargs["caudalUnitarioMasico"]:
self.tipoFlujo = 1
elif self.kwargs["caudalUnitarioMolar"]:
self.tipoFlujo = 2
elif self.kwargs["caudalMasico"] and self.kwargs["fraccionMolar"]:
self.tipoFlujo = 3
elif self.kwargs["caudalMasico"] and self.kwargs["fraccionMasica"]:
self.tipoFlujo = 4
elif self.kwargs["caudalMolar"] and self.kwargs["fraccionMolar"]:
self.tipoFlujo = 5
elif self.kwargs["caudalMolar"] and self.kwargs["fraccionMasica"]:
self.tipoFlujo = 6
elif self.kwargs["caudalVolumetrico"] and self.kwargs["fraccionMolar"]:
self.kwargs["caudalMolar"] = 1
self.tipoFlujo = 5
elif self.kwargs["caudalVolumetrico"] and \
self.kwargs["fraccionMasica"]:
self.kwargs["caudalMolar"] = 1
self.tipoFlujo = 6
elif self.kwargs["mezcla"]:
self.tipoFlujo = 7
# Solid definition
self.tipoSolido = 0
if sum(self.kwargs["caudalSolido"]) > 0:
if self.kwargs["distribucion_fraccion"] and \
self.kwargs["distribucion_diametro"]:
self.tipoSolido = 2
elif self.kwargs["diametroMedio"]:
self.tipoSolido = 1
if self.kwargs["solido"]:
self.tipoSolido = self.kwargs["solido"].status
return self.tipoTermodinamica and self.tipoFlujo
def calculo(self):
Config = config.getMainWindowConfig()
if self.kwargs["mezcla"]:
self.mezcla = self.kwargs["mezcla"]
else:
self.mezcla = Mezcla(self.tipoFlujo, **self.kwargs)
self.ids = self.mezcla.ids
self.componente = self.mezcla.componente
self.fraccion = self.mezcla.fraccion
self.caudalmasico = self.mezcla.caudalmasico
self.caudalmolar = self.mezcla.caudalmolar
self.fraccion_masica = self.mezcla.fraccion_masica
self.caudalunitariomasico = self.mezcla.caudalunitariomasico
self.caudalunitariomolar = self.mezcla.caudalunitariomolar
T = unidades.Temperature(self.kwargs.get("T", None))
P = unidades.Pressure(self.kwargs.get("P", None))
x = self.kwargs.get("x", None)
self._method()
setData = True
if self._thermo == "freesteam":
compuesto = freeSteam.Freesteam(**self.kwargs)
elif self._thermo == "iapws":
compuesto = iapws97.IAPWS97(**self.kwargs)
elif self._thermo == "refprop":
if not self.kwargs["ids"]:
self.kwargs["ids"] = self.ids
# Avoid overwrite refprop H parameter
kwargs = self.kwargs.copy()
del kwargs["H"]
compuesto = refProp.RefProp(**kwargs)
elif self._thermo == "gerg":
ids = []
for id in self.ids:
ids.append(gerg.id_GERG.index(id))
kwargs = self.kwargs
kwargs["mezcla"] = self.mezcla
compuesto = gerg.GERG(componente=ids, fraccion=self.fraccion, **kwargs)
elif self._thermo == "coolprop":
if not self.kwargs["ids"]:
self.kwargs["ids"] = self.ids
compuesto = coolProp.CoolProp(**self.kwargs)
elif self._thermo == "meos":
if self.tipoTermodinamica == "TP":
compuesto = mEoS.__all__[mEoS.id_mEoS.index(self.ids[0])](T=T, P=P)
elif self.tipoTermodinamica == "Tx":
compuesto = mEoS.__all__[mEoS.id_mEoS.index(self.ids[0])](T=T, x=x)
elif self.tipoTermodinamica == "Px":
compuesto = mEoS.__all__[mEoS.id_mEoS.index(self.ids[0])](P=P, x=x)
elif self._thermo == "eos":
if self.kwargs["K"]:
index = K_name.index(self.kwargs["K"])
K = EoS.K[index]
print(K)
else:
K = EoS.K[Config.getint("Thermo","K")]
if self.kwargs["H"]:
index = H_name.index(self.kwargs["H"])
H = EoS.H[index]
print(H)
else:
H = EoS.H[Config.getint("Thermo","H")]
setData = False
self.M = unidades.Dimensionless(self.mezcla.M)
self.Tc = self.mezcla.Tc
self.Pc = self.mezcla.Pc
self.SG = unidades.Dimensionless(self.mezcla.SG)
if self.tipoTermodinamica == "TP":
self.T = unidades.Temperature(T)
self.P = unidades.Pressure(P)
eos = K(self.T, self.P.atm, self.mezcla)
self.eos = eos
self.x = unidades.Dimensionless(eos.x)
else:
self.x = unidades.Dimensionless(x)
# self.mezcla.recallZeros(eos.xi)
# self.mezcla.recallZeros(eos.yi)
# self.mezcla.recallZeros(eos.Ki, 1.)
if 0. < self.x < 1.:
self.Liquido = Mezcla(tipo=5, fraccionMolar=eos.xi, caudalMolar=self.caudalmolar*(1-self.x))
self.Gas = Mezcla(tipo=5, fraccionMolar=eos.yi, caudalMolar=self.caudalmolar*self.x)
elif self.x <= 0:
self.Liquido = self.mezcla
self.Gas = Mezcla()
else:
self.Liquido = Mezcla()
self.Gas = self.mezcla
self.Gas.Z = unidades.Dimensionless(float(eos.Z[0]))
self.Liquido.Z = unidades.Dimensionless(float(eos.Z[1]))
if H == K:
eosH = eos
else:
eosH = H(self.T, self.P.atm, self.mezcla)
self.H_exc = eosH.H_exc
self.Liquido.Q = unidades.VolFlow(0)
self.Gas.Q = unidades.VolFlow(0)
self.Liquido.h = unidades.Power(0)
self.Gas.h = unidades.Power(0)
if self.x < 1:
# There is liquid phase
Hl = (self.Liquido._Ho(self.T).Jg-self.Liquido.Hv_DIPPR(self.T).Jg)*self.Liquido.caudalmasico.gh
self.Liquido.h = unidades.Power(Hl-R*self.T/self.M*self.H_exc[1]*(1-self.x)*self.Liquido.caudalmasico.gh, "Jh")
self.Liquido.cp = self.Liquido.Cp_Liquido(T)
self.Liquido.rho = self.Liquido.RhoL(T, self.P)
self.Liquido.mu = self.Liquido.Mu_Liquido(T, self.P.atm)
self.Liquido.k = self.Liquido.ThCond_Liquido(T, self.P.atm, self.Liquido.rho)
self.Liquido.sigma = self.Liquido.Tension(T)
self.Liquido.Q = unidades.VolFlow(self.Liquido.caudalmasico/self.Liquido.rho)
self.Liquido.Prandt = self.Liquido.cp*self.Liquido.mu/self.Liquido.k
if self.x > 0:
# There is gas phase
Hg = self.Gas._Ho(self.T).Jg*self.Gas.caudalmasico.gh
self.Gas.h = unidades.Power(Hg-R*self.T/self.M*self.H_exc[0]*self.x*self.Gas.caudalmasico.gh, "Jh")
self.Gas.cp = self.Gas.Cp_Gas(T, self.P.atm)
self.Gas.rho = unidades.Density(self.P.atm/self.Gas.Z/R_atml/self.T*self.M, "gl")
self.Gas.rhoSd = unidades.Density(1./self.Gas.Z/R_atml/298.15*self.M, "gl")
self.Gas.mu = self.Gas.Mu_Gas(T, self.P.atm, self.Gas.rho)
self.Gas.k = self.Gas.ThCond_Gas(T, self.P.atm, self.Gas.rho)
self.Gas.Q = unidades.VolFlow(self.Gas.caudalmasico/self.Gas.rho)
self.Gas.Prandt = self.Gas.cp*self.Gas.mu/self.Gas.k
self.Q = unidades.VolFlow(self.Liquido.Q+self.Gas.Q)
self.h = unidades.Power(self.Liquido.h+self.Gas.h)
self.Molaridad = [caudal/self.Q.m3h for caudal in self.caudalunitariomolar]
# TODO:
self.cp_cv = 0.5
self.cp_cv_ideal = 0.5
self.s = 0
self.rho = 0
if setData:
# Asignación de valores comun
self.cmp = compuesto
self.T = compuesto.T
self.P = compuesto.P
self.x = compuesto.x
self.M = unidades.Dimensionless(compuesto.M)
self.Tc = compuesto.Tc
self.Pc = compuesto.Pc
self.h = unidades.Power(compuesto.h*self.caudalmasico)
self.s = unidades.Entropy(compuesto.s*self.caudalmasico)
self.rho = compuesto.rho
self.Q = unidades.VolFlow(compuesto.v*self.caudalmasico)
# harcoded thermo derived global properties in corriente for avoid
# losing data
if self._thermo != "eos":
compuesto._fillCorriente(self)
# if self.__class__!=mEoS.H2O:
# agua=mEoS.H2O(T=self.T, P=self.P.MPa)
# self.SG=unidades.Dimensionless(self.rho/agua.rho)
# else:
self.SG = unidades.Dimensionless(1.)
self.Liquido = compuesto.Liquido
self.Gas = compuesto.Gas
# if self.x<1: #Fase líquida
# self.Liquido=compuesto.Liquido
# if self.x>0: #Fase gaseosa
# self.Gas=compuesto
# elif 0<self.x<1: #Ambas fases
# self.Liquido=compuesto.Liquido
# self.Gas=compuesto.Gas
# self.QLiquido=self.Q*(1-self.x)
# elif self.x==1: #Fase gaseosa
# self.Gas=compuesto
# self.Liquido=None
# self.QLiquido=unidades.VolFlow(0)
# self.Gas.rhoSd=unidades.Density(1./R_atml/298.15*self.M, "gl")
if self.x > 0:
self.Gas.Q = unidades.VolFlow(self.Q*(1-self.x))
self.Gas.caudalmasico = unidades.MassFlow(self.caudalmasico*self.x)
self.Gas.caudalmolar = unidades.MolarFlow(self.caudalmolar*self.x)
kw = mix_molarflow_molarfraction(
self.Gas.caudalmolar, self.Gas.fraccion, self.componente)
self.Gas.caudalunitariomasico = [
unidades.MassFlow(f) for f in kw["unitMassFlow"]]
self.Gas.caudalunitariomolar = [
unidades.MolarFlow(f) for f in kw["unitMolarFlow"]]
self.Gas.ids = self.ids
if self.x < 1:
self.Liquido.Q = unidades.VolFlow(self.Q*(1-self.x))
self.Liquido.caudalmasico = unidades.MassFlow(self.caudalmasico*(1-self.x))
self.Liquido.caudalmolar = unidades.MolarFlow(self.caudalmolar*(1-self.x))
kw = mix_molarflow_molarfraction(
self.Liquido.caudalmolar, self.Liquido.fraccion, self.componente)
self.Liquido.caudalunitariomasico = [
unidades.MassFlow(f) for f in kw["unitMassFlow"]]
self.Liquido.caudalunitariomolar = [
unidades.MolarFlow(f) for f in kw["unitMolarFlow"]]
self.Liquido.sigma = compuesto.sigma
self.Liquido.ids = self.ids
if Config.get("Components", "Solids"):
if self.kwargs["solido"]:
self.solido = self.kwargs["solido"]
else:
self.solido = Solid(**self.kwargs)
if self.solido.status:
self.solido.RhoS(T)
else:
self.solido = None
if self.kwargs["caudalVolumetrico"]:
self.kwargs["caudalMolar"] *= self.kwargs["caudalVolumetrico"]/self.Q
Q = self.kwargs["caudalVolumetrico"]
self.kwargs["caudalVolumetrico"] = None
self.calculo()
self.kwargs["caudalVolumetrico"] = Q
self.kwargs["caudalMolar"] = None
def _method(self):
"""Find the thermodynamic method to use"""
Config = config.getMainWindowConfig()
# MEoS availability,
if self.kwargs["MEoS"] is not None:
_meos = self.kwargs["MEoS"]
else:
_meos = Config.getboolean("Thermo", "MEoS")
mEoS_available = self.ids[0] in mEoS.id_mEoS
MEoS = _meos and len(self.ids) == 1 and mEoS_available
# iapws availability
if self.kwargs["iapws"] is not None:
_iapws = self.kwargs["iapws"]
else:
_iapws = Config.getboolean("Thermo", "iapws")
IAPWS = _iapws and len(self.ids) == 1 and self.ids[0] == 62
# freesteam availability
if self.kwargs["freesteam"] is not None:
_freesteam = self.kwargs["freesteam"]
else:
_freesteam = Config.getboolean("Thermo", "freesteam")
FREESTEAM = _freesteam and len(self.ids) == 1 and \
self.ids[0] == 62 and os.environ["freesteam"]
COOLPROP_available = True
GERG_available = True
REFPROP_available = True
for id in self.ids:
if id not in coolProp.__all__:
COOLPROP_available = False
if id not in refProp.__all__:
REFPROP_available = False
if id not in gerg.id_GERG:
GERG_available = False
# coolprop availability
if self.kwargs["coolProp"] is not None:
_coolprop = self.kwargs["coolProp"]
else:
_coolprop = Config.getboolean("Thermo", "coolprop")
COOLPROP = _coolprop and os.environ["CoolProp"] and COOLPROP_available
# refprop availability
if self.kwargs["refprop"] is not None:
_refprop = self.kwargs["refprop"]
else:
_refprop = Config.getboolean("Thermo", "refprop")
REFPROP = _refprop and os.environ["refprop"] and REFPROP_available
# GERG availability
if self.kwargs["GERG"] is not None:
_gerg = self.kwargs["GERG"]
else:
_gerg = Config.getboolean("Thermo", "GERG")
GERG = _gerg and GERG_available
# Final selection
if IAPWS and FREESTEAM:
self._thermo = "freesteam"
self._dependence = "freesteam"
elif IAPWS:
self._thermo = "iapws"
elif _meos and REFPROP:
self._thermo = "refprop"
self._dependence = "refprop"
elif _meos and COOLPROP:
self._thermo = "coolprop"
self._dependence = "CoolProp"
elif MEoS and GERG:
self._thermo = "gerg"
elif MEoS:
self._thermo = "meos"
else:
self._thermo = "eos"
def setSolid(self, solid):
self.solido = solid
@property
def psystream(self):
xw = self.fraccion_masica[self.ids.index(62)]
xa = self.fraccion_masica[self.ids.index(475)]
psystream = PsyStream(caudalMasico=self.caudalMasico, P=self.P,
tdb=self.T, w=xw/xa)
return psystream
def clone(self, **kwargs):
"""Create a new stream instance with change only kwags new values"""
old_kwargs = self.kwargs.copy()
if "split" in kwargs:
split = kwargs["split"]
del kwargs["split"]
if self.kwargs["caudalUnitarioMasico"]:
kwargs["caudalUnitarioMasico"] = []
for caudal in self.kwargs["caudalUnitarioMasico"]:
kwargs["caudalUnitarioMasico"].append(split*caudal)
if self.kwargs["caudalUnitarioMolar"]:
kwargs["caudalUnitarioMolar"] = []
for caudal in self.kwargs["caudalUnitarioMolar"]:
kwargs["caudalUnitarioMolar"].append(split*caudal)
if self.kwargs["caudalMasico"]:
kwargs["caudalMasico"] = split*self.kwargs["caudalMasico"]
if self.kwargs["caudalMolar"]:
kwargs["caudalMolar"] = split*self.kwargs["caudalMolar"]
if "x" in kwargs:
del old_kwargs["T"]
if "mezcla" in kwargs:
old_kwargs.update(kwargs["mezcla"].kwargs)
del kwargs["mezcla"]
old_kwargs.update(kwargs)
return Corriente(**old_kwargs)
def __repr__(self):
if self.status:
return "Corriente at %0.2fK and %0.2fatm" % (self.T, self.P.atm)
else:
return "%s empty" % (self.__class__)
def txt(self):
txt = str(self.notasPlain)+os.linesep+os.linesep
txt += "#---------------"
txt += QApplication.translate("pychemqt", "Input properties")
txt += "-----------------#"+os.linesep
for key, value in list(self.kwargs.items()):
if value:
txt += key+": "+str(value)+os.linesep
if self.calculable:
txt += os.linesep + "#---------------"
txt += QApplication.translate("pychemqt", "Global stream")
txt += "-------------------#"+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Temperature"),
self.T.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Pressure"),
self.P.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Vapor Fraction"),
self.x.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Molar Flow"),
self.caudalmasico.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Mass Flow"),
self.caudalmolar.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Volumetric Flow"),
self.Q.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Enthalpy"),
self.h.str)+os.linesep
txt += "%-25s\t%s" % ("Tc", self.Tc.str)+os.linesep
txt += "%-25s\t%s" % ("Pc", self.Pc.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "SG, water=1"),
self.SG.str)+os.linesep
txt += os.linesep+"%-25s\t%s" % (
QApplication.translate("pychemqt", "Molecular weight"),
self.M.str)+os.linesep
txt += "#"+QApplication.translate("pychemqt", "Molar Composition")
txt += os.linesep
for cmp, xi in zip(self.componente, self.fraccion):
txt += "%-25s\t %0.4f" % (cmp.nombre, xi)+os.linesep
if self.x > 0:
txt += os.linesep+"#---------------"
txt += QApplication.translate("pychemqt", "Vapor Only")
txt += "--------------------#"+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Molar Flow"),
self.Gas.caudalmasico.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Mass Flow"),
self.Gas.caudalmolar.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Volumetric Flow"),
self.Gas.Q.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Molecular weight"),
self.Gas.M.str)+os.linesep
txt += os.linesep+"#"
txt += QApplication.translate("pychemqt", "Molar Composition")
txt += os.linesep
for cmp, xi in zip(self.componente, self.Gas.fraccion):
txt += "%-25s\t %0.4f" % (cmp.nombre, xi)+os.linesep
txt += os.linesep+"%-25s\t%s" % (
QApplication.translate("pychemqt", "Density"),
self.Gas.rho.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Compresibility"),
self.Gas.Z.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Enthalpy"),
self.Gas.h.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Heat Capacity"),
self.Gas.cp.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Viscosity"),
self.Gas.mu.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Thermal conductivity"),
self.Gas.k.str)+os.linesep
if self.x < 1:
txt += os.linesep+"#---------------"
txt += QApplication.translate("pychemqt", "Liquid Only")
txt += "-------------------#"+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Molar Flow"),
self.Liquido.caudalmasico.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Mass Flow"),
self.Liquido.caudalmolar.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Volumetric Flow"),
self.Liquido.Q.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Molecular weight"),
self.Liquido.M.str)+os.linesep
txt += os.linesep+"#"
txt += QApplication.translate("pychemqt", "Molar Composition")
txt += os.linesep
for cmp, xi in zip(self.componente, self.Liquido.fraccion):
txt += "%-25s\t %0.4f" % (cmp.nombre, xi)+os.linesep
txt += os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Density"),
self.Liquido.rho.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Compresibility"),
self.Liquido.Z.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Enthalpy"),
self.Liquido.h.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Heat Capacity"),
self.Liquido.cp.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Viscosity"),
self.Liquido.mu.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Thermal Conductivity"),
self.Liquido.k.str)+os.linesep
txt += "%-25s\t%s" % (
QApplication.translate("pychemqt", "Surface Tension"),
self.Liquido.sigma.str)+os.linesep
else:
txt += os.linesep+"#---------------"
txt += QApplication.translate("pychemqt", "No Fluid Stream")
txt += "-------------------#"+os.linesep
if self.solido.status:
txt += os.linesep+"#---------------"
txt += QApplication.translate("pychemqt", "Solid")
txt += "-------------------#"+os.linesep
for cmp, G in zip(self.solido.componente, self.solido.caudalUnitario):
txt += "%-25s\t%s" % (cmp.nombre, G.str)+os.linesep
txt += os.linesep
txt += "%-25s\t%s" % (QApplication.translate("pychemqt", "Density"),
self.solido.rho.str)+os.linesep
txt += "%-25s\t%s" % (QApplication.translate("pychemqt", "Mean Diameter"),
self.solido.diametro_medio.str)+os.linesep
if self.solido.diametros:
txt += os.linesep + "#"
txt += QApplication.translate("pychemqt",
"Particle Size Distribution")
txt += os.linesep
txt += "%s, %s \t%s" % (QApplication.translate("pychemqt", "Diameter"), unidades.Length.text("ParticleDiameter"), QApplication.translate("pychemqt", "Fraction"))+os.linesep
for di, xi in zip(self.solido.diametros, self.solido.fracciones):
txt += "%10.4f\t%0.4f\t" % (di.config("ParticleDiameter"),
xi)+os.linesep
# Add advanced properties to report
if self.calculable and self._thermo != "eos":
doc = self._doc()
if 0 < self.x < 1:
param = "%-40s\t%-20s\t%-20s"
else:
param = "%-40s\t%s"
if self.x == 0:
txtphases = "%60s" % QApplication.translate("pychemqt", "Liquid")+os.linesep
phases = [self.Liquido]
elif self.x == 1:
txtphases = "%60s" % QApplication.translate("pychemqt", "Gas")+os.linesep
phases = [self.Gas]
else:
txtphases = "%60s\t%20s" % (QApplication.translate("pychemqt", "Liquid"),
QApplication.translate("pychemqt", "Gas"))+os.linesep
phases = [self.Liquido, self.Gas]
complejos = ""
data = self.cmp.properties()
for propiedad, key, unit in data:
if key == "sigma":
if self.x < 1:
complejos += "%-40s\t%s" % (propiedad, self.Liquido.sigma.str)
complejos += os.linesep
elif key in ["f", "fi"]:
complejos += propiedad + os.linesep
for i, cmp in enumerate(self.componente):
values = [" " + cmp.nombre]
for phase in phases:
values.append(phase.__getattribute__(key)[i].str)
complejos += param % tuple(values) + os.linesep
elif key in ["K", "csat", "dpdt_sat", "cv2p", "chempot"]:
complejos += propiedad + os.linesep
for i, cmp in enumerate(self.componente):
values = [" " + cmp.nombre]
values.append(self.__getattribute__(key)[i].str)
complejos += "%-40s\t%s" % tuple(values)
complejos += os.linesep
elif key in self.Gas.__dict__ or key in self.Liquido.__dict__:
values = [propiedad]
for phase in phases:
values.append(phase.__getattribute__(key).str)
complejos += param % tuple(values) + os.linesep
else:
complejos += "%-40s\t%s" % (propiedad, self.__getattribute__(key).str)
complejos += os.linesep
txt += doc + os.linesep + txtphases + complejos
return txt
def _doc(self):
"""Return a text repr of class with all properties"""
if self._thermo == "meos":
title = QApplication.translate("pychemqt", "Advanced MEoS properties")
doc_param = [self.cmp._constants["__doi__"]]
else:
title = QApplication.translate("pychemqt", "Advanced thermo properties")
doc_param = self.cmp.__doi__
doc = ""
for doi in doc_param:
doc += doi["autor"] + "; " + doi["title"] + "; " + doi["ref"]
doc += os.linesep
txt = os.linesep + os.linesep + "#---------------"
txt += title + "-------------------#" + os.linesep
txt += doc
return txt
@classmethod
def propertiesNames(cls):
list = [
(QApplication.translate("pychemqt", "Temperature"), "T", unidades.Temperature),
(QApplication.translate("pychemqt", "Pressure"), "P", unidades.Pressure),
(QApplication.translate("pychemqt", "Vapor Fraction"), "x", unidades.Dimensionless),
(QApplication.translate("pychemqt", "Molar Flow"), "caudalmolar", unidades.MolarFlow),
(QApplication.translate("pychemqt", "Mass Flow"), "caudalmasico", unidades.MassFlow),
(QApplication.translate("pychemqt", "Volumetric Flow"), "Q", unidades.VolFlow),
(QApplication.translate("pychemqt", "Enthalpy"), "h", unidades.Enthalpy),
(QApplication.translate("pychemqt", "Critic Temperature"), "Tc", unidades.Temperature),
(QApplication.translate("pychemqt", "Critic Pressure"), "Pc", unidades.Pressure),
(QApplication.translate("pychemqt", "SG, water=1"), "SG", unidades.Dimensionless),
(QApplication.translate("pychemqt", "Molecular weight"), "M", unidades.Dimensionless),
(QApplication.translate("pychemqt", "Molar Composition"), "fraccion", unidades.Dimensionless),
(QApplication.translate("pychemqt", "Mass Composition"), "fraccion_masica", unidades.Dimensionless),
(QApplication.translate("pychemqt", "Molar Component Flow"), "caudalunitariomolar", unidades.MolarFlow),
(QApplication.translate("pychemqt", "Mass Component Flow"), "caudalunitariomasico", unidades.MassFlow),
(QApplication.translate("pychemqt", "Notes"), "notasPlain", str)]
return list
def propertiesListTitle(self, index):
"""Define los titulos para los popup de listas"""
lista = [comp.nombre for comp in self.componente]
return lista
def writeToJSON(self, data):
"""Read entity from file"""
config.Entity.writeToJSON(self, data)
# Solid state, can be defined without a thermo status
solid = {}
if self.solido is not None:
self.solido.writeStatetoJSON(solid)
data["solid"] = solid
def writeStatetoJSON(self, state):
state["thermo"] = self._thermo
state["bool"] = self._bool
state["thermoType"] = self.tipoTermodinamica
state["T"] = self.T
state["P"] = self.P
state["x"] = self.x
state["M"] = self.M
state["Tc"] = self.Tc
state["Pc"] = self.Pc
state["h"] = self.h
state["s"] = self.s
state["rho"] = self.rho
state["Q"] = self.Q
state["SG"] = self.SG
if self._thermo != "eos":
self.cmp._writeGlobalState(self, state)
state["fluxType"] = self.tipoFlujo
self.mezcla.writeStatetoJSON(state)
self.Liquido.writeStatetoJSON(state, "liquid")
self.Gas.writeStatetoJSON(state, "gas")
if state["liquid"]:
state["liquid"]["sigma"] = self.Liquido.sigma
if self._thermo == "meos":
state["meos_eq"] = self.cmp.kwargs["eq"]
def readFromJSON(self, data):
"""Read entity from file"""
config.Entity.readFromJSON(self, data)
# Read solid
self.solido = Solid()
self.solido.readStatefromJSON(data["solid"])
def readStatefromJSON(self, state):
self._thermo = state["thermo"]
self._bool = state["bool"]
self.tipoTermodinamica = state["thermoType"]
self.T = unidades.Temperature(state["T"])
self.P = unidades.Pressure(state["P"])
self.x = unidades.Dimensionless(state["x"])
self.M = unidades.Dimensionless(state["M"])
self.Tc = unidades.Temperature(state["Tc"])
self.Pc = unidades.Pressure(state["Pc"])
self.h = unidades.Power(state["h"])
self.s = unidades.Entropy(state["s"])
self.rho = unidades.Density(state["rho"])
self.Q = unidades.VolFlow(state["Q"])
self.SG = unidades.Dimensionless(state["SG"])
self.tipoFlujo = state["fluxType"]
self.mezcla = Mezcla()
self.mezcla.readStatefromJSON(state["mezcla"])
self.ids = self.mezcla.ids
self.componente = self.mezcla.componente
self.fraccion = self.mezcla.fraccion
self.caudalmasico = self.mezcla.caudalmasico
self.caudalmolar = self.mezcla.caudalmolar
self.fraccion_masica = self.mezcla.fraccion_masica
self.caudalunitariomasico = self.mezcla.caudalunitariomasico
self.caudalunitariomolar = self.mezcla.caudalunitariomolar
# Define the advanced thermo component if is necessary
if self._thermo == "freesteam":
self.cmp = freeSteam.Freesteam()
elif self._thermo == "iapws":
self.cmp = iapws97.IAPWS97()
elif self._thermo == "refprop":
self.cmp = refProp.RefProp(ids=self.ids)
elif self._thermo == "coolprop":
self.cmp = coolProp.CoolProp(ids=self.ids)
elif self._thermo == "meos":
eq = state["meos_eq"]
self.cmp = mEoS.__all__[mEoS.id_mEoS.index(self.ids[0])](eq=eq)
# Load advanced properties from global phase
if self._thermo != "eos":
self.cmp._readGlobalState(self, state)
# Load state for phases
if self._thermo in ["iapws", "freesteam"]:
self.Liquido = ThermoWater()
self.Gas = ThermoWater()
elif self._thermo in ["coolprop", "meos"]:
self.Liquido = ThermoAdvanced()
self.Gas = ThermoAdvanced()
elif self._thermo == "refprop":
self.Liquido = ThermoRefProp()
self.Gas = ThermoRefProp()
else:
self.Liquido = Mezcla()
self.Gas = Mezcla()
self.Liquido.readStatefromJSON(state["liquid"])
self.Gas.readStatefromJSON(state["gas"])
if state["liquid"]:
self.Liquido.sigma = unidades.Tension(state["liquid"]["sigma"])
def calculo(self):
Config = config.getMainWindowConfig()
if self.kwargs["mezcla"]:
self.mezcla = self.kwargs["mezcla"]
else:
self.mezcla = Mezcla(self.tipoFlujo, **self.kwargs)
self.ids = self.mezcla.ids
self.componente = self.mezcla.componente
self.fraccion = self.mezcla.fraccion
self.caudalmasico = self.mezcla.caudalmasico
self.caudalmolar = self.mezcla.caudalmolar
self.fraccion_masica = self.mezcla.fraccion_masica
self.caudalunitariomasico = self.mezcla.caudalunitariomasico
self.caudalunitariomolar = self.mezcla.caudalunitariomolar
T = unidades.Temperature(self.kwargs.get("T", None))
P = unidades.Pressure(self.kwargs.get("P", None))
x = self.kwargs.get("x", None)
self._method()
setData = True
if self._thermo == "freesteam":
compuesto = freeSteam.Freesteam(**self.kwargs)
elif self._thermo == "iapws":
compuesto = iapws97.IAPWS97(**self.kwargs)
elif self._thermo == "refprop":
if not self.kwargs["ids"]:
self.kwargs["ids"] = self.ids
# Avoid overwrite refprop H parameter
kwargs = self.kwargs.copy()
del kwargs["H"]
compuesto = refProp.RefProp(**kwargs)
elif self._thermo == "gerg":
ids = []
for id in self.ids:
ids.append(gerg.id_GERG.index(id))
kwargs = self.kwargs
kwargs["mezcla"] = self.mezcla
compuesto = gerg.GERG(componente=ids, fraccion=self.fraccion, **kwargs)
elif self._thermo == "coolprop":
if not self.kwargs["ids"]:
self.kwargs["ids"] = self.ids
compuesto = coolProp.CoolProp(**self.kwargs)
elif self._thermo == "meos":
if self.tipoTermodinamica == "TP":
compuesto = mEoS.__all__[mEoS.id_mEoS.index(self.ids[0])](T=T, P=P)
elif self.tipoTermodinamica == "Tx":
compuesto = mEoS.__all__[mEoS.id_mEoS.index(self.ids[0])](T=T, x=x)
elif self.tipoTermodinamica == "Px":
compuesto = mEoS.__all__[mEoS.id_mEoS.index(self.ids[0])](P=P, x=x)
elif self._thermo == "eos":
if self.kwargs["K"]:
index = K_name.index(self.kwargs["K"])
K = EoS.K[index]
print(K)
else:
K = EoS.K[Config.getint("Thermo","K")]
if self.kwargs["H"]:
index = H_name.index(self.kwargs["H"])
H = EoS.H[index]
print(H)
else:
H = EoS.H[Config.getint("Thermo","H")]
setData = False
self.M = unidades.Dimensionless(self.mezcla.M)
self.Tc = self.mezcla.Tc
self.Pc = self.mezcla.Pc
self.SG = unidades.Dimensionless(self.mezcla.SG)
if self.tipoTermodinamica == "TP":
self.T = unidades.Temperature(T)
self.P = unidades.Pressure(P)
eos = K(self.T, self.P.atm, self.mezcla)
self.eos = eos
self.x = unidades.Dimensionless(eos.x)
else:
self.x = unidades.Dimensionless(x)
# self.mezcla.recallZeros(eos.xi)
# self.mezcla.recallZeros(eos.yi)
# self.mezcla.recallZeros(eos.Ki, 1.)
if 0. < self.x < 1.:
self.Liquido = Mezcla(tipo=5, fraccionMolar=eos.xi, caudalMolar=self.caudalmolar*(1-self.x))
self.Gas = Mezcla(tipo=5, fraccionMolar=eos.yi, caudalMolar=self.caudalmolar*self.x)
elif self.x <= 0:
self.Liquido = self.mezcla
self.Gas = Mezcla()
else:
self.Liquido = Mezcla()
self.Gas = self.mezcla
self.Gas.Z = unidades.Dimensionless(float(eos.Z[0]))
self.Liquido.Z = unidades.Dimensionless(float(eos.Z[1]))
if H == K:
eosH = eos
else:
eosH = H(self.T, self.P.atm, self.mezcla)
self.H_exc = eosH.H_exc
self.Liquido.Q = unidades.VolFlow(0)
self.Gas.Q = unidades.VolFlow(0)
self.Liquido.h = unidades.Power(0)
self.Gas.h = unidades.Power(0)
if self.x < 1:
# There is liquid phase
Hl = (self.Liquido._Ho(self.T).Jg-self.Liquido.Hv_DIPPR(self.T).Jg)*self.Liquido.caudalmasico.gh
self.Liquido.h = unidades.Power(Hl-R*self.T/self.M*self.H_exc[1]*(1-self.x)*self.Liquido.caudalmasico.gh, "Jh")
self.Liquido.cp = self.Liquido.Cp_Liquido(T)
self.Liquido.rho = self.Liquido.RhoL(T, self.P)
self.Liquido.mu = self.Liquido.Mu_Liquido(T, self.P.atm)
self.Liquido.k = self.Liquido.ThCond_Liquido(T, self.P.atm, self.Liquido.rho)
self.Liquido.sigma = self.Liquido.Tension(T)
self.Liquido.Q = unidades.VolFlow(self.Liquido.caudalmasico/self.Liquido.rho)
self.Liquido.Prandt = self.Liquido.cp*self.Liquido.mu/self.Liquido.k
if self.x > 0:
# There is gas phase
Hg = self.Gas._Ho(self.T).Jg*self.Gas.caudalmasico.gh
self.Gas.h = unidades.Power(Hg-R*self.T/self.M*self.H_exc[0]*self.x*self.Gas.caudalmasico.gh, "Jh")
self.Gas.cp = self.Gas.Cp_Gas(T, self.P.atm)
self.Gas.rho = unidades.Density(self.P.atm/self.Gas.Z/R_atml/self.T*self.M, "gl")
self.Gas.rhoSd = unidades.Density(1./self.Gas.Z/R_atml/298.15*self.M, "gl")
self.Gas.mu = self.Gas.Mu_Gas(T, self.P.atm, self.Gas.rho)
self.Gas.k = self.Gas.ThCond_Gas(T, self.P.atm, self.Gas.rho)
self.Gas.Q = unidades.VolFlow(self.Gas.caudalmasico/self.Gas.rho)
self.Gas.Prandt = self.Gas.cp*self.Gas.mu/self.Gas.k
self.Q = unidades.VolFlow(self.Liquido.Q+self.Gas.Q)
self.h = unidades.Power(self.Liquido.h+self.Gas.h)
self.Molaridad = [caudal/self.Q.m3h for caudal in self.caudalunitariomolar]
# TODO:
self.cp_cv = 0.5
self.cp_cv_ideal = 0.5
self.s = 0
self.rho = 0
if setData:
# Asignación de valores comun
self.cmp = compuesto
self.T = compuesto.T
self.P = compuesto.P
self.x = compuesto.x
self.M = unidades.Dimensionless(compuesto.M)
self.Tc = compuesto.Tc
self.Pc = compuesto.Pc
self.h = unidades.Power(compuesto.h*self.caudalmasico)
self.s = unidades.Entropy(compuesto.s*self.caudalmasico)
self.rho = compuesto.rho
self.Q = unidades.VolFlow(compuesto.v*self.caudalmasico)
# harcoded thermo derived global properties in corriente for avoid
# losing data
if self._thermo != "eos":
compuesto._fillCorriente(self)
# if self.__class__!=mEoS.H2O:
# agua=mEoS.H2O(T=self.T, P=self.P.MPa)
# self.SG=unidades.Dimensionless(self.rho/agua.rho)
# else:
self.SG = unidades.Dimensionless(1.)
self.Liquido = compuesto.Liquido
self.Gas = compuesto.Gas
# if self.x<1: #Fase líquida
# self.Liquido=compuesto.Liquido
# if self.x>0: #Fase gaseosa
# self.Gas=compuesto
# elif 0<self.x<1: #Ambas fases
# self.Liquido=compuesto.Liquido
# self.Gas=compuesto.Gas
# self.QLiquido=self.Q*(1-self.x)
# elif self.x==1: #Fase gaseosa
# self.Gas=compuesto
# self.Liquido=None
# self.QLiquido=unidades.VolFlow(0)
# self.Gas.rhoSd=unidades.Density(1./R_atml/298.15*self.M, "gl")
if self.x > 0:
self.Gas.Q = unidades.VolFlow(self.Q*(1-self.x))
self.Gas.caudalmasico = unidades.MassFlow(self.caudalmasico*self.x)
self.Gas.caudalmolar = unidades.MolarFlow(self.caudalmolar*self.x)
kw = mix_molarflow_molarfraction(
self.Gas.caudalmolar, self.Gas.fraccion, self.componente)
self.Gas.caudalunitariomasico = [
unidades.MassFlow(f) for f in kw["unitMassFlow"]]
self.Gas.caudalunitariomolar = [
unidades.MolarFlow(f) for f in kw["unitMolarFlow"]]
self.Gas.ids = self.ids
if self.x < 1:
self.Liquido.Q = unidades.VolFlow(self.Q*(1-self.x))
self.Liquido.caudalmasico = unidades.MassFlow(self.caudalmasico*(1-self.x))
self.Liquido.caudalmolar = unidades.MolarFlow(self.caudalmolar*(1-self.x))
kw = mix_molarflow_molarfraction(
self.Liquido.caudalmolar, self.Liquido.fraccion, self.componente)
self.Liquido.caudalunitariomasico = [
unidades.MassFlow(f) for f in kw["unitMassFlow"]]
self.Liquido.caudalunitariomolar = [
unidades.MolarFlow(f) for f in kw["unitMolarFlow"]]
self.Liquido.sigma = compuesto.sigma
self.Liquido.ids = self.ids
if Config.get("Components", "Solids"):
if self.kwargs["solido"]:
self.solido = self.kwargs["solido"]
else:
self.solido = Solid(**self.kwargs)
if self.solido.status:
self.solido.RhoS(T)
else:
self.solido = None
if self.kwargs["caudalVolumetrico"]:
self.kwargs["caudalMolar"] *= self.kwargs["caudalVolumetrico"]/self.Q
Q = self.kwargs["caudalVolumetrico"]
self.kwargs["caudalVolumetrico"] = None
self.calculo()
self.kwargs["caudalVolumetrico"] = Q
self.kwargs["caudalMolar"] = None
def __call__(self, **kwargs):
if kwargs.get("mezcla", None):
kwargs.update(kwargs["mezcla"].kwargs)
if kwargs.get("solido", None):
kwargs.update(kwargs["solido"].kwargs)
if kwargs.get("caudalUnitarioMasico", []):
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMolar"] = []
self.kwargs["fraccionMasica"] = []
self.kwargs["caudalMasico"] = None
self.kwargs["caudalVolumetrico"] = None
self.kwargs["caudalMolar"] = None
elif kwargs.get("caudalUnitarioMolar", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["fraccionMolar"] = []
self.kwargs["fraccionMasica"] = []
self.kwargs["caudalMasico"] = None
self.kwargs["caudalVolumetrico"] = None
self.kwargs["caudalMolar"] = None
elif kwargs.get("caudalMasico", None) and \
kwargs.get("fraccionMolar", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMasica"] = []
self.kwargs["caudalVolumetrico"] = None
self.kwargs["caudalMolar"] = None
elif kwargs.get("caudalMasico", None) and \
kwargs.get("fraccionMasica", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMolar"] = []
self.kwargs["caudalVolumetrico"] = None
self.kwargs["caudalMolar"] = None
elif kwargs.get("caudalMasico", None):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["caudalVolumetrico"] = None
self.kwargs["caudalMolar"] = None
elif kwargs.get("caudalMolar", None) and \
kwargs.get("fraccionMolar", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMasica"] = []
self.kwargs["caudalMasico"] = None
self.kwargs["caudalVolumetrico"] = None
elif kwargs.get("caudalMolar", None) and \
kwargs.get("fraccionMasica", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMolar"] = []
self.kwargs["caudalMasico"] = None
self.kwargs["caudalVolumetrico"] = None
elif kwargs.get("caudalMolar", None):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["caudalMasico"] = None
self.kwargs["caudalVolumetrico"] = None
elif kwargs.get("caudalVolumetrico", None) and \
kwargs.get("fraccionMolar", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMasica"] = []
self.kwargs["caudalMasico"] = None
elif kwargs.get("caudalVolumetrico") and \
kwargs.get("fraccionMasica", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMolar"] = []
self.kwargs["caudalMasico"] = None
elif kwargs.get("caudalVolumetrico", None):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["caudalMolar"] = []
self.kwargs["caudalMasico"] = None
elif kwargs.get("fraccionMasica", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMolar"] = []
elif kwargs.get("fraccionMolar", []):
self.kwargs["caudalUnitarioMasico"] = []
self.kwargs["caudalUnitarioMolar"] = []
self.kwargs["fraccionMasica"] = []
elif kwargs.get("x", None) and self.kwargs["T"] and self.kwargs["P"]:
self.kwargs["T"] = 0.0
elif kwargs.get("T", 0.0) and self.kwargs["x"] and self.kwargs["P"]:
self.kwargs["P"] = 0.0
elif kwargs.get("P", 0.0) and self.kwargs["T"] and self.kwargs["x"]:
self.kwargs["x"] = None
self.kwargs.update(kwargs)
for key, value in list(self.kwargs.items()):
if value:
self._bool = True
break
logging.info('Calculate STREAM')
kw_new = {}
for key, value in list(kwargs.items()):
if self.__class__.kwargs[key] != value:
kw_new[key] = value
logging.debug('kwarg; %s' % kw_new)
if self.calculable:
statusmsg = (
QApplication.translate("pychemqt", "Underspecified"),
QApplication.translate("pychemqt", "Solved"),
QApplication.translate("pychemqt", "Ignored"),
QApplication.translate("pychemqt", "Warning"),
QApplication.translate("pychemqt", "Calculating..."),
QApplication.translate("pychemqt", "Error"))
status = statusmsg[self.status]
logging.debug('%s %s' % (status, self.msg))
QApplication.processEvents()
self.status = 1
self.calculo()
self.msg = ""
elif self.tipoFlujo:
if self.kwargs["mezcla"]:
self.mezcla = self.kwargs["mezcla"]
else:
self.mezcla = Mezcla(self.tipoFlujo, **self.kwargs)
elif self.tipoSolido:
if self.kwargs["solido"]:
self.solido = self.kwargs["solido"]
else:
self.solido = Solid(**self.kwargs)
if self.solido:
self.solido.RhoS(self.kwargs["T"])
# Eq 21
bptr = 1 - 0.1519*cmp.f_acent - 3.9462*cmp.f_acent**2 + \
7.0539*cmp.f_acent**3
bpt = bptr * bc
mb = bptr - 1 # Eq 17
apt = 29.7056 * bpt * R * Tpt # Eq 19
b = bc * (1 + mb * (1 - Tita)) # Eq 16
ma = (apt / ac)**0.5 - 1 # Eq 11
a = ac * (1 + ma * (1 - Tita**0.5))**2 # Eq 10
return a, b
if __name__ == "__main__":
from lib.mezcla import Mezcla
mix = Mezcla(5, ids=[4], caudalMolar=1, fraccionMolar=[1])
eq = Nasrifar(300, 9.9742e5, mix)
print('%0.1f' % (eq.Vl.ccmol))
eq = Nasrifar(300, 42.477e5, mix)
print('%0.1f' % (eq.Vg.ccmol))
mix = Mezcla(5,
ids=[46, 2],
caudalMolar=1,
fraccionMolar=[0.04752, 0.95248])
eq = Nasrifar(105, 5e7, mix)
print(eq.rhoL.kmolm3)