Calculer le rendement d'une machine vapeur complexe

[nonoceb]

Bonjour, je cherche à evaluer la modification du surchauffeur sur un cycle de hirn à resurchauffe, comme celui ci :
Pièce jointe 472902

En gros je cherche à ajouter une surchauffe de 285°c 70bar à 500°c 70bar, uniquement sur la branche qui alimente la turbine HP, sans surchauffer SH, la branche du surchauffeur BP qui reste à 285°c.
J'ai biensur les enthalpies pour fabriquer ma vapeur surchauffée à 500°c mais j'ai du mal à evaluer comment cette vapeur va travailler et son influence sur le rendement, idealement j'aimerai pouvoir trouver un logiciel facile pour implementer ce genre de modification.
La librairie pyromat.org a implementer une demo pour le cycle de hirn mais on a pas les sources pour le modifier...et le get's started est trop light pour mon faible niveau.

Je cherche surtout à obtenir le rendement, apres la modification. Pas facile mais si vous avez une idée, merci
-----

[gienas]

Bonjour nonoceb et tout le groupe et bonne année 2023

Bienvenue sur le forum.

… je cherche à evaluer la modification …

La pièce jointe n’est pas passée.

Elle doit être présentée à nouveau dans une nouvelle réponse.

[nonoceb]

Bonjour et bonne année!
Je repost le diagramme

[nonoceb]

J'ai tenté ce script python (adriaan-van-niekerk sur medium com)mais je n'arrive pas à surchauffer de la vapeur, je reste en vapeur saturée

import pyromat as pm 
pm.config["unit_pressure"] = "bar"
pm.config["def_p"] = 100 
mp_water = pm.get("mp.H2O") 

#saturated liquid, thus x = 0
p1 = 10 # ←given
T1 = mp_water.Ts(p=p1)[0]
s1 = mp_water.ss(p=p1)[0]
p2 = 70 # ←given and converted to kPa
s2= s1
v = 1/mp_water.ds(p=p1)[0]
w_p = v*(p2-p1)
print(f"Work required by pump: {round(float(w_p),1)} kJ/kg")
h1 = mp_water.hs(p=p1)[0]
h2 = h1+w_p
T2 = mp_water.T_h(p=p2,h=h2)
print(f"h2 = {round(float(h2),1)} kJ/kg")

# steam leaves the boiler as saturated vapor, thus x = 1
p3 = p2
T3 = mp_water.Ts(p=p3)
h3 = mp_water.hs(p=p3)[1]
s3dash = mp_water.ss(p=p3)[0]
T3dash = T3
s3 = mp_water.ss(p=p3)[1]
q_H = h3-h2
print(f"Heat input by boiler: {round(float(q_H),1)} kJ/kg")
print(f"Temperature boiler: {T3} k")
print(f"Pressure boiler: {p3} bar")

p4 = p1
s4 = s3
T4, x = mp_water.T_s(s=s4,p=p4, quality=True)
h4 = mp_water.h(p=p4,x=x)
w_t = h3-h4
print(f"Quality of low pressure steam: {round(float(x),4)}")
print(f"Temperature : {T4} k")
print(f"Pressure : {p4} bar")
print(f"Work generated by turbine: {round(float(w_t),1)} kJ/kg")

q_L = h4-h1
print(f"Heat rejected by the condenser: {round(float(q_L),1)} kJ/kg")

eta_th = (w_t-w_p)/q_H*100
print(f"Thermal efficiency is: {round(float(eta_th),1)}%")

[A voir en vidéo sur Futura]