{Prova 3 - Termodinamica - Lucas Paladini}

W_dot_liq = 150000 [kW]

eta_Turbina = 0,8
eta_Turbina = (h_5-h_6)/(h_5-h_6s)
eta_Turbina = (h_5-h_7)/(h_5-h_7s)

s_5 = s_6s
s_5 = s_7s

eta_Bomba = 0,75
eta_Bomba = (h_3-h_4s)/(h_3-h_4)
eta_Bomba = (h_1-h_2s)/(h_1-h_2)

s_3 = s_4s
s_1 = s_2s

T_F_turbina = 350 [k]
T_F_bomba = 350 [K]
T_F_caldeira = 1300 [K]
T_F_condensador = 300 [K]

0 = -10/T_F_turbina + m_dot_5*s_5 - m_dot_6*s_6 + sigma_T1
0 = -20/T_F_turbina + m_dot_5*s_5 - m_dot_7*s_7 + sigma_T2
0 = Q_dot_cond/T_F_condensador + m_dot_7*s_7 - m_dot_1*s_1 + sigma_Cond
0 = -1/T_F_bomba + m_dot_1*(s_1-s_2) + sigma_B1
0 = 0 + (m_dot_2*s_2 + m_dot_6*s_6) - (m_dot_3*s_3) + sigma_Reg
0 = -1/T_F_bomba + m_dot_3*s_3 - m_dot_4*s_4 + sigma_B2
0 = Q_dot_caldeira/T_F_caldeira + m_dot_4*s_4 - m_dot_5*s_5 + sigma_Cald

0 = -10 - W_dot_T1 + m_dot_5*h_5 - m_dot_6*h_6
0 = -20 - W_dot_T2 + m_dot_5*h_5 - m_dot_7*h_7

0 = 0 + (m_dot_2*h_2 + m_dot_6*h_6) - (m_dot_3*h_3)

0 = Q_dot_cond + m_dot_7*h_7 - m_dot_1*h_1
0 = -1 - W_dot_B1 + m_dot_1*h_1 - m_dot_2*h_2
0 = -1 - W_dot_B2 + m_dot_3*h_3 - m_dot_4*h_4
0 = Q_dot_caldeira + m_dot_4*h_4 - m_dot_5*h_5

W_dot_T1 + W_dot_T2 + W_dot_B1 + W_dot_B2 = W_dot_liq

P_5 = 10000 [kPa]
T_5 = 773 [K]
s_5 = Entropy(Water;p=p_5;t=t_5)
h_5 = Enthalpy(Water;p=p_5;t=t_5)
v_5 = Volume(Water;p=p_5;t=t_5)

P_6 = 500 [kPa]
h_6s = Enthalpy(Water;p=p_6;s=s_6s)
s_6 = Entropy(Water;p=p_6;h=h_6)

P_7 = 10 [kPa]
h_7s = Enthalpy(Water;p=p_7;s=s_7s)
s_7 = Entropy(Water;p=p_7;h=h_7)

x_1 = 0
p_1 = p_7
h_1 = Enthalpy(Water;p=p_1;x=x_1)
s_1 = Entropy(Water;p=p_1;x=x_1)

s_2 = Entropy(Water;p=p_2;h=h_2)

x_3 = 0
p_3 = p_6
p_3 = p_2
h_3 = Enthalpy(Water;x=x_3;p=p_3)
s_3 = Entropy(Water;x=x_3;p=p_3)

p_4 = p_5

m_dot_5 = m_dot_6 + m_dot_7
m_dot_6 + m_dot_2 = m_dot_3

m_dot_4 = m_dot_5
m_dot_3 = m_dot_4
m_dot_7 = m_dot_1
m_dot_1 = m_dot_2