> restart:#"m17_p05"

Se quiere diseñar una central de ciclo combinado Brayton-Rankine para una potencia total de 10 MW. El compresor de la turbina de gas tiene un rendimiento de 0,85. La presión y temperatura de entrada a la turbina de la turbina de gas, cuyo rendimiento es 0,85, son de 1 MPa y 1300 K. Los gases de escape salen de la caldera de vapor a 250 °C. La presión y temperatura de entrada a la turbina de vapor, cuyo rendimiento es 0,80, son de 3,5 MPa y 350 °C, y la presión de salida de 10 kPa. Se pide:

a) Hacer un esquema de la instalación.

b) Potencia específica de la turbina de gas.

c) Gasto de agua y gasto de aire.

d) Calor recibido del combustible, calor intercambiado en la caldera y rendimiento energético global.

> read`../therm_eq.m`:read`../therm_const.m`:read`../therm_proc.m`:with(therm_proc):

> su1:="Aire":su2:="H2O":dat:=[Pe=10e6*W_,eta[CG]=0.85,eta[TG]=0.85,T3TG=1300*K_,p3TG=1e6*Pa_,T5TG=(250+273)*K_,T3TV=(350+273)*K_,p3TV=3.5e6*Pa_,eta[TV]=0.8,p4TV=10e3*Pa_];

Typesetting:-mprintslash([dat := [Pe = `+`(`*`(0.10e8, `*`(W_))), eta[CG] = .85, eta[TG] = .85, T3TG = `+`(`*`(1300, `*`(K_))), p3TG = `+`(`*`(0.1e7, `*`(Pa_))), T5TG = `+`(`*`(523, `*`(K_))), T3TV = ...
Typesetting:-mprintslash([dat := [Pe = `+`(`*`(0.10e8, `*`(W_))), eta[CG] = .85, eta[TG] = .85, T3TG = `+`(`*`(1300, `*`(K_))), p3TG = `+`(`*`(0.1e7, `*`(Pa_))), T5TG = `+`(`*`(523, `*`(K_))), T3TV = ...

Image

> eqBE:=eq5_43;Wgdat:=get_gas_data(su2):Wldat:=get_liq_data(su2):get_pv_data(su2):Adat:=get_gas_data(su1):dat:=op(dat),Wgdat,Wldat,Const,SI2,SI1:

Typesetting:-mprintslash([eqBE := Dh[t] = `+`(w, q)], [Dh[t] = `+`(w, q)])

a) Hacer un esquema de la instalación.
(arriba.)

b) Potencia específica de la turbina de gas

> wTG:=c[p]*(T3-T4)-c[p]*(T2-T1);eq5_59;T2_:=subs(eta=eta[CG],T1=T0,pi[12]=p3TG/p0,Adat,dat,solve(eq5_59,T2)):'T2'=evalf(%,3);eq5_60;T4_:=subs(eta=eta[TG],T3=T3TG,p3=p3TG,p4=p0,Adat,dat,solve(eq5_60,T4)):'T4'=evalf(%,3);T4TG:=T4_:wTG_:=subs(Adat,dat,c[p]*(T3TG-T4_)-c[p]*(T2_-T0)):'wTG'=evalf(%,3);qTG_:=subs(Adat,dat,c[p]*(T3TG-T2_)):'qTG'=evalf(%,3);

Typesetting:-mprintslash([wTG := `+`(`*`(c[p], `*`(`+`(T3, `-`(T4)))), `-`(`*`(c[p], `*`(`+`(T2, `-`(T1))))))], [`+`(`*`(c[p], `*`(`+`(T3, `-`(T4)))), `-`(`*`(c[p], `*`(`+`(T2, `-`(T1))))))])
eta = `/`(`*`(`+`(`^`(pi[12], `/`(`*`(`+`(gamma, `-`(1))), `*`(gamma))), `-`(1))), `*`(`+`(`/`(`*`(T2), `*`(T1)), `-`(1))))
T2 = `+`(`*`(603., `*`(K_)))
eta = `/`(`*`(`+`(1, `-`(`/`(`*`(T4), `*`(T3))))), `*`(`+`(1, `-`(`^`(`/`(`*`(p4), `*`(p3)), `/`(`*`(`+`(gamma, `-`(1))), `*`(gamma)))))))
T4 = `+`(`*`(768., `*`(K_)))
wTG = `+`(`/`(`*`(0.218e6, `*`(J_)), `*`(kg_)))
qTG = `+`(`/`(`*`(0.699e6, `*`(J_)), `*`(kg_)))

i.e. la potencia unitaria neta es de 218 kJ/kg, y calor recibido 700 kJ/kg.

c) Gasto de agua y gasto de air

TV:

> T1_:=subs(dat,evalf(fsolve(subs(dat,SI0,p4TV=pv(T)),T=100..1000)))*K_:'T1'=evalf(%,3);h1_:=subs(dat,T=T1_,dat,hl(T)):'h1'=evalf(%,3);h3_:=subs(dat,T=T3TV,dat,hv(T)):'h3'=evalf(%,3);wet_vap(subs(dat,T3TV),subs(dat,p3TV),subs(dat,p4TV));h4_:=subs(dat,h3_-(h3_-h4s_)*eta[TV]):'h4'=evalf(%,3);wTV_:=h3_-h4_:'wTV'=evalf(%,3);qTV_:=h3_-h1_:'qTV'=evalf(%,3);eqBEcarga:=ma*wTG+mw*wTV=P;eqBEHX:=ma*c[p]*(T4-T5)=mw*qTV;eqBEcarga_:=subs(dat,ma*wTG_+mw*wTV_=Pe):'eqBEcarga'=evalf(%,2);eqBEHX_:=subs(Adat,dat,ma*c[p]*(T4TG-T5TG)=mw*qTV_):'eqBEHX'=evalf(%,2);sol1_:=subs(dat,solve({eqBEcarga_,eqBEHX_},{ma,mw})):evalf(%,2);

T1 = `+`(`*`(319., `*`(K_)))
h1 = `+`(`/`(`*`(0.192e6, `*`(J_)), `*`(kg_)))
h3 = `+`(`/`(`*`(0.319e7, `*`(J_)), `*`(kg_)))
`Point 4s inside wet bell`, [x4s_ = .82, h4s_ = `+`(`/`(`*`(0.21e7, `*`(J_)), `*`(kg_))), s3_ = `+`(`/`(`*`(0.68e4, `*`(J_)), `*`(kg_, `*`(K_))))]
h4 = `+`(`/`(`*`(0.235e7, `*`(J_)), `*`(kg_)))
wTV = `+`(`/`(`*`(0.838e6, `*`(J_)), `*`(kg_)))
qTV = `+`(`/`(`*`(0.300e7, `*`(J_)), `*`(kg_)))
Typesetting:-mprintslash([eqBEcarga := `+`(`*`(ma, `*`(`+`(`*`(c[p], `*`(`+`(T3, `-`(T4)))), `-`(`*`(c[p], `*`(`+`(T2, `-`(T1)))))))), `*`(mw, `*`(wTV))) = P], [`+`(`*`(ma, `*`(`+`(`*`(c[p], `*`(`+`(T...
Typesetting:-mprintslash([eqBEHX := `*`(ma, `*`(c[p], `*`(`+`(T4, `-`(T5))))) = `*`(mw, `*`(qTV))], [`*`(ma, `*`(c[p], `*`(`+`(T4, `-`(T5))))) = `*`(mw, `*`(qTV))])
eqBEcarga = (`+`(`/`(`*`(0.22e6, `*`(ma, `*`(`^`(m_, 2)))), `*`(`^`(s_, 2))), `/`(`*`(0.84e6, `*`(mw, `*`(`^`(m_, 2)))), `*`(`^`(s_, 2)))) = `+`(`*`(0.10e8, `*`(W_))))
eqBEHX = (`+`(`/`(`*`(0.25e6, `*`(ma, `*`(`^`(m_, 2)))), `*`(`^`(s_, 2)))) = `+`(`/`(`*`(0.30e7, `*`(mw, `*`(`^`(m_, 2)))), `*`(`^`(s_, 2)))))
{ma = `+`(`/`(`*`(35., `*`(kg_)), `*`(s_))), mw = `+`(`/`(`*`(2.9, `*`(kg_)), `*`(s_)))}

i.e. han de circular 35 kg/s de aire por la TG y 2,9 kg/s de agua por la TV.

d) Calor recibido del combustible, calor intercambiado en la caldera y rendimiento energético global.

> Qcomb:=ma*qTG;Qcomb_:=subs(sol1_,dat,ma*qTG_):'Qcomb'=evalf(%,2);QHX:=mw*qTV;QHX_:=subs(sol1_,dat,mw*qTV_):'QHX'=evalf(%,2);eta:='Pe/Qcomb';eta_:=subs(dat,Pe/Qcomb_):'eta'=evalf(%,2);

Typesetting:-mprintslash([Qcomb := `*`(ma, `*`(qTG))], [`*`(ma, `*`(qTG))])
Qcomb = `+`(`*`(0.24e8, `*`(W_)))
Typesetting:-mprintslash([QHX := `*`(mw, `*`(qTV))], [`*`(mw, `*`(qTV))])
QHX = `+`(`*`(0.86e7, `*`(W_)))
Typesetting:-mprintslash([eta := `/`(`*`(Pe), `*`(Qcomb))], [`/`(`*`(Pe), `*`(Qcomb))])
eta = .41

i.e. en la cámara de combustión se reciben 24 MW, y en el intercambiador de los ciclos se aprovechan 8,6 MW- El rendimiento global es del 41 %.

>