I want to build a secondary loop system for a nuclear power plant. But I don't know how to set the initial parameters.

[WangJ501]

This system includes a steam generator, a high-pressure steam turbine, a steam water separator reheater, a low-pressure steam turbine, a condenser, six heat recovery heaters , and a deaerator.This will look a bit complicated.

This is the code.
second.py

from tespy.components import(SimpleHeatExchanger, SteamTurbine, Condenser, Pump, CycleCloser, Splitter, Separator, Merge,HeatExchanger,Sink,Source,Drum,Node,DropletSeparator)

cc=CycleCloser('CycleCloser')
sg=SimpleHeatExchanger('Steam Generator')
hptu1=SteamTurbine('HPTU1')
hptu2=SteamTurbine('HPTU2')
hptu3=SteamTurbine('HPTU3')
hptu4=SteamTurbine('HPTU4')
hptu5=SteamTurbine('HPTU5')
hptu6=SteamTurbine('HPTU6')
hp1=Splitter('hp1',num_out=2)
hp21=Splitter('hp21',num_out=2)
hp22=Splitter('hp22',num_out=2)
hp3=Splitter('hp3',num_out=2)
hp4=Splitter('hp4',num_out=2)
m1=Merge('m1',num_in=2)
split0=Splitter('split0',num_out=2)
from tespy.connections import Connection
#steam generator+high-pressure turbine
c1=Connection(sg,'out1',cc,'in1',label='1')
c2=Connection(hp1,'out1',hptu1,'in1',label='2')
c3=Connection(hp1,'out2',hptu2,'in1',label='3')
c4=Connection(hptu1,'out1',hp3,'in1',label='4')
c5=Connection(hptu2,'out1',hp4,'in1',label='5')
c6=Connection(hp3,'out1',hptu3,'in1',label='6')
c7=Connection(hp4,'out1',hptu4,'in1',label='7')
c8=Connection(hptu3,'out1',hp21,'in1',label='8')
c9=Connection(hptu4,'out1',hp22,'in1',label='9')
c10=Connection(hp21,'out1',hptu5,'in1',label='10')
c11=Connection(hp22,'out1',hptu6,'in1',label='11')
c12=Connection(hptu5,'out1',m1,'in1',label='12')
c13=Connection(hptu6,'out1',m1,'in2',label='13')
c91=Connection(split0,'out1',hp1,'in1',label='91')

#steam water separator reheater
rh1=HeatExchanger('RH1')
rh2=HeatExchanger('RH2')
sp=DropletSeparator('sp')
hp41=Splitter('hp41',num_out=2)
c14=Connection(m1,'out1',sp,'in1',label='14')
c15=Connection(sp,'out2',rh1,'in1',label='15')
c16=Connection(rh1,'out1',rh2,'in1',label='16')
c17=Connection(hp4,'out2',hp41,'in1',label='17')
c18=Connection(hp41,'out1',rh1,'in2',label='18')
c19=Connection(split0,'out2',rh2,'in2',label='19')

#low-pressure turbine
m2=Merge('m2',num_in=2)
lp1=Splitter('lp1',num_out=2)
lp2=Splitter('lp2',num_out=2)
lp3=Splitter('lp3',num_out=2)
lp4=Splitter('lp4',num_out=2)
lp5=Splitter('lp5',num_out=2)
lp6=Splitter('lp6',num_out=2)
lp7=Splitter('lp7',num_out=2)
lp8=Splitter('lp8',num_out=2)
lp9=Splitter('lp9',num_out=2)
lptu1=SteamTurbine('LPTU1')
lptu2=SteamTurbine('LPTU2')
lptu3=SteamTurbine('LPTU3')
lptu4=SteamTurbine('LPTU4')
lptu5=SteamTurbine('LPTU5')
lptu6=SteamTurbine('LPTU6')
lptu7=SteamTurbine('LPTU7')
lptu8=SteamTurbine('LPTU8')
lptu9=SteamTurbine('LPTU9')
lptu10=SteamTurbine('LPTU10')
c20=Connection(rh2,'out1',lp1,'in1',label='20')
c21=Connection(lp1,'out1',lptu1,'in1',label='21')
c22=Connection(lp1,'out2',lptu2,'in1',label='22')
c23=Connection(lptu1,'out1',lp2,'in1',label='23')
c24=Connection(lptu2,'out1',lp3,'in1',label='24')
c25=Connection(lp2,'out1',lptu3,'in1',label='25')
c26=Connection(lp3,'out1',lptu4,'in1',label='26')
c27=Connection(lptu3,'out1',lp4,'in1',label='27')
c28=Connection(lptu4,'out1',lp5,'in1',label='28')
c29=Connection(lp4,'out1',lptu5,'in1',label='29')
c30=Connection(lp5,'out1',lptu6,'in1',label='30')
c31=Connection(lptu5,'out1',lp6,'in1',label='31')
c32=Connection(lptu6,'out1',lp7,'in1',label='32')
c33=Connection(lp6,'out1',lptu7,'in1',label='33')
c34=Connection(lp7,'out1',lptu8,'in1',label='34')
c35=Connection(lptu7,'out1',lp8,'in1',label='35')
c36=Connection(lp8,'out1',lptu9,'in1',label='36')
c37=Connection(lp9,'out1',lptu10,'in1',label='37')
c38=Connection(lptu9,'out1',m2,'in1',label='38')
c39=Connection(lptu10,'out1',m2,'in2',label='39')
c93=Connection(lptu8,'out1',lp9,'in1',label='93')

#Condenser
cd=Condenser('Condenser')
hotwell=Node('Hot Well',num_in=2,num_out=1)#热井
so=Source('Source')
si=Sink('Sink')
c40=Connection(m2,'out1',cd,'in1',label='40')
c41=Connection(cd,'out1',hotwell,'in1',label='41')
c42=Connection(so,'out1',cd,'in2',label='42')
c43=Connection(cd,'out2',si,'in1',label='43')

#Condensate pump
cp=Pump('Cooling Pump')
c44=Connection(hotwell,'out1',cp,'in1',label='44')

#CycleCloser
c45=Connection(cc,'out1',split0,'in1',label='45')
#Shaft seal cooler
gc=SimpleHeatExchanger('GC')
c46=Connection(cp,'out1',gc,'in1',label='46')

#low-pressure heat recovery heaters
htr1=HeatExchanger('HTR1')
htr2=HeatExchanger('HTR2')
htr3=HeatExchanger('HTR3')
htr4=HeatExchanger('HTR4')
p3=Pump('HTR3 Pump')
m3=Merge('M3',num_in=2)
m4=Merge('M4',num_in=2)
m5=Merge('M5',num_in=2)
m6=Merge('M6',num_in=2)
m7=Merge('M7',num_in=2)
m8=Merge('M8',num_in=2)
d1=Node('HTR1 Drum',num_in=2,num_out=1)
d2=Node('HTR2 Drum',num_in=1,num_out=1)
d4=Node('HTR4 Drum',num_in=1,num_out=1)
c47=Connection(gc,'out1',htr1,'in2',label='47')
c48=Connection(htr1,'out2',htr2,'in2',label='48')
c49=Connection(htr2,'out2',htr3,'in2',label='49')
c50=Connection(htr3,'out2',m8,'in1',label='50')
c51=Connection(lp6,'out2',m3,'in1',label='51')
c52=Connection(lp7,'out2',m3,'in2',label='52')
c53=Connection(lp8,'out2',m4,'in1',label='53')
c54=Connection(lp9,'out2',m4,'in2',label='54')
c55=Connection(m3,'out1',m5,'in1',label='55')
c56=Connection(m4,'out1',m5,'in2',label='56')
c57=Connection(m5,'out1',htr1,'in1',label='57')
c58=Connection(htr1,'out1',d1,'in1',label='58')
c59=Connection(d1,'out1',hotwell,'in2',label='59')
c60=Connection(lp4,'out2',m6,'in1',label='60')
c61=Connection(lp5,'out2',m6,'in2',label='61')
c62=Connection(m6,'out1',htr2,'in1',label='62')
c63=Connection(htr2,'out1',d2,'in1',label='63')
c64=Connection(d2,'out1',d1,'in2',label='64')
c65=Connection(lp3,'out2',htr3,'in1',label='65')
c66=Connection(htr3,'out1',m7,'in1',label='66')
c67=Connection(m7,'out1',p3,'in1',label='67')
c68=Connection(p3,'out1',m8,'in2',label='68')
c69=Connection(m8,'out1',htr4,'in2',label='69')
c70=Connection(lp2,'out2',htr4,'in1',label='70')
c71=Connection(htr4,'out1',d4,'in1',label='71')
c92=Connection(d4,'out1',m7,'in2',label='92')

#deaerator
dea1=Node('Dea1',num_in=4,num_out=1)
m10=Merge('M10',num_in=2)
c72=Connection(hp21,'out2',m10,'in1',label='72')
c73=Connection(hp22,'out2',m10,'in2',label='73')
c74=Connection(m10,'out1',dea1,'in1',label='74')
c75=Connection(htr4,'out2',dea1,'in2',label='75')
c76=Connection(sp,'out1',dea1,'in3',label='76')


#feed pump
fp=Pump('Feed Pump')
c77=Connection(dea1,'out1',fp,'in1',label='77')

#high-pressure heat recovery heaters
htr6=HeatExchanger('HTR6')
htr7=HeatExchanger('HTR7')
m11=Merge('M11',num_in=2)
m12=Merge('M12',num_in=2)
d6=Node('Drum6',num_in=2,num_out=1)
d7=Node('Drum7',num_in=1,num_out=1)
c78=Connection(fp,'out1',htr6,'in2',label='78')
c79=Connection(htr6,'out2',htr7,'in2',label='79')
c80=Connection(hp3,'out2',m11,'in1',label='80')
c81=Connection(rh2,'out2',m11,'in2',label='81')
c82=Connection(m11,'out1',htr6,'in1',label='82')
c83=Connection(htr6,'out1',d6,'in1',label='83')
c84=Connection(d6,'out1',dea1,'in4',label='84')
c85=Connection(hp41,'out2',m12,'in1',label='85')
c86=Connection(rh1,'out2',m12,'in2',label='86')
c87=Connection(m12,'out1',htr7,'in1',label='87')
c88=Connection(htr7,'out1',d7,'in1',label='88')
c89=Connection(d7,'out1',d6,'in2',label='89')
c90=Connection(htr7,'out2',sg,'in1',label='90')

my_plant.add_conns(c1,c2,c3,c4,c5,c6,c7,c8,c9,c10,
c11,c12,c13,c14,c15,c16,c17,c18,c19,c20,
c21,c22,c23,c24,c25,c26,c27,c28,c29,c30,
c31,c32,c33,c34,c35,c36,c37,c38,c39,c40,
c41,c42,c43,c44,c45,c46,c47,c48,c49,c50,
c51,c52,c53,c54,c55,c56,c57,c58,c59,c60,
c61,c62,c63,c64,c65,c66,c67,c68,c69,c70,c71,c72,c73,c74,c75,c76,c77,c78,c79,c80,
c81,c82,c83,c84,c85,c86,c87,c88,c89,c90,
c91,c92,c93)


c1.set_attr(T=650,m=600,p=21,fluid={'H2O':1.0})
c42.set_attr(T=25,p=0.15, fluid={'water': 1.0})
c43.set_attr(T=40)
cd.set_attr(pr1=1, pr2=0.98)
sg.set_attr(pr=0.9)
fp.set_attr(eta_s=0.75)
cp.set_attr(eta_s=0.75)
p3.set_attr(eta_s=0.75)
#cd.set_attr(ttd_u=15)
hptu1.set_attr(eta_s=0.90)
hptu2.set_attr(eta_s=0.90)
hptu3.set_attr(eta_s=0.90)
hptu4.set_attr(eta_s=0.90)
hptu5.set_attr(eta_s=0.90)
hptu6.set_attr(eta_s=0.90)
lptu1.set_attr(eta_s=0.90)
lptu2.set_attr(eta_s=0.90)
lptu3.set_attr(eta_s=0.90)
lptu4.set_attr(eta_s=0.90)
lptu5.set_attr(eta_s=0.90)
lptu6.set_attr(eta_s=0.90)
lptu7.set_attr(eta_s=0.90)
lptu8.set_attr(eta_s=0.90)
lptu9.set_attr(eta_s=0.90)
lptu10.set_attr(eta_s=0.90)
c14.set_attr(fluid = {'H2O': 1.0},p = 15,x=0.9)
c15.set_attr(x=1)
c4.set_attr(p=20,T=620)
c5.set_attr(p=20,T=620)
c91.set_attr(m=500)
c2.set_attr(m=250)
c6.set_attr(m=225)
c7.set_attr(m=225)
c10.set_attr(m=200)
c11.set_attr(m=200)
c18.set_attr(m=15)
c21.set_attr(m=190)
c25.set_attr(m=180)
c26.set_attr(m=180)
c29.set_attr(m=170)
c30.set_attr(m=170)
c33.set_attr(m=160)
c34.set_attr(m=160)
c36.set_attr(m=150)
c37.set_attr(m=150)
gc.set_attr(pr=0.9)
htr1.set_attr(pr2=0.9,pr1=0.9)
rh1.set_attr(pr1=0.9)
rh2.set_attr(pr1=0.9)
htr2.set_attr(pr2=0.9,pr1=0.9)
htr3.set_attr(pr2=0.9,pr1=0.9)
htr4.set_attr(pr2=0.9,pr1=0.9)
htr6.set_attr(pr2=0.9)
htr7.set_attr(pr2=0.9)
c40.set_attr(p=0.2,T=120)
c20.set_attr(T=650)
c90.set_attr(T=300)
c75.set_attr(T=150)
c41.set_attr(T=45)



my_plant.set_attr(iterinfo=True)
my_plant.solve('design')
my_plant.print_results()

Please help me improve the parameter settings. The final data doesn't matter, it just needs to be runnable.Thank you so much！

[fwitte]

Hi @WangJ501,

this is asking for a bit much honestly as the model is quite complex. I'd advise you to build it incrementally to understand what specifications you have to make.

Best

[WangJ501]

Okay, thank you.