Untitled

p_b=150 [bar]

T_max= 500 [c]

p_reheat=30.5 [bar]

m_OFWH = 0.1658 

p_c = 0.1 [bar]

p_OFWH=5.876[bar]

{p_CFWH=1.2[bar]}

eta_turbine=0.85
eta_pump=0.95

{point 1} 
p[1]=p_c 
x[1]=0
h[1]=enthalpy(water,p=p[1],x=x[1])
s[1]=entropy(water,p=p[1],x=x[1])

{point 2} 
p[2]=p_OFWH
{isentropic}
s_s[2]=s[1]
h_s[2]=enthalpy(water,p=p[2],s=s_s[2])
{Actual} 
eta_pump=(h_s[2]-h[1])/(h[2]-h[1])


{point 3}
{mass balance}
{h[3]=h[12]}

{point 4}
p[4]=p_OFWH
x[4]=0
h[4]=enthalpy(water,p=p[4],x=x[4])
s[4]=entropy(water,p=p[4],x=x[4])


{point 5}
p[5]=p_b
s_s[5]=s[4]
h_s[5]=enthalpy(water,p=p[5],s=s_s[5])
{Actual} 
eta_pump=(h_s[5]-h[4])/(h[5]-h[4])

{point 6} 
p[6] =p_b
T[6]=T_max 
h[6]=enthalpy(water,p=p[6],T=T[6])
s[6]=entropy(water,p=p[6],T=T[6])

{point 7}
p[7]=p_reheat
s_s[7]=s[6]
h_s[7]=enthalpy(water,p=p[7],s=s_s[7])
{Actual} 
 eta_turbine= (h[6]-h[7])/(h[6]-h_s[7])


{point 8} 
p[8] =30.5
T[8]=500
h[8]=enthalpy(water,p=p[8],T=T[8])
s[8]=entropy(water,p=p[8],T=T[8])


{point 9}
p[9]=p_OFWH
s_s[9]=s[8]
h_s[9]=enthalpy(water,p=p[9],s=s_s[9])
s[9]=entropy(water,p=p[9],h=h[9])
{Actual} 
 eta_turbine=(h[8]-h[9])/(h[8]-h_s[9])


{point 10} 
p[10] =p_CFWH
s_s[10]=s[9]
h_s[10]=enthalpy(water,p=p[10],s=s_s[10])
s[10]=entropy(water,p=p[10],h=h[10])
{Actual} 
eta_turbine=(h[9]-h[10])/(h[9]-h_s[10 ])


{point 11} 
p[11] =p_c
s_s[11]=s[10]
h_s[11]=enthalpy(water,p=p[11],s=s_s[11])
{Actual} 
 eta_turbine=(h[10]-h[11 ])/(h[10 ]-h_s[11 ])

{point 12}
p[12]=p_CFWH
x[12]=0
h[12]=enthalpy(water,p=p[12],x=x[12])

{point 13}
h[12]=h[13]

{mass balance of OFWH}

(m_OFWH)*(h[9])+(1-m_OFWH)*(h[3])=h[4]

{mass balance of closed feed water}
((1-m_OFWH)*(h[2]))+(m_CFWH*h[10])=(h[3]*1-m_OFWH)+(h[12]*m_CFWH)


{Work net}

w_net = q_add - q_rej
q_add = (h[6]-h[5])+(h[8]-h[7])
q_rej = ((h[13]-h[1])*(m_CFWH))+((h[11]-h[1])*(1-m_OFWH-m_CFWH))

{efficiency}

eta_th= 1-(q_rej/q_add)