#!/usr/bin/python3
# -*- coding: utf-8 -*-
"""Pychemqt, Chemical Engineering Process simulator
Copyright (C) 2009-2025, Juan José Gómez Romera <jjgomera@gmail.com>
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program. If not, see <http://www.gnu.org/licenses/>."""
###############################################################################
# Module to define gas pressure equipments:
# - Compressor
# - Turbine
###############################################################################
from math import log, exp
import os
from tools.qt import translate
from scipy.constants import R
from scipy.optimize import fsolve
from lib.unidades import (DeltaT, DeltaP, Temperature, Pressure, MassFlow,
Power, Currency, Dimensionless)
from equipment.parents import equipment
[docs]
class Compressor(equipment):
"""Class to model a gas compressor
Parameters:
entrada: Instance of class Corriente to define the imput stream
metodo: Value to specified the calculate type
0 - Output pressure and efficiency
1 - Compression ratio and efficiency
2 - Real work and efficiency
3 - Output pressure and real work
4 - Compression ratio and real work
5 - Calculate flow ot stream to compress known output pressure,
work and efficiency
6 - Compressor curve (unimplemented)
termodicnamica: Index with the thermodynamic model
0 - Adiabatic
1 - Polytropic
2 - Isothermic
Pout: Output pressure
razon: Compression ratio
rendimeinto: Compressor efficiency
etapas: Compressor etapas
trabajo: Compressor work
Coste
f_install: instalation factor
base: Base index
actual: Current index
compresor: Compressor type
0 - Centrifugal compressor
1 - Reciprocating compressor
2 - Screw compressor
transmision: Transmission type
0 - Belt drive coupling
1 - Chain drive coupling
2 - Variable speed drive coupling
motor: Motor type
0 - Open drip-proof
1 - Totally enclosed, fan-cooled
2 - Explosion-proof
rpm
0 - 3600 rpm
1 - 1800 rpm
2 - 1200 rpm
>>> from lib.corriente import Corriente
>>> corriente=Corriente(T=400, P=101325, caudalMasico=0.1, \
fraccionMasica=[1., 0, 0, 0])
>>> compresor=Compressor(entrada=corriente, metodo=1, termodinamica=0, \
razon=3, rendimiento=0.75, etapas=1)
>>> print(compresor.power.kW)
31.1007420914
>>> compresor(compresor=2, transmision=1, motor=0, rpm=1)
>>> print(compresor.C_inst)
60464.196881
"""
title = translate("equipment", "Compressor")
help = ""
kwargs = {"entrada": None,
"metodo": 0,
"termodinamica": 0,
"Pout": 0.0,
"razon": 0.0,
"rendimiento": 0.0,
"etapas": 0,
"trabajo": 0.0,
"f_install": 1.3,
"Base_index": 0.0,
"Current_index": 0.0,
"compresor": 0,
"transmision": 0,
"motor": 0,
"rpm": 0}
kwargsInput = ("entrada", )
kwargsValue = ("Pout", "razon", "rendimiento", "etapas", "trabajo")
kwargsList = ("metodo", "termodinamica", "compresor", "transmision",
"motor", "rpm")
calculateValue = ("power", "cp_cv", "razonCalculada",
"rendimientoCalculado")
calculateCostos = ("C_comp", "C_motor", "C_trans", "C_adq", "C_inst")
indiceCostos = 7
TEXT_METODO = [
translate("equipment", "Specify out pressure and efficiency"),
translate("equipment", "Specify actual power and efficiency"),
translate("equipment", "Specify out pressure and actual power"),
translate("equipment", "Specify pressure ratio and actual power"),
translate("equipment", "Calculate input flowrate")]
TEXT_TERMODINAMICA = [translate("equipment", "Adiabatic"),
translate("equipment", "Polytropic"),
translate("equipment", "Isothermic")]
TEXT_COMPRESOR = [
translate("equipment", "Centrifugal compressor"),
translate("equipment", "Reciprocating compressor"),
translate("equipment", "Screw compressor")]
TEXT_TRANSMISION = [
translate("equipment", "Belt drive coupling"),
translate("equipment", "Chain drive coupling"),
translate("equipment", "Variable speed drive coupling")]
TEXT_MOTOR = [
translate("equipment", "Open drip-proof"),
translate("equipment", "Totally enclosed, fan-cooled"),
translate("equipment", "Explosion-proof")]
TEXT_RPM = ["3600 RPM", "1800 RPM", "1200 RPM"]
@property
def isCalculable(self):
if self.kwargs["f_install"] and self.kwargs["Base_index"] and \
self.kwargs["Current_index"]:
self.statusCoste = True
else:
self.statusCoste = False
if not self.kwargs["etapas"]:
self.msg = translate("equipment", "undefined variables")
self.status = 0
if not self.kwargs["entrada"]:
if self.kwargs["metodo"] == 5:
if (self.kwargs["razon"] or self.kwargs["Pout"]) and \
self.kwargs["trabajo"] and self.kwargs["rendimiento"]:
self.status = 1
self.msg = ""
return True
else:
self.msg = translate("equipment", "undefined variables")
self.status = 0
else:
self.msg = translate("equipment", "undefined input")
self.status = 0
else:
if self.kwargs["metodo"] == 0:
valores = self.kwargs["Pout"] and self.kwargs["rendimiento"]
elif self.kwargs["metodo"] == 1:
valores = self.kwargs["razon"] and self.kwargs["rendimiento"]
elif self.kwargs["metodo"] == 2:
valores = self.kwargs["trabajo"] and self.kwargs["rendimiento"]
elif self.kwargs["metodo"] == 3:
valores = self.kwargs["Pout"] and self.kwargs["trabajo"]
elif self.kwargs["metodo"] == 4:
valores = self.kwargs["razon"] and self.kwargs["trabajo"]
elif self.kwargs["metodo"] == 5:
valores = self.kwargs["razon"] and self.kwargs["rendimiento"]
if valores:
self.status = 1
self.msg = ""
return True
else:
self.msg = translate("equipment", "undefined variables")
self.status = 0
[docs]
def calculo(self):
self.entrada = self.kwargs["entrada"]
metodo = self.kwargs["metodo"]
self.Pout = Pressure(self.kwargs["Pout"])
razon = self.kwargs["razon"]
self.rendimientoCalculado = Dimensionless(self.kwargs["rendimiento"])
if self.kwargs["etapas"]:
self.etapas = self.kwargs["etapas"]
else:
self.etapas = 1.
self.power = Power(self.kwargs["trabajo"])
def f(Pout, rendimiento):
W_ideal = self.__Wideal(Pout)
power = W_ideal*self.entrada.caudalmasico.gs/rendimiento
return power
if metodo in [0, 3] or (metodo == 5 and self.Pout):
if self.etapas == 1:
razon = self.Pout.atm/self.entrada.P.atm
else:
razon = (self.Pout.atm/self.entrada.P.atm)**(1./self.etapas)
elif metodo in [1, 4] or (metodo == 5 and razon):
if self.etapas == 1:
self.Pout = Pressure(self.entrada.P*razon)
else:
self.Pout = Pressure(razon**self.etapas*self.entrada.P)
if metodo in [0, 1]:
Wid = self.__Wideal(self.Pout.atm)
power = Wid*self.entrada.caudalmasico.gs/self.rendimientoCalculado
self.power = Power(power*self.etapas)
elif metodo == 2:
def funcion(P):
return f(P, self.rendimientoCalculado)-self.power
self.Pout = Pressure(fsolve(funcion, self.entrada.P+1))
if self.etapas == 1:
razon = self.Pout/self.entrada.P
else:
razon = (self.Pout/self.entrada.P)**(1./self.etapas)
elif metodo in [3, 4]:
def funcion(rendimiento):
return f(self.Pout.atm, rendimiento)-self.power
self.rendimientoCalculado = Dimensionless(fsolve(funcion, 0.5))
elif metodo == 5:
Wideal = self.__Wideal(self.Pout.atm)
G = MassFlow(self.power*self.rendimientoCalculado/Wideal, "gs")
self.Tout = self.__Tout(Wideal)
self.entrada = self.entrada.clone(caudalMasico=G)
Wideal = self.__Wideal(self.Pout.atm)
self.Tout = Temperature(self.__Tout(Wideal))
self.salida = [self.entrada.clone(T=self.Tout, P=self.Pout)]
self.razonCalculada = Dimensionless(self.Pout/self.entrada.P)
self.deltaT = DeltaT(self.salida[0].T-self.entrada.T)
self.deltaP = DeltaP(self.salida[0].P-self.entrada.P)
self.cp_cv = self.entrada.Gas.cp_cv
self.Pin = self.entrada.P
self.Tin = self.entrada.T
def __Wideal(self, Pout):
"""Calculate the ideal work"""
if self.kwargs["termodinamica"] == 0:
cpv = self.entrada.Gas.cp_cv
else:
cpv = 1.40388
if self.kwargs["termodinamica"] == 2:
W_ideal = R*self.entrada.T/self.entrada.M*log(Pout/self.entrada.P)
else:
W_ideal = R*self.entrada.T*cpv/self.entrada.M/(cpv-1) *\
((Pout/self.entrada.P.atm)**((cpv-1)/cpv)-1)
return W_ideal
def __Tout(self, Wid):
"""Calculate the output temperature"""
if self.kwargs["termodinamica"] == 0:
cpv = self.entrada.Gas.cp_cv
else:
cpv = 1.40388
if self.kwargs["termodinamica"] == 2:
Tout = self.entrada.T
else:
eta = self.rendimientoCalculado
Tout = (self.Pout.atm/self.entrada.P.atm)**((cpv-1)/cpv) *\
self.entrada.T+(1-eta)/eta*(Wid/self.entrada.Gas.cv)
return Tout
[docs]
def coste(self):
HP = self.power.hp/self.etapas
CI = self.kwargs["Current_index"]
BI = self.kwargs["Base_index"]
if self.kwargs["compresor"] == 0:
# Centrifugal compressor
C = 6.49*HP**0.62*1000
elif self.kwargs["compresor"] == 1:
# Reciprocating compressor
C = 5.96*HP**0.61*1000
elif self.kwargs["compresor"] == 2:
# Screw compressor
C = 1.49*HP**0.71*1000
C_comp = self.etapas*C*CI/BI
if self.kwargs["motor"] == 0: # Open, drip-proof
if self.kwargs["rpm"] == 0 and HP <= 7.5:
a1, a2, a3 = 4.8314, 0.0966, 0.10960
elif self.kwargs["rpm"] == 0 and HP <= 250. and HP > 7.5:
a1, a2, a3 = 4.1514, 0.5347, 0.05252
elif self.kwargs["rpm"] == 0 and HP > 250.:
a1, a2, a3 = 4.2432, 1.03251, -0.03595
elif self.kwargs["rpm"] == 1 and HP <= 7.5:
a1, a2, a3 = 4.7075, -0.01511, 0.22888
elif self.kwargs["rpm"] == 1 and HP <= 250. and HP > 7.5:
a1, a2, a3 = 4.5212, 0.47242, 0.04820
elif self.kwargs["rpm"] == 1 and HP > 250.:
a1, a2, a3 = 7.4044, -0.06464, 0.05448
elif self.kwargs["rpm"] == 2 and HP <= 7.5:
a1, a2, a3 = 4.9298, 0.30118, 0.12630
elif self.kwargs["rpm"] == 2 and HP <= 250. and HP > 7.5:
a1, a2, a3 = 5.0999, 0.35861, 0.06052
elif self.kwargs["rpm"] == 2 and HP > 250.:
a1, a2, a3 = 4.6163, 0.88531, -0.02188
elif self.kwargs["motor"] == 1: # Totally enclosed, fan-cooled
if self.kwargs["rpm"] == 0 and HP <= 7.5:
a1, a2, a3 = 5.1058, 0.03316, 0.15374
elif self.kwargs["rpm"] == 0 and HP <= 250. and HP > 7.5:
a1, a2, a3 = 3.8544, 0.83311, 0.02399
elif self.kwargs["rpm"] == 0 and HP > 250.:
a1, a2, a3 = 5.3182, 1.08470, -0.05695
elif self.kwargs["rpm"] == 1 and HP <= 7.5:
a1, a2, a3 = 4.9687, -0.00930, 0.22616
elif self.kwargs["rpm"] == 1 and HP > 7.5:
a1, a2, a3 = 4.5347, 0.57065, 0.04609
elif self.kwargs["rpm"] == 2 and HP <= 7.5:
a1, a2, a3 = 5.1532, 0.28931, 0.14357
elif self.kwargs["rpm"] == 2 and HP > 7.5:
a1, a2, a3 = 5.3858, 0.31004, 0.07406
elif self.kwargs["motor"] == 2: # Explosion-proof
if self.kwargs["rpm"] == 0 and HP <= 7.5:
a1, a2, a3 = 5.3934, -0.00333, 0.15475
elif self.kwargs["rpm"] == 0 and HP > 7.5:
a1, a2, a3 = 4.4442, 0.60820, 0.05202
elif self.kwargs["rpm"] == 1 and HP <= 7.5:
a1, a2, a3 = 5.2851, 0.00048, 0.19949
elif self.kwargs["rpm"] == 1 and HP > 7.5:
a1, a2, a3 = 4.8178, 0.51086, 0.05293
elif self.kwargs["rpm"] == 2 and HP <= 7.5:
a1, a2, a3 = 5.4166, 0.31216, 0.10573
elif self.kwargs["rpm"] == 2 and HP > 7.5:
a1, a2, a3 = 5.5655, 0.31284, 0.07212
LnHP = log(HP)
C_motor = self.etapas*1.2 * exp(a1 + a2 * LnHP + a3 * LnHP**2)
if self.kwargs["compresor"] == 2:
C_trans = 0
elif self.kwargs["transmision"] == 0:
C_trans = 1.2*exp(3.689+0.8917*log(HP))
elif self.kwargs["transmision"] == 1:
C_trans = 1.2*exp(5.329+0.5048*log(HP))
elif self.kwargs["transmision"] == 2:
C_trans = 12000/(1.562+7.877/HP)
C_trans = self.etapas*C_trans*CI/BI
C_adq = C_comp+C_motor+C_trans
C_inst = C_adq * self.kwargs["f_install"]
self.C_comp = Currency(C_comp)
self.C_motor = Currency(C_motor)
self.C_trans = Currency(C_trans)
self.C_adq = Currency(C_adq)
self.C_inst = Currency(C_inst)
[docs]
def propTxt(self):
txt = "#---------------"
txt += translate("equipment", "Calculate properties")
txt += "-----------------#"+os.linesep
txt += self.propertiesToText(range(11))
if self.statusCoste:
txt += os.linesep+"#---------------"
txt += translate("equipment",
"Preliminary Cost Estimation")
txt += "-----------------#" + os.linesep
txt += self.propertiesToText(range(11, 23))
return txt
[docs]
@classmethod
def propertiesEquipment(cls):
l = [(translate("equipment", "Input Temperature"), "Tin", Temperature),
(translate("equipment", "Input Pressure"), "Pin", Pressure),
(translate("equipment", "Temperature increase"), "deltaT", DeltaT),
(translate("equipment", "Output Temperature"), "Tout", Temperature),
(translate("equipment", "Output Pressure"), "Pout", Pressure),
(translate("equipment", "Pressure increase"), "deltaP", DeltaP),
(translate("equipment", "Pressure ratio"), "razonCalculada",
Dimensionless),
(translate("equipment", "Thermodinamic mode"),
("TEXT_TERMODINAMICA", "termodinamica"), str),
(translate("equipment", "Power"), "power", Power),
(translate("equipment", "Efficiency"), "rendimientoCalculado",
Dimensionless),
(translate("equipment", "Especific capacities ratio"), "cp_cv",
Dimensionless),
(translate("equipment", "Base index"), "Base_index", float),
(translate("equipment", "Current index"), "Current_index", float),
(translate("equipment", "Install factor"), "f_install", float),
(translate("equipment", "Compressor Type"),
("TEXT_COMPRESOR", "compresor"), str),
(translate("equipment", "Transmission Type"),
("TEXT_TRANSMISION", "transmision"), str),
(translate("equipment", "Motor Type"), ("TEXT_MOTOR", "motor"), str),
(translate("equipment", "Motor RPM"), ("TEXT_RPM", "rpm"), str),
(translate("equipment", "Cost compressor"), "C_comp", Currency),
(translate("equipment", "Cost Transmission"), "C_trans", Currency),
(translate("equipment", "Cost motor"), "C_motor", Currency),
(translate("equipment", "Purchase Cost"), "C_adq", Currency),
(translate("equipment", "Installed Cost"), "C_inst", Currency)]
return l
[docs]
def writeStatetoJSON(self, state):
"""Write instance parameter to file"""
state["Pout"] = self.Pout
state["Tout"] = self.Tout
state["rendimientoCalculado"] = self.rendimientoCalculado
state["etapas"] = self.etapas
state["power"] = self.power
state["razonCalculada"] = self.razonCalculada
state["deltaT"] = self.deltaT
state["deltaP"] = self.deltaP
state["cp_cv"] = self.cp_cv
state["Pin"] = self.Pin
state["Tin"] = self.Tin
state["statusCoste"] = self.statusCoste
if self.statusCoste:
state["C_comp"] = self.C_comp
state["C_motor"] = self.C_motor
state["C_trans"] = self.C_trans
state["C_adq"] = self.C_adq
state["C_inst"] = self.C_inst
[docs]
def readStatefromJSON(self, state):
"""Load instance parameter from saved file"""
self.Pout = Pressure(state["Pout"])
self.Tout = Temperature(state["Tout"])
self.rendimientoCalculado = Dimensionless(state["rendimientoCalculado"])
self.etapas = state["etapas"]
self.power = Power(state["power"])
self.razonCalculada = Dimensionless(state["razonCalculada"])
self.deltaT = DeltaT(state["deltaT"])
self.deltaP = DeltaP(state["deltaP"])
self.cp_cv = Dimensionless(state["cp_cv"])
self.Pin = Pressure(state["Pin"])
self.Tin = Temperature(state["Tin"])
self.statusCoste = state["statusCoste"]
if self.statusCoste:
self.C_comp = Currency(state["C_comp"])
self.C_motor = Currency(state["C_motor"])
self.C_trans = Currency(state["C_trans"])
self.C_adq = Currency(state["C_adq"])
self.C_inst = Currency(state["C_inst"])
self.salida = [None]
[docs]
class Turbine(equipment):
"""Class to model a gas expander, turbine
Parameters:
entrada: Instance of class Corriente to define the imput stream
metodo: Value to specified the calculate type
0 - Output pressure and efficiency
1 - Compression ratio and efficiency
2 - Real work and efficiency
3 - Output pressure and real work
4 - Compression ratio and real work
5 - Calculate flow ot stream to compress known output pressure,
work and efficiency
termodinamica: Index with the thermodynamic model
0 - Adiabatic
1 - Polytropic
2 - Isothermic
Pout: Output pressure
razon: Pressures ratio
rendimeinto: Turbine efficiency
trabajo: Turbine work
Coste
f_install: instalation factor
base: Base index
actual: Current index
>>> from lib.corriente import Corriente
>>> corriente=Corriente(T=400, P=101325, caudalMasico=0.1, \
fraccionMasica=[1., 0, 0, 0])
>>> turbina=Turbine(entrada=corriente, metodo=1, razon=0.3, rendimiento=1.)
>>> print(turbina.power.MJh)
-69.138664576
>>> print(turbina.C_inst)
30713.7301133
"""
title = translate("equipment", "Turbine")
help = ""
kwargs = {
"entrada": None,
"metodo": 0,
"termodinamica": 0,
"Pout": 0.0,
"razon": 0.0,
"rendimiento": 0.0,
"trabajo": 0.0,
"f_install": 1.5,
"Base_index": None,
"Current_index": None}
kwargsInput = ("entrada", )
kwargsValue = ("Pout", "razon", "rendimiento", "trabajo")
kwargsList = ("metodo", "termodinamica")
calculateValue = ("power", "cp_cv", "razonCalculada",
"rendimientoCalculado")
calculateCostos = ("C_adq", "C_inst")
indiceCostos = 2
TEXT_METODO = [
translate("equipment", "Specify out pressure and efficiency"),
translate("equipment", "Specify pressure ratio and efficiency"),
translate("equipment", "Specify actual power and efficiency"),
translate("equipment", "Specify out pressure and actual power"),
translate("equipment", "Specify pressure ratio and actual power"),
translate("equipment", "Calculate input flowrate")]
TEXT_TERMODINAMICA = [translate("equipment", "Adiabatic"),
translate("equipment", "Polytropic"),
translate("equipment", "Isotermic")]
@property
def isCalculable(self):
if self.kwargs["f_install"] and self.kwargs["Base_index"] and \
self.kwargs["Current_index"]:
self.statusCoste = True
else:
self.statusCoste = False
if not self.kwargs["entrada"]:
if self.kwargs["metodo"] == 5:
if (self.kwargs["razon"] or self.kwargs["Pout"]) and \
self.kwargs["trabajo"] and self.kwargs["rendimiento"]:
self.status = 1
self.msg = ""
return True
else:
self.msg = translate("equipment", "undefined variables")
self.status = 0
else:
self.msg = translate("equipment", "undefined input")
self.status = 0
else:
if self.kwargs["metodo"] == 0:
valores = self.kwargs["Pout"] and self.kwargs["rendimiento"]
elif self.kwargs["metodo"] == 1:
valores = self.kwargs["razon"] and self.kwargs["rendimiento"]
elif self.kwargs["metodo"] == 2:
valores = self.kwargs["trabajo"] and self.kwargs["rendimiento"]
elif self.kwargs["metodo"] == 3:
valores = self.kwargs["Pout"] and self.kwargs["trabajo"]
elif self.kwargs["metodo"] == 4:
valores = self.kwargs["razon"] and self.kwargs["trabajo"]
elif self.kwargs["metodo"] == 5:
valores = self.kwargs["razon"] and self.kwargs["rendimiento"]
if valores:
self.status = 1
self.msg = ""
return True
else:
self.msg = translate("equipment", "undefined variables")
self.status = 0
[docs]
def calculo(self):
self.entrada = self.kwargs["entrada"]
self.Pout = Pressure(self.kwargs["Pout"])
self.razon = Dimensionless(self.kwargs["razon"])
self.rendimientoCalculado = Dimensionless(self.kwargs["rendimiento"])
self.power = Power(-abs(self.kwargs["trabajo"]))
def f(Pout, rendimiento):
W_ideal = self.__Wideal(Pout)
power = W_ideal*self.entrada.caudalmasico.gs*rendimiento
return power
self.cp_cv = self.entrada.Gas.cp_cv
if self.kwargs["metodo"] in [0, 3] or \
(self.kwargs["metodo"] == 5 and self.Pout):
self.razon = Dimensionless(self.Pout/self.entrada.P)
elif self.kwargs["metodo"] in [1, 4] or \
(self.kwargs["metodo"] == 5 and self.razon):
self.Pout = Pressure(self.entrada.P*self.razon)
if self.kwargs["metodo"] in [0, 1]:
Wideal = self.__Wideal(self.Pout)
G = self.entrada.caudalmasico.gs
self.power = Power(Wideal*G*self.rendimientoCalculado)
elif self.kwargs["metodo"] == 2:
def function(P):
return f(P, self.rendimientoCalculado)-self.power
self.Pout = Pressure(fsolve(function, self.entrada.P.atm+1), "atm")
self.razon = Dimensionless(self.Pout/self.entrada.P)
elif self.kwargs["metodo"] in [3, 4]:
def function(rendimiento):
return f(self.Pout.atm, rendimiento)-self.power
self.rendimientoCalculado = Dimensionless(fsolve(function, 0.5))
elif self.kwargs["metodo"] == 5:
Wideal = self.__Wideal(self.Pout)
G = MassFlow(self.power/self.rendimientoCalculado/Wideal, "gs")
self.Tout = self.__Tout(Wideal)
self.entrada = self.entrada.clone(caudalMasico=G)
Wideal = self.__Wideal(self.Pout)
self.Tout = Temperature(self.__Tout(Wideal))
self.razonCalculada = Dimensionless(self.Pout/self.entrada.P)
self.salida = [self.entrada.clone(T=self.Tout, P=self.Pout)]
self.deltaT = DeltaT(self.salida[0].T-self.entrada.T)
self.deltaP = DeltaP(self.salida[0].P-self.entrada.P)
self.Pin = self.entrada.P
self.Tin = self.entrada.T
def __Wideal(self, Pout):
"""Calculo del trabajo ideal"""
if self.kwargs["termodinamica"] == 0:
cpv = self.cp_cv
else:
cpv = 1.40388
if self.kwargs["termodinamica"] == 2:
W_ideal = R*self.entrada.T/self.entrada.M*log(Pout/self.entrada.P)
else:
W_ideal = R*self.entrada.T*cpv/self.entrada.M/(cpv-1) * \
((Pout/self.entrada.P)**((cpv-1)/cpv)-1)
return W_ideal
def __Tout(self, W_ideal):
"""Cálculo de la temperatura de salida"""
if self.kwargs["termodinamica"] == 0:
cpv = self.cp_cv
else:
cpv = 1.40388
if self.kwargs["termodinamica"] == 2:
Tout = self.entrada.T
else:
Tout = (self.Pout/self.entrada.P)**((cpv-1)/cpv)*self.entrada.T + \
(1-self.rendimientoCalculado)/self.rendimientoCalculado * \
(W_ideal/self.entrada.Gas.cv)
return Tout
[docs]
def coste(self):
HP = abs(self.power.hp)
CI = self.kwargs["Current_index"]
BI = self.kwargs["Base_index"]
if self.salida[0].P.psi >= 14.696:
C = 0.31*HP**0.81*1000
else:
C = 0.69*HP**0.81*1000
self.C_adq = Currency(C * CI / BI)
self.C_inst = Currency(self.C_adq * self.kwargs["f_install"])
[docs]
def propTxt(self):
txt = "#---------------"
txt += translate("equipment", "Calculate properties")
txt += "-----------------#"+os.linesep
txt += self.propertiesToText(range(11))
if self.statusCoste:
txt += os.linesep+"#---------------"
txt += translate("equipment", "Preliminary Cost Estimation")
txt += "-----------------#" + os.linesep
txt += self.propertiesToText(range(11, 16))
return txt
[docs]
@classmethod
def propertiesEquipment(cls):
l = [(translate("equipment", "Input Temperature"), "Tin", Temperature),
(translate("equipment", "Output Temperature"), "Tout", Temperature),
(translate("equipment", "Temperature increase"), "deltaT", DeltaT),
(translate("equipment", "Input Pressure"), "Pin", Pressure),
(translate("equipment", "Output Pressure"), "Pout", Pressure),
(translate("equipment", "Pressure increase"), "deltaP", DeltaP),
(translate("equipment", "Pressure ratio"), "razonCalculada",
Dimensionless),
(translate("equipment", "Thermodinamic mode"),
("TEXT_TERMODINAMICA", "termodinamica"), str),
(translate("equipment", "Power"), "power", Power),
(translate("equipment", "Efficiency"), "rendimientoCalculado",
Dimensionless),
("Cp/Cv", "cp_cv", Dimensionless),
(translate("equipment", "Base index"), "Base_index", float),
(translate("equipment", "Current index"), "Current_index", float),
(translate("equipment", "Install factor"), "f_install", float),
(translate("equipment", "Purchase Cost"), "C_adq", Currency),
(translate("equipment", "Installed Cost"), "C_inst", Currency)]
return l
[docs]
def writeStatetoJSON(self, state):
"""Write instance parameter to file"""
state["Pout"] = self.Pout
state["Tout"] = self.Tout
state["rendimientoCalculado"] = self.rendimientoCalculado
state["power"] = self.power
state["razonCalculada"] = self.razonCalculada
state["razon"] = self.razon
state["deltaT"] = self.deltaT
state["deltaP"] = self.deltaP
state["cp_cv"] = self.cp_cv
state["Pin"] = self.Pin
state["Tin"] = self.Tin
state["cp_cv"] = self.cp_cv
state["statusCoste"] = self.statusCoste
if self.statusCoste:
state["C_adq"] = self.C_adq
state["C_inst"] = self.C_inst
[docs]
def readStatefromJSON(self, state):
"""Load instance parameter from saved file"""
self.Pout = Pressure(state["Pout"])
self.Tout = Temperature(state["Tout"])
self.rendimientoCalculado = Dimensionless(state["rendimientoCalculado"])
self.power = Power(state["power"])
self.razonCalculada = Dimensionless(state["razonCalculada"])
self.razon = Dimensionless(state["razon"])
self.deltaT = DeltaT(state["deltaT"])
self.deltaP = DeltaP(state["deltaP"])
self.cp_cv = Dimensionless(state["cp_cv"])
self.Pin = Pressure(state["Pin"])
self.Tin = Temperature(state["Tin"])
self.cp_cv = Dimensionless(state["cp_cv"])
self.statusCoste = state["statusCoste"]
if self.statusCoste:
self.C_adq = Currency(state["C_adq"])
self.C_inst = Currency(state["C_inst"])
self.salida = [None]
if __name__ == '__main__':
import doctest
doctest.testmod()
