"Rankine cycle example with reheat and regeneration"
 
"Given"
P_cond = 10 
P_boil = 10000
P_RH = 6000
 
T_3 = 700
m_dot_B = 1.25 
 
"Open feedwater heater"
0 = m_dot_C + m_dot_9 - m_dot_B
0 = m_dot_C*h_8 + m_dot_9*h_9 - m_dot_B*h_1
 
 
h_1= enthalpy(Steam, P=P_RH, x=0)
h_9 = enthalpy(Steam, P=P_RH, s=s_3)
s_3= entropy(Steam, P=P_boil, T=T_3)
h_8 = enthalpy(Steam, P=P_RH, s=s_7)
s_7= entropy(Steam, P=P_cond, x=0)
 
"Pump 1"
W_dot_p_1 = m_dot_B*(h_2-h_1)
s_1= entropy(Steam, P=P_RH, x=0)
h_2= enthalpy(Steam, P=P_boil, s=s_1)
  
"Boiler"
Q_dot_b = m_dot_B*(h_3 - h_2)
h_3 = enthalpy(Steam, T=T_3,P=P_boil)
 
"HP Turbine"
W_dot_HPT = m_dot_B*(h_3 - h_4)
h_4= enthalpy(Steam,P=P_RH, s=s_3)
 
"Reheat"
Q_dot_RH = m_dot_C*(h_5 - h_4)
h_5 = enthalpy(Steam, T=T_3,P=P_RH)
 
"LP Turbine"
W_dot_LPT = m_dot_C*(h_5 - h_6)
h_6= enthalpy(Steam,P=P_cond, s=s_5)
s_5= entropy(Steam, P=P_RH, T=T_3)
 
"Condensor"
Q_dot_cond = m_dot_C*(h_6-h_7)
h_7= enthalpy(Steam, P=P_cond, x=0)
 
"Pump 2"
W_dot_p_2 = m_dot_C*(h_8-h_7)
 
"Cycle efficiency" 
eta_th = (W_dot_HPT+W_dot_LPT-W_dot_p_1-W_dot_p_2)/(Q_dot_b+Q_dot_RH)