Untitled

p_b=150 [bar]
T_max= 500 [c]

{p_reheat=10 [bar]}
T_reheat= 500 [c]

p_c = 0.1 [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_b
{isentropic}
s_s[2]=s[1]
h_s[2]=enthalpy(water,p=p[2],s=s_s[2])
{Actual} 
"h[2]=( (h_s[2]-h[1]) / eta_pump ) + h[1]"
eta_pump = ( h_s[2] - h[1] ) / ( h[2] - h[1] )

{point 3}
p[3]=p_b
T[3]=T_max

h[3]=enthalpy(water,p=p[3],T=T[3])
S[3]=entropy(water,p=p[3],T=T[3])


{point 4}
p[4]=p_reheat
{isentropic}
s_s[4]=s[3]
h_s[4]=enthalpy(water,p=p[4],s=s[3])
{Actual} 
"h[4]=h[3]-((eta_turbine)/h[3]-h_s[4])"
eta_turbine = ( h[5] - h[4] ) / ( h[5] - h_s[4] )

{point 5}
p[5]=p_reheat
T[5]=T[3]
h[5]=enthalpy(water,p=p[5],T=T[5]) 
s[5]=entropy(water,p=p[5],T=T[5])


{point 6} 
p[6] =p_c
{isentropic} 
s_s[6]=s[5]
h_s[6]=enthalpy(water,p=p[6],s=s_s[6])
{Actual} 
"h[6]=h[5]-((eta_turbine)/h[5]-h_s[6])"
eta_turbine = ( h[5] - h[6] ) / ( h[5] - h_s[6] )


{Work net}
w_net=((h[3]-h[4])+(h[5]-h[6]))-(h[2]-h[1])

{efficiency}
{eta_th=w_net / q_add}
q_add=(h[3]-h[2])+(h[5]-h[4])
eta_th=w_net / q_add