"Ciclo de Carnot"
Temp_H=300 [C]
Temp_L=30 [C]
"Ponto 1"
T[1]=Temp_H
x[1]=1
Pressao[1]=Pressure(Steam_IAPWS;T=T[1];x=1)
s[1]=Entropy(Steam_IAPWS;T=T[1];x=x[1])
h[1]=Enthalpy(Steam_IAPWS;T=T[1];x=x[1])
"Ponto 2"
T[2]=Temp_L
s[2]=s[1]
Pressao[2]=Pressure(Steam_IAPWS;T=T[2];s=s[2])
h[2]=Enthalpy(Steam_IAPWS;T=T[2];s=s[2])
x[2]=Quality(Steam_IAPWS;T=T[2];s=s[2])
"Ponto 3"
T[3]=Temp_L
s[3]=s[4]
Pressao[3]=Pressure(Steam_IAPWS;T=T[3];s=s[3])
h[3]=Enthalpy(Steam_IAPWS;T=T[3];s=s[3])
x[3]=Quality(Steam_IAPWS;T=T[3];s=s[3])
"Ponto 4"
T[4]=Temp_H
s[4]=Entropy(Steam_IAPWS;T=T[4];x=x[4])
Pressao[4]=Pressure(Steam_IAPWS;T=T[4];s=s[4])
h[4]=Enthalpy(Steam_IAPWS;T=T[4];s=s[4])
x[4]=0
"Eff_carnot_max=1-((TL+273)/(TH+273))"
eff_ciclo=Wciclo/QH
"Primeira Lei caldeira "
h[1]=QH-h[4]
"Primeira Lei turbina"
WT+h[2]=h[1]
"Primeira Lei condensador"
h[3]=Ql+h[2]
"Primeira Lei bomba"
WB + h[4] = h[3]
Trabalho_ciclo=QH-abs(QL)
Trabalho_ciclo2 = WT-abs(WB)