fluid_list = ['O2', 'H2O', 'H2']
nw = network(fluids=fluid_list,
T_unit='C',
p_unit='bar',
v_unit='l / s',
iterinfo=False)
# %% components
fw = source('feed water')
oxy = sink('oxygen sink')
hydro = sink('hydrogen sink')
cw_cold = source('cooling water source')
cw_hot = sink('cooling water sink')
comp = compressor('compressor', eta_s=0.9)
el = water_electrolyzer('electrolyzer')
# %% connections
fw_el = connection(fw, "out1", el, "in2")
el_comp = connection(el, 'out3', comp, 'in1')
comp_hydro = connection(comp, 'out1', hydro, 'in1')
el_oxy = connection(el, 'out2', oxy, 'in1')
cw_cold_el = connection(cw_cold, 'out1', el, 'in1')
el_cw_hot = connection(el, 'out1', cw_hot, 'in1')
nw.add_conns(fw_el, el_comp, comp_hydro, el_oxy, cw_cold_el, el_cw_hot)
# %% busses
h_unit='kJ / kg',
m_unit='kg / s')
# %% components
# sources & sinks
cc = cycle_closer('coolant cycle closer')
cb = source('consumer back flow')
cf = sink('consumer feed flow')
amb = source('ambient air')
amb_out1 = sink('sink ambient 1')
amb_out2 = sink('sink ambient 2')
# ambient air system
sp = splitter('splitter')
fan = compressor('fan')
# consumer system
cd = condenser('condenser')
dhp = pump('district heating pump')
cons = heat_exchanger_simple('consumer')
# evaporator system
ves = valve('valve')
dr = drum('drum')
ev = heat_exchanger('evaporator')
su = heat_exchanger('superheater')
erp = pump('evaporator reciculation pump')
cd = condenser('condenser')
dhp = pump('district heating pump')
cons = heat_exchanger_simple('consumer')
# evaporator system
ves = valve('valve')
dr = drum('drum')
ev = heat_exchanger('evaporator')
su = heat_exchanger('superheater')
erp = pump('evaporator reciculation pump')
# compressor-system
cp1 = compressor('compressor 1')
cp2 = compressor('compressor 2')
ic = heat_exchanger('intercooler')
# %% connections
# consumer system
c_in_cd = connection(cc, 'out1', cd, 'in1')
cb_dhp = connection(cb, 'out1', dhp, 'in1')
dhp_cd = connection(dhp, 'out1', cd, 'in2')
cd_cons = connection(cd, 'out2', cons, 'in1')
cons_cf = connection(cons, 'out1', cf, 'in1')
nw.add_conns(c_in_cd, cb_dhp, dhp_cd, cd_cons, cons_cf)
# consumer system
cd = condenser('condenser')
dhp = pump('district heating pump')
cons = heat_exchanger_simple('consumer')
# evaporator system
va = valve('valve')
dr = drum('drum')
ev = heat_exchanger('evaporator')
erp = pump('evaporator reciculation pump')
# compressor-system
cp = compressor('compressor')
# %% connections
# consumer system
c_in_cd = connection(cc, 'out1', cd, 'in1')
cb_dhp = connection(cb, 'out1', dhp, 'in1')
dhp_cd = connection(dhp, 'out1', cd, 'in2')
cd_cons = connection(cd, 'out2', cons, 'in1')
cons_cf = connection(cons, 'out1', cf, 'in1')
nw.add_conns(c_in_cd, cb_dhp, dhp_cd, cd_cons, cons_cf)
# connection condenser - evaporator system
nw = network(fluids=fluid_list,
p_unit='bar',
T_unit='C',
h_unit='kJ / kg',
v_unit='l / s',
iterinfo=False)
t_dh_in = 50
t_dh_out = 124 # might change due to 4GDH
# %% components
# gas turbine part
air_source = source('air source')
air_compressor = compressor('air compressor')
combustion_chamber = combustion_chamber('combustion chamber')
fuel_source = source('fuel source')
gas_turbine = turbine('gas turbine')
steam_generator = heat_exchanger('steamgenerator')
exhaust_gas = sink('exhaust gas')
# steam turbine part
fs_source = source('steam source')
steam_turbine = turbine('steam turbine')
dh_heat_exchanger = heat_exchanger('district heating heat exchanger')
pump = pump('pump')
fs_sink = sink('steam sink')
class HeatPump(object):
# define the structure of heat pump
# %% network
nw = network(fluids=['water', 'NH3', 'air'],
T_unit='C',
p_unit='bar',
h_unit='kJ / kg',
m_unit='kg / s')
# %% components
# sources & sinks
cc = cycle_closer('coolant cycle closer')
cb = source('consumer back flow')
cf = sink('consumer feed flow')
amb = source('ambient air')
amb_out1 = sink('sink ambient 1')
amb_out2 = sink('sink ambient 2')
# ambient air system
sp = splitter('splitter')
pu = pump('pump')
# consumer system
cd = condenser('condenser')
dhp = pump('district heating pump')
cons = heat_exchanger_simple('consumer')
# evaporator system
ves = valve('valve')
dr = drum('drum')
ev = heat_exchanger('evaporator')
su = heat_exchanger('superheater')
erp = pump('evaporator reciculation pump')
# compressor-system
cp1 = compressor('compressor 1')
cp2 = compressor('compressor 2')
ic = heat_exchanger('intercooler')
# %% connections
# consumer system
c_in_cd = connection(cc, 'out1', cd, 'in1')
cb_dhp = connection(cb, 'out1', dhp, 'in1')
dhp_cd = connection(dhp, 'out1', cd, 'in2')
cd_cons = connection(cd, 'out2', cons, 'in1')
cons_cf = connection(cons, 'out1', cf, 'in1')
nw.add_conns(c_in_cd, cb_dhp, dhp_cd, cd_cons, cons_cf)
# connection condenser - evaporator system
cd_ves = connection(cd, 'out1', ves, 'in1')
nw.add_conns(cd_ves)
# evaporator system
ves_dr = connection(ves, 'out1', dr, 'in1')
dr_erp = connection(dr, 'out1', erp, 'in1')
erp_ev = connection(erp, 'out1', ev, 'in2')
ev_dr = connection(ev, 'out2', dr, 'in2')
dr_su = connection(dr, 'out2', su, 'in2')
nw.add_conns(ves_dr, dr_erp, erp_ev, ev_dr, dr_su)
amb_p = connection(amb, 'out1', pu, 'in1')
p_sp = connection(pu, 'out1', sp, 'in1')
sp_su = connection(sp, 'out1', su, 'in1')
su_ev = connection(su, 'out1', ev, 'in1')
ev_amb_out = connection(ev, 'out1', amb_out1, 'in1')
nw.add_conns(amb_p, p_sp, sp_su, su_ev, ev_amb_out)
# connection evaporator system - compressor system
su_cp1 = connection(su, 'out2', cp1, 'in1')
nw.add_conns(su_cp1)
# compressor-system
cp1_he = connection(cp1, 'out1', ic, 'in1')
he_cp2 = connection(ic, 'out1', cp2, 'in1')
cp2_c_out = connection(cp2, 'out1', cc, 'in1')
sp_ic = connection(sp, 'out2', ic, 'in2')
ic_out = connection(ic, 'out2', amb_out2, 'in1')
nw.add_conns(cp1_he, he_cp2, sp_ic, ic_out, cp2_c_out)
def __init__(self, q, eff, Temp):
r"""
:param Temp:
:param q: q output
:param eff: efficient of each part in pump
"""
self.q = q
self.eff = eff
self.Temp = Temp
def caculation(self):
self.set_attr()
HeatPump.nw.solve('design')
P = [
HeatPump.cp1.P.val, HeatPump.cp2.P.val, HeatPump.erp.P.val,
HeatPump.pu.P.val
]
P_total = sum(map(abs, P))
P = list(map(abs, P))
COP = self.q / P_total
# T = [HeatPump.su_cp1.T.val, HeatPump.cp2_c_out.T.val, HeatPump.cd_ves.T.val, HeatPump.su_ev.T.val]
# p = [HeatPump.su_cp1.p.val, HeatPump.cp2_c_out.p.val, HeatPump.cd_ves.p.val, HeatPump.su_ev.p.val,
# HeatPump.cp1_he.p.val]
return P, P_total, COP
def set_attr(self):
r"""
# %% set the attribution of the heat pump
:return: heat output of the heat pump
"""
HeatPump.cd.set_attr(pr1=0.99,
pr2=0.99,
ttd_u=15,
design=['pr2', 'ttd_u'],
offdesign=['zeta2', 'kA'])
HeatPump.dhp.set_attr(eta_s=self.eff,
design=['eta_s'],
offdesign=['eta_s_char'])
HeatPump.cons.set_attr(pr=0.99, design=['pr'], offdesign=['zeta'])
# water pump
HeatPump.pu.set_attr(eta_s=self.eff,
design=['eta_s'],
offdesign=['eta_s_char'])
# evaporator system
kA_char1 = ldc('heat exchanger', 'kA_char1', 'DEFAULT', char_line)
kA_char2 = ldc('heat exchanger', 'kA_char2', 'EVAPORATING FLUID',
char_line)
HeatPump.ev.set_attr(pr1=0.98,
pr2=0.99,
ttd_l=5,
kA_char1=kA_char1,
kA_char2=kA_char2,
design=['pr1', 'ttd_l'],
offdesign=['zeta1', 'kA'])
HeatPump.su.set_attr(pr1=0.98,
pr2=0.99,
ttd_u=2,
design=['pr1', 'pr2', 'ttd_u'],
offdesign=['zeta1', 'zeta2', 'kA'])
HeatPump.erp.set_attr(eta_s=self.eff,
design=['eta_s'],
offdesign=['eta_s_char'])
# compressor system
HeatPump.cp1.set_attr(eta_s=self.eff,
design=['eta_s'],
offdesign=['eta_s_char'])
HeatPump.cp2.set_attr(eta_s=self.eff,
pr=3,
design=['eta_s'],
offdesign=['eta_s_char'])
HeatPump.ic.set_attr(pr1=0.99,
pr2=0.98,
design=['pr1', 'pr2'],
offdesign=['zeta1', 'zeta2', 'kA'])
# %% connection parametrization
# condenser system
HeatPump.c_in_cd.set_attr(fluid={'air': 0, 'NH3': 1, 'water': 0})
HeatPump.cb_dhp.set_attr(T=20,
p=10,
fluid={
'air': 0,
'NH3': 0,
'water': 1
})
HeatPump.cd_cons.set_attr(T=self.Temp)
HeatPump.cons_cf.set_attr(h=ref(HeatPump.cb_dhp, 1, 0),
p=ref(HeatPump.cb_dhp, 1, 0))
# evaporator system cold side
HeatPump.erp_ev.set_attr(m=ref(HeatPump.ves_dr, 1.25, 0), p0=5)
HeatPump.su_cp1.set_attr(p0=5, h0=1700)
# evaporator system hot side
# pumping at constant rate in partload
HeatPump.amb_p.set_attr(T=12,
p=2,
fluid={
'air': 0,
'NH3': 0,
'water': 1
},
offdesign=['v'])
HeatPump.sp_su.set_attr(offdesign=['v'])
HeatPump.ev_amb_out.set_attr(p=2, T=9, design=['T'])
# compressor-system
HeatPump.he_cp2.set_attr(Td_bp=5, p0=20, design=['Td_bp'])
HeatPump.ic_out.set_attr(T=15, design=['T'])
# %% key paramter
HeatPump.cons.set_attr(Q=self.q)
fluid_list = ['Ar', 'N2', 'O2', 'CO2', 'CH4', 'H2O']
nw = network(fluids=fluid_list,
p_unit='bar',
T_unit='C',
h_unit='kJ / kg',
p_range=[1, 15],
T_range=[10, 1200],
h_range=[500, 4000])
t_dh_in = 50
t_dh_out = 124 # might change due to 4GDH
# %% components
# gas turbine part
compressor_gtp = compressor('compressor')
combustion_chamber = combustion_chamber('combustion')
gas_turbine = turbine('gas turbine')
fuel_source = source('fuel source')
combustion_air_source = source('ambient air')
rauchgas = sink('Rauchgas')
steam_generator_gas = heat_exchanger('Abhitzekessel')
# steam turbine part
fs_source = source('Frischdampf')
hp_turbine = turbine('Hochdruck Turbine')
splitter_1 = splitter('splitter hp-extraction')
lp_turbine = turbine('low pressure turbine')