> restart:#"m18_p20"

Una máquina refrigerante que utiliza R-12 se usa para enfriar una corriente de aire atmosférico de 0,1 kg/s desde las condiciones ambiente de 94 kPa, 30 ºC y 40% de humedad hasta 5 ºC. Se pide:

a) Hacer un esquema de la instalación y de las evoluciones del aire en el diagrama h-w y del R-12 en el p-h.
b) Determinar los estados del ciclo del R-12 sabiendo que los saltos térmicos a través de los intercambiadores han de ser como mínimo de 5 ºC.
c) Calcular la potencia necesaria para el compresor, así como el gasto circulante de R-12.
Datos:

> read"../therm_eq.m":read"../therm_proc.m":with(therm_proc):

> su1:="Aire":su2:="H2O":su3:="CCl2F2":dat:=[m=0.1*kg_/s_,p0=94e3*Pa_,T0=(30+273)*K_,phi0=0.4,Ts=(5+273)*K_,DT=5*K_,eta[C]=1];

`:=`(dat, [m = `+`(`/`(`*`(.1, `*`(kg_)), `*`(s_))), p0 = `+`(`*`(0.94e5, `*`(Pa_))), T0 = `+`(`*`(303, `*`(K_))), phi0 = .4, Ts = `+`(`*`(278, `*`(K_))), DT = `+`(`*`(5, `*`(K_))), eta[C] = 1])

Esquema:

> `:=`(Sistemas, [comp, conden, valv, evap])

> `:=`(Componentes, [Air, H2O, CCl2F2])

> `:=`(Estados, [1, 2, 3, 4, s])

> Adat:=get_gas_data(su1):Adat:=subs(c[p]=c[pa],R=R[a],M=M[a],T[b]=nada,[Adat]):dat:=op(dat),Const,SI2,SI1:

> Wgdat:=get_gas_data(su2):Wgdat:=subs(c[p]=c[pv],R=R[v],M=M[v],[Wgdat]):Wldat:=get_liq_data(su2):Wdat:=op(Wgdat),Wldat:

> R12gdat:=get_gas_data(su3):R12ldat:=get_liq_data(su3):R12dat:=R12gdat,R12ldat:get_pv_data(su3):pvR12(T):=pv(T):get_pv_data(su2):subs([c[pa],c[pv],c,T[b],h[lv0]])=subs(Adat,Wdat,[c[pa],c[pv],c,T[b],h[lv0]]);subs([c[p],c,T[b],h[lv0]])=subs(R12dat,[c[p],c,T[b],h[lv0]]);

[c[pa], c[pv], c, T[b], h[lv0]] = [`+`(`/`(`*`(1004., `*`(J_)), `*`(kg_, `*`(K_)))), `+`(`/`(`*`(1900., `*`(J_)), `*`(kg_, `*`(K_)))), `+`(`/`(`*`(4180., `*`(J_)), `*`(kg_, `*`(K_)))), `+`(`*`(373.2, ...

[c[p], c, T[b], h[lv0]] = [`+`(`/`(`*`(573., `*`(J_)), `*`(kg_, `*`(K_)))), `+`(`/`(`*`(966., `*`(J_)), `*`(kg_, `*`(K_)))), `+`(`*`(243.0, `*`(K_))), `+`(`/`(`*`(0.165e6, `*`(J_)), `*`(kg_)))]

Image

a) Hacer un esquema de la instalación y de las evoluciones del aire en el diagrama h-w y del R-12 en el p-h.

b) Determinar los estados del ciclo del R-12 sabiendo que los saltos térmicos a través de los intercambiadores han de ser como mínimo de 5 ºC.

> T1:=Ts-DT;T1_:=subs(dat,T1);p1_:=evalf(subs(T=T1_,pvR12(T))):'p1'=evalf(%,2);T3:=T0+DT;T3_:=subs(dat,T3);p2_:=evalf(subs(T=T3_,pvR12(T))):'p2'=evalf(%,2);h1_:=subs(R12dat,T=T1_,hv(T)):'h1'=evalf(%,2);;T2s_:=subs(dat,R12dat,T1_*(p2_/p1_)^((gamma-1)/gamma)):'T2s'=evalf(%,3);h2s_:=subs(dat,R12dat,T=T2s_,hv(T)):'h2s'=evalf(%,2);h2_:=subs(dat,h1_+(h2s_-h1_)/eta[C]):'h2'=evalf(%,2);h3_:=subs(dat,R12dat,T=T3_,hl(T)):'h3'=evalf(%,2);

`:=`(T1, `+`(Ts, `-`(DT)))

`:=`(T1_, `+`(`*`(273, `*`(K_))))

p1 = `+`(`*`(0.31e6, `*`(Pa_)))

`:=`(T3, `+`(T0, DT))

`:=`(T3_, `+`(`*`(308, `*`(K_))))

p2 = `+`(`*`(0.85e6, `*`(Pa_)))

h1 = `+`(`/`(`*`(0.42e6, `*`(J_)), `*`(kg_)))

T2s = `+`(`*`(308., `*`(K_)))

h2s = `+`(`/`(`*`(0.44e6, `*`(J_)), `*`(kg_)))

h2 = `+`(`/`(`*`(0.44e6, `*`(J_)), `*`(kg_)))

h3 = `+`(`/`(`*`(0.30e6, `*`(J_)), `*`(kg_)))

c) Calcular la potencia necesaria para el compresor, así como el gasto circulante de R-12.

> P[C]:=subs(dat,mR12*(h2-h1));eqBEevap:=mR12*(h1-h4)=m*(h0-hs);w0_:=evalf(subs(dat,Wdat,w(phi0,T0,p0))):'w0'=evalf(%,2);h0_:=subs(dat,Adat,Wdat,T=T0,dat,h(T,w0_)):'h0'=evalf(%,2);ws_:=evalf(subs(dat,Wdat,w(1,Ts,p0))):'ws'=evalf(%,2);hs_:=subs(dat,Adat,Wdat,T=Ts,dat,h(T,ws_)):'hs'=evalf(%,2);mR12_:=subs(dat,solve(subs(dat,mR12*(h1_-h3_)=m*(h0_-hs_)),mR12)):'mR12'=evalf(%,2);P[C]:=subs(dat,mR12_*(h2_-h1_)):'P[C]'=evalf(%,2);

`:=`(P[C], `*`(mR12, `*`(`+`(h2, `-`(h1)))))

`:=`(eqBEevap, `*`(mR12, `*`(`+`(h1, `-`(h4)))) = `*`(m, `*`(`+`(h0, `-`(hs)))))

w0 = 0.11e-1

h0 = `+`(`/`(`*`(0.59e5, `*`(J_)), `*`(kg_)))

ws = 0.59e-2

hs = `+`(`/`(`*`(0.20e5, `*`(J_)), `*`(kg_)))

mR12 = `+`(`/`(`*`(0.33e-1, `*`(kg_)), `*`(s_)))

P[C] = `+`(`*`(0.67e3, `*`(W_)))

i.e. el compresor consume 670 W.

>