>	restart:#"m08_p21"

Determinar el calor a extraer de una corriente de aire húmedo de 0,5 kg/s a 30 C, 95 kPa y 90% de humedad, para que salga a 10 C.

Datos:

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

>	su1:="Aire":su2:="H2O":dat:=[m=0.5*kg_/s_,T0=(30+273)*K_,p0=95e3*Pa_,phi0=0.9,T2=(10+273)*K_];

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) Determinar el calor a extraer de una corriente de aire húmedo de 0,5 kg/s a 30 C, 95 kPa y 90% de humedad, para que salga a 10 C.

Previsiblemente condensará, luego:

>	w0_:=evalf(subs(dat,w(phi0,T0,p0)));w1:=w0;h0_:=subs(dat,Adat,Wdat,T=T0,dat,h(T,w0_));h1:=h0;T0dew_:=evalf(subs(dat,solve(pv(T)=phi0*pv(T0),T)));T0dew=TKC(%);phi2:=1;w2_:=evalf(subs(dat,w(phi2,T2,p0)));h2_:=subs(dat,Adat,Wdat,T=T2,dat,h(T,w2_));Q_:=subs(dat,m*(h2_-h0_));

Efectivamente, comienza a condensar enseguida, a Tdew=28 ºC, y hay que refrigerar con 33 kW (nótese que si no condensara, bastarían Q=m·cpa·DT=0,5·1·20=10 kW.

>