> restart:#"m17_p49"

En una central de regasificación de gas natural licuado, se va a usar un ciclo Brayton de nitrógeno en circuito cerrado para generar 1 MW de potencia mecánica. El compresor, de rendimiento isentrópico 0,8 y relación de presiones 4, toma el nitrógeno a 4 MPa y -100 ºC, y lo hace pasar por un recuperador de calor con los gases de salida de la turbina a contracorriente. A la salida del intercambiador, el gas recibe calor de la combustión exterior de gas natural hasta alcanzar los 800 ºC antes de entrar en la turbina, de rendimiento isentrópico 0,9. Tras la turbina, el gas pasa por el recuperador, y finalmente por un cambiador de calor donde una corriente de gas natural licuado lo enfría hasta las condiciones de entrada al compresor. Suponiendo para el intercambiador regenerativo un rendimiento del 70%, s pide:

a) Temperaturas de salida del compresor y de la turbina.

b) Gasto másico circulante, y temperatura a la que se calienta el fluido de trabajo en el intercambiador.

c) Gasto de gas natural necesario para alcanzar las condiciones de entrada a la turbina.

d) Calor transmitido al gas natural licuado.

e) Rendimiento térmico del ciclo motor.

Datos:

> read`../therm_eq.m`:read`../therm_const.m`:read`../therm_proc.m`:with(therm_proc):Digits:=5:with(plots):

> su1:="N2":dat:=[Wdot=1e6*W_,eta[C]=0.8,pi[12]=4,p1=4e6*Pa_,T1=(-100+273.15)*K_,T3=(800+273.15)*K_,eta[T]=0.9,eta[R]=0.7];

[Wdot = `+`(`*`(0.1e7, `*`(W_))), eta[C] = .8, pi[12] = 4, p1 = `+`(`*`(0.4e7, `*`(Pa_))), T1 = `+`(`*`(173.15, `*`(K_))), T3 = `+`(`*`(1073.2, `*`(K_))), eta[T] = .9, eta[R] = .7]
[Wdot = `+`(`*`(0.1e7, `*`(W_))), eta[C] = .8, pi[12] = 4, p1 = `+`(`*`(0.4e7, `*`(Pa_))), T1 = `+`(`*`(173.15, `*`(K_))), T3 = `+`(`*`(1073.2, `*`(K_))), eta[T] = .9, eta[R] = .7]

> gdat:=get_gas_data(su1):ldat:=get_liq_data(su1):dat:=op(dat),gdat,ldat,Const,SI2,SI1:get_pv_data(su1):

Image

a) Temperaturas de salida del compresor y de la turbina.

> eq5_59;T2_:=solve(subs(eta=eta[C],dat,eq5_59),T2);eq5_60;T4_:=solve(subs(eta=eta[T],dat,eq5_60),T4);

eta = `/`(`*`(`+`(`^`(pi[12], `/`(`*`(`+`(gamma, `-`(1))), `*`(gamma))), `-`(1))), `*`(`+`(`/`(`*`(T2), `*`(T1)), `-`(1))))
`+`(`*`(278.25, `*`(K_)))
eta = `/`(`*`(`+`(1, `-`(`/`(`*`(T4), `*`(T3))))), `*`(`+`(1, `-`(`^`(`/`(1, `*`(pi[12])), `/`(`*`(`+`(gamma, `-`(1))), `*`(gamma)))))))
`+`(`*`(757.50, `*`(K_)))

i.e. del compresor sale a 278 K, y de la turbina  a 758 K.

b) Gasto másico circulante, y temperatura a la que se calienta el fluido de trabajo en el intercambiador.

> eqW:=Wdot=mdot*c[p]*(T3-T4-T2+T1);mdot_:=subs(dat,solve(subs(T2=T2_,T4=T4_,dat,eqW),mdot));eqR:=eta=(T2R-T2)/(T4-T2);T2R_:=solve(subs(eta=eta[R],T2=T2_,T4=T4_,dat,eqR),T2R);eqBER:=T2R-T2=T4-T4R;T4R_:=T4_-T2R_+T2_;

Wdot = `*`(mdot, `*`(c[p], `*`(`+`(T3, `-`(T4), `-`(T2), T1))))
`+`(`/`(`*`(4.5668, `*`(kg_)), `*`(s_)))
eta = `/`(`*`(`+`(T2R, `-`(T2))), `*`(`+`(T4, `-`(T2))))
`+`(`*`(613.72, `*`(K_)))
`+`(T2R, `-`(T2)) = `+`(T4, `-`(T4R))
`+`(`*`(422.08, `*`(K_)))

i.e. hay que hacer circular 4.6 kg/s de nitrógeno, que en el recuperador, por la parte de alta pasa de 278 K a 614 K., y por la parte de baja, pasa de 758 K a 422 K.

c) Gasto de gas natural necesario para alcanzar las condiciones de entrada a la turbina.

Tomamos para el gas natural un PCI=50 MJ/kg, que puede obtenerse de Fuel_data: PCS=55,5 MJ/kg=890 kJ/mol; CH4+2O2=CO2+2H2O, PCI=PCS-2·hLV=890-2·44=800 kJ/mol=50 MJ/kg.

> Qdot2R3:=mdot*c[p]*(T3-T2R);Qdot2R3_:=subs(T2R=T2R_,mdot=mdot_,dat,%);mdotGN:='Qdot2R3'/PCI;PCI:=50e6*J_/kg_;mdotGN_:=subs(dat,Qdot2R3_/PCI);

`*`(mdot, `*`(c[p], `*`(`+`(T3, `-`(T2R)))))
`+`(`*`(0.21823e7, `*`(W_)))
`/`(`*`(Qdot2R3), `*`(PCI))
`+`(`/`(`*`(0.50e8, `*`(J_)), `*`(kg_)))
`+`(`/`(`*`(0.43646e-1, `*`(kg_)), `*`(s_)))

i.e. se consumen 0,044 kg/s de gas natural para aportar los 2,2 MW tras el recuperador.

d) Calor transmitido al gas natural licuado.

> Qdot4R1:=mdot*c[p]*(T4R-T1);Qdot4R1_:=subs(mdot=mdot_,T4R=T4R_,dat,%);

`*`(mdot, `*`(c[p], `*`(`+`(T4R, `-`(T1)))))
`+`(`*`(0.11823e7, `*`(W_)))

i.e. se comunica 1,2 MW al gas natural licuado.

e) Rendimiento térmico del ciclo motor.

Si no contabilizamos la refrigeración, por ser una planta de regasificación, será:

> eqeta:=eta=Wdot/'Qdot2R3';eqeta_:=subs(dat,eta=Wdot/Qdot2R3_);

eta = `/`(`*`(Wdot), `*`(Qdot2R3))
eta = .45823

Es un rendimiento muy alto (los rendimientos prácticos de estos ciclos son más bien del 30% en vez del 46%) porque se ha supuesto que la refrigeración es gratis.

>