Need help with the values

[fuzz08]

Hello. I have recently started to start learning TESPy. I have formulated a model on my own. It generates an output.
However, I am unsure about the results. It would be helpful if I could get some help (with the values).

I tried running the same code in Visual Studio, but it's not generating the output. I am getting so many errors. But with Spyder, it's generating an output.

The code -

# %%[sec_1]
from tespy.components.basics.cycle_closer import CycleCloser
from tespy.networks import Network
from tespy.components import (
    CycleCloser, Pipe, Pump, Valve, SimpleHeatExchanger,Splitter,Merge
)
from tespy.connections import Connection

fluid_list = ['INCOMP::Water']
nw = Network(fluids=fluid_list)
nw.set_attr(T_unit='C', p_unit='bar', h_unit='kJ / kg')

# central heating plant
hs = SimpleHeatExchanger('heat source')
cc = CycleCloser('cycle closer')
pu_1 = Pump('feed pump_1')
pu_2 = Pump('feed pump_2')

# consumer_1
cons_1 = SimpleHeatExchanger('consumer_1')
val_1 = Valve('control valve_1')

# consumer_2
cons_2 = SimpleHeatExchanger('consumer_2')
val_2 = Valve('control valve_2')


#splitterand Merge
sp = Splitter("splitter")
mg = Merge("merge")

# pipes
pipe_feed_1 = Pipe('feed pipe_1')
pipe_return_1 = Pipe('return pipe_1')

pipe_feed_2 = Pipe('feed pipe_2')
pipe_return_2 = Pipe('return pipe_2')


# connections
c0 = Connection(cc, "out1", hs, "in1", label="0")
c1 = Connection(hs, "out1", pu_1, "in1", label="1")
c2 = Connection(pu_1, "out1", pipe_feed_1, "in1", label="2")

c3 = Connection(pipe_feed_1, "out1", sp, "in1", label="3")
c11 = Connection(sp, "out1", cons_1, "in1", label="11")
c13 = Connection(sp, "out2",  pu_2, "in1", label="13")


c10 = Connection(pu_2,"out1", pipe_feed_2,"in1",label="10")
c9 =  Connection(pipe_feed_2, "out1", cons_2, "in1", label="9")
c8 = Connection(cons_2, "out1", val_2, "in1", label="8")
c7 = Connection(val_2, "out1", pipe_return_2, "in1", label="7")

c14 = Connection(pipe_return_2, "out1", mg, "in1", label="14")
c12 = Connection(cons_1, "out1", mg, "in2", label="12")

c4 = Connection(mg, "out1", val_1, "in1", label="4")
c5 = Connection(val_1, "out1", pipe_return_1, "in1", label="5")
c6 = Connection(pipe_return_1, "out1", cc, "in1", label="6")

nw.add_conns(c0, c1, c2, c3, c4, c5, c6,c7,c8,c9,c10,c11,c12,c13,c14)


# %%[sec_2]
cons_1.set_attr(Q=-10000, pr=0.98)
cons_2.set_attr(Q=-9000, pr=0.98)
hs.set_attr(pr=1)
pu_1.set_attr(eta_s=0.75)
pipe_feed_1.set_attr(Q=-250, pr=0.98)
pipe_feed_2.set_attr(Q=-200, pr=0.98)
pipe_return_1.set_attr(Q=-180, pr=0.98)
pipe_return_2.set_attr(Q=-160, pr=0.98)

c1.set_attr(T=90, p=10, fluid={'Water': 1})
c2.set_attr(p=13)
c4.set_attr(T=65)
c8.set_attr(T=60)
c9.set_attr(T=40)
c10.set_attr(T=45)


nw.solve(mode="design")
nw.print_results()

# %%[sec_3]
pipe_feed_1.set_attr(
    ks=0.0005,  # pipe's roughness in meters
    L=100,  # length in m
    D="var",  # diameter in m
)
pipe_return_1.set_attr(
    ks=0.0004,  # pipe's roughness in meters
    L=90,  # length in m
    D="var",  # diameter in m
)
pipe_feed_2.set_attr(
    ks=0.0005,  # pipe's roughness in meters
    L=80,  # length in m
    D="var",  # diameter in m
)
pipe_return_2.set_attr(
    ks=0.0004,  # pipe's roughness in meters
    L=70,  # length in m
    D="var",  # diameter in m
)

nw.solve(mode="design")
nw.print_results()


# %%[sec_4]
pipe_feed_1.set_attr(D=pipe_feed_1.D.val, pr=None)
pipe_return_1.set_attr(D=pipe_return_1.D.val, pr=None)
pipe_feed_2.set_attr(D=pipe_feed_2.D.val, pr=None)
pipe_return_2.set_attr(D=pipe_return_2.D.val, pr=None)

pipe_feed_1.set_attr(
    Tamb=0,  # ambient temperature level in network's temperature unit
    kA="var"  # area independent heat transfer coefficient
)
pipe_return_1.set_attr(
    Tamb=0,  # ambient temperature level in network's temperature unit
    kA="var"  # area independent heat transfer coefficient
)

pipe_feed_2.set_attr(
    Tamb=0,  # ambient temperature level in network's temperature unit
    kA="var"  # area independent heat transfer coefficient
)
pipe_return_2.set_attr(
    Tamb=0,  # ambient temperature level in network's temperature unit
    kA="var"  # area independent heat transfer coefficient
)
nw.solve(mode="design")
nw.print_results()

[fwitte]

@fuzz08,

thanks for reaching out. I am not quite sure, what the differences between spyder and vscode are here. In any case, the parameter specifcations as you provided do not seem reasonable to me, i.e. temperature 9 below temperature 8, temperature 4 above the temperatures of the inflows. Together with the other specifications these seem fishy not only from context, but e.g. pressure cannot be determined correcty. I suggest, you build up your system component by component and always try to find a solution before you add new components. This way, it becomes much easier to understand, what specifications are required for what reason.

I can send you a working "solution", but I if you want to learn the software, you should try the other approach first.

Best

Francesco

[fuzz08]

Hey @fwitte,
Thanks for the reply !

I made some adjustments to the District Model Heating example by adding more components of mine. If you could provide me a "working solution" that would help a lot. It would help me understand the system a bit better. Thanks !

[fwitte]

#[sec_1]

from tespy.components.basics.cycle_closer import CycleCloser
from tespy.networks import Network
from tespy.components import (
CycleCloser, Pipe, Pump, Valve, SimpleHeatExchanger,Splitter,Merge
)
from tespy.connections import Connection

nw = Network()
nw.set_attr(T_unit='C', p_unit='bar', h_unit='kJ / kg')
# central heating plant

hs = SimpleHeatExchanger('heat source')
cc = CycleCloser('cycle closer')
pu_1 = Pump('feed pump_1')
pu_2 = Pump('feed pump_2')
# consumer_1

cons_1 = SimpleHeatExchanger('consumer_1')
val_1 = Valve('control valve_1')
# consumer_2

cons_2 = SimpleHeatExchanger('consumer_2')
val_2 = Valve('control valve_2')

#splitterand Merge
sp = Splitter("splitter")
mg = Merge("merge")
# pipes

pipe_feed_1 = Pipe('feed pipe_1')
pipe_return_1 = Pipe('return pipe_1')

pipe_feed_2 = Pipe('feed pipe_2')
pipe_return_2 = Pipe('return pipe_2')
# connections

c0 = Connection(cc, "out1", hs, "in1", label="0")
c1 = Connection(hs, "out1", pu_1, "in1", label="1")
c2 = Connection(pu_1, "out1", pipe_feed_1, "in1", label="2")

c3 = Connection(pipe_feed_1, "out1", sp, "in1", label="3")
c11 = Connection(sp, "out1", cons_1, "in1", label="11")
c13 = Connection(sp, "out2", pu_2, "in1", label="13")

c10 = Connection(pu_2,"out1", pipe_feed_2,"in1",label="10")
c9 = Connection(pipe_feed_2, "out1", cons_2, "in1", label="9")
c8 = Connection(cons_2, "out1", val_2, "in1", label="8")
c7 = Connection(val_2, "out1", pipe_return_2, "in1", label="7")

c14 = Connection(pipe_return_2, "out1", mg, "in1", label="14")
c12 = Connection(cons_1, "out1", mg, "in2", label="12")

c4 = Connection(mg, "out1", val_1, "in1", label="4")
c5 = Connection(val_1, "out1", pipe_return_1, "in1", label="5")
c6 = Connection(pipe_return_1, "out1", cc, "in1", label="6")

nw.add_conns(c0, c1, c2, c3, c4, c5, c6,c7,c8,c9,c10,c11,c12,c13,c14)
#[sec_2]

cons_1.set_attr(Q=-10000, pr=0.98)
cons_2.set_attr(Q=-9000, pr=0.98)
hs.set_attr(pr=1)
pu_1.set_attr(eta_s=0.75)
pu_2.set_attr(eta_s=0.75)
pipe_feed_1.set_attr(Q=-250, pr=0.98)
pipe_feed_2.set_attr(Q=-200, pr=0.98)
pipe_return_1.set_attr(Q=-180, pr=0.98)
pipe_return_2.set_attr(Q=-160, pr=0.98)
val_2.set_attr(pr=1)

for c in nw.conns["object"]:
    c.set_attr(p0=10)

c1.set_attr(T=90, p=10, fluid={"INCOMP::Water": 1})
c2.set_attr(p=13)
c4.set_attr()
c8.set_attr(T=60)
c12.set_attr(T=60)
c9.set_attr()
c10.set_attr()

nw.solve(mode="design")
nw.print_results()
#[sec_3]

pipe_feed_1.set_attr(
ks=0.0005, # pipe's roughness in meters
L=100, # length in m
D="var", # diameter in m
)
pipe_return_1.set_attr(
ks=0.0004, # pipe's roughness in meters
L=90, # length in m
D="var", # diameter in m
)
pipe_feed_2.set_attr(
ks=0.0005, # pipe's roughness in meters
L=80, # length in m
D="var", # diameter in m
)
pipe_return_2.set_attr(
ks=0.0004, # pipe's roughness in meters
L=70, # length in m
D="var", # diameter in m
)

nw.solve(mode="design")
nw.print_results()
#[sec_4]

pipe_feed_1.set_attr(D=pipe_feed_1.D.val, pr=None)
pipe_return_1.set_attr(D=pipe_return_1.D.val, pr=None)
pipe_feed_2.set_attr(D=pipe_feed_2.D.val, pr=None)
pipe_return_2.set_attr(D=pipe_return_2.D.val, pr=None)

pipe_feed_1.set_attr(
Tamb=0, # ambient temperature level in network's temperature unit
kA="var" # area independent heat transfer coefficient
)
pipe_return_1.set_attr(
Tamb=0, # ambient temperature level in network's temperature unit
kA="var" # area independent heat transfer coefficient
)

pipe_feed_2.set_attr(
Tamb=0, # ambient temperature level in network's temperature unit
kA="var" # area independent heat transfer coefficient
)
pipe_return_2.set_attr(
Tamb=0, # ambient temperature level in network's temperature unit
kA="var" # area independent heat transfer coefficient
)
nw.solve(mode="design")
nw.print_results()