Untitled

p_b=150 [bar]

T_max= 500 [c]

p_reheat=8.806 [bar]

T_reheat= 500 [c]

p_c = 0.1 [bar]

p_OFWH=5[bar]

p_CFWH=3[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} 
 h[2]=((h_s[2]-h[1])/eta_pump)+h[1]

{point 3}
p[3]=p_CFWH
{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_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} 
 h[5]=((h_s[5]-h[4])/eta_pump)+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 8}
p[8]=p_reheat
s_s[8]=s[6]
h_s[8]=enthalpy(water,p=p[8],s=s_s[8])
{Actual} 
 h[8]=((h_s[8]-h[6])/eta_pump)+h[6]


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

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

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


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


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

h[3] = h[12]

{point 13} 
p[13]=p_c
h[13]=h[3]

{mass balance of OFWH}

(x)*(h[7])+(1-x)*(h[3])=h[4]


{mass balance of CFWH}

(y)*(h[10])+(1-x)*(h[2])=(y)*(h[12])+(1-x)*(h[3])


{Work net}
w_net=(h[6]-h[5])+(h[9]-h[8])-(1-x-y)*(h[11])+(y)*(h[13])+(1-x)*(h[1])


{efficiency}
{eta_th=w_net / q_add}
q_add=(h[6]-h[5])+(h[9]-h[8])*(1-x)


eta_th=w_net / q_add