> | restart:#"m08_p63" |
Para un acondicionamiento de aire en verano, se quiere pasar 1 kg/s de aire ambiente, que está a 36 ºC y 30% HR, hasta 20 ºC y 60% HR. Para ello se va a dividir el flujo de aire en dos partes, enfriando una de ellas con una máquina frigorÃfica y mezclándolas después. Se pide:
a) Hacer un esquema de los procesos que sufre el aire, en un diagrama h-w.
b) Determinar las temperaturas de rocÃo y de bulbo húmedo del ambiente.
c) Determinar la humedad absoluta y la entalpÃa del aire antes y después de acondicionarlo.
d) Calcular la temperatura a la que hay que enfriar el aire en el refrigerador y el gasto másico de aire que ha de pasar por él.
e) Determinar el calor a evacuar en el refrigerador.
Datos:
read"../therm_eq.m":read"../therm_proc.m":with(therm_proc):Digits:=4:
> | su1:="Aire":su2:="H2O":dat:=[m=1*kg_/s_,T0=(36+273)*K_,phi0=0.3,T3=(20+273)*K_,phi3=0.6]; |
![]() |
Eqs. const.:
> | Adat:=get_gas_data(su1):Adat:=subs(c[p]=c[pa],R=R[a],M=M[a],T[b]=nada,[Adat]):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:get_pv_data(su2):dat:=op(dat),Const,SI2,SI1: |
a) Hacer un esquema de los procesos que sufre el aire, en un diagrama h-w.
(Ver arriba)
b) Determinar las temperatura de rocÃo y de bulbo húmedo del ambiente.
> | eq8_9;pv0:=subs(dat,evalf(subs(dat,pv(T0))));Tdew__:=solve(pv(Td)=phi0*pv(T0),Td):Tdew_:=subs(dat,evalf(subs(dat,Tdew__[2])));'Tdew'=TKC(%);eq8_10;Twet_:=fsolve(subs(dat,Adat,Wdat,SI0,c[pa]*(T0-T0)+w(phi0,T0,p0)*(h[lv0]-(c[pv]-c)*(T[b]-T[f]))=c[pa]*(T-T0)+w(1,T,p0)*(h[lv0]-(c[pv]-c)*(T[b]-T[f]))),T=200..400)*K_;'Twet'=TKC(%); |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
i.e. TrocÃo=15 ºC y Thúmeda=22 ºC.
c) Determinar la humedad absoluta y la entalpÃa del aire antes y después de acondicionarlo.
> | eq8_8;pv0_:=subs(dat,evalf(subs(dat,pv(T0))));pv2_:=subs(dat,evalf(subs(dat,pv(T2))));w0_:=subs(dat,evalf(subs(dat,w(phi0,T0,p0)))):'w0'=evalf(%,2);w3_:=subs(dat,evalf(subs(dat,w(phi3,T3,p0)))):'w3'=evalf(%,2);eq8_11;h0_:=subs(Adat,Wdat,T=T0,dat,h(T,w0_));h3_:=subs(Adat,Wdat,T=T3,dat,h(T,w3_)); |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
i.e. entra con 11 g/kg y 65 kJ/kg, y debe salir con 9 g/kg y 42 kJ/kg.
d) Calcular la temperatura a la que hay que enfriar el aire en el refrigerador.
A la salida estará saturado porque estará condensando.
Sea lambda la fracción refrigerada (junto a T2 son las incógnitas, que se obtienen de los balance en el mezclador):
> | eqBMw:=lambda*ma*w2+(1-lambda)*ma*w0=ma*w3;eqBE:=lambda*ma*h2+(1-lambda)*ma*h0=ma*h3;w2:=w(1,T2,p0);h2:=subs(Adat,Wdat,T=T2,dat,h(T,w2));sol:=fsolve(subs(w0=w0_,w3=w3_,h0=h0_,h3=h3_,dat,SI0,expand({eqBMw/ma,eqBE/ma})),{lambda,T2},{lambda=0..1,T2=250..350}):T2_:=subs(sol,T2)*K_;T2=TKC(%);lambda_:=evalf(subs(sol,lambda),2);h2__:=subs(dat,evalf(subs(T2=T2_,dat,h2))); |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
i.e. hay que refrigerar hasta 282 K (9 ºC) una fracción del 59% (0,59 kg/s).
Con el diagrama h-w bastarÃa haber trazado la recta que une los puntos 1 y 3; el corte con phi=1 da directamente el punto 2 (w2 y h2)´. Lambada se calcularÃa con la regla de la palanca entre los puntos 1, 3 y 2..
Se puede hacer por aproximaciones en función de T2 (con T2<T3), hasta que los dos balances den la misma lambda:
> | N:=10:arr:=array(1..N,1..3):for i from 1 to N do T2:=subs(dat,T3-(T3-T[tr])*(i/N));w2_:=subs(dat,w(1,T2,p0));h2_:=subs(Adat,Wdat,T=T2,dat,h(T,w2));lambda1:=(w3_-w0_)/(w2_-w0_);lambda2:=(h3_-h0_)/(h2_-h0_);print(`T=`,evalf(T2,3),`lambda=`,evalf(lambda1,2),`Delta=`,evalf(lambda1-lambda2,2));arr[i,1]:=T2;arr[i,2]:=lambda1;arr[i,3]:=lambda2;od:plot(subs(SI0,[[seq([arr[i,1],arr[i,2]],i=1..N)],[seq([arr[i,1],arr[i,3]],i=1..N)]]),273..293,0..1); |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
e) Determinar el calor a evacuar en el refrigerador.
> | h2:='h2':eqBE:=QR=ma*lambda*(h0-h2);eqBE_:=subs(ma=m,h0=h0_,h2=h2__,sol,dat,%); |
![]() |
![]() |
i.e. hay que extraer 23 kW. Nótese que este resultado podrÃa haberse obtenido directamente de los datos, sin resolver el problema de mezclado, haciendo el balance energético global:
> | eqBE:=QR=ma*(h0-h3);eqBE_:=subs(ma=m,h0=h0_,h3=h3_,dat,%); |
![]() |
![]() |
Todo el problema se podrÃa haber resuelto gráficamente mediante el diagrama de propiedades del aire húmedo (e.g. la fracción de gasto a refrigerar serÃa el cociente de distancias <13>/<12>
> |