>	restart:#"m08_p39"

Se trata del anlisis energtico de un secadero de malta (cebada germinada a remojo, donde el almidón se ha transformado en azcar, que, una vez secada en caliente, i.e. tostada, se usa en la fabricación de cerveza y whisky). Aunque realmente suele hacerse el secado en tandas con bandejas perforadas, supóngase aqu que se hace en continuo soplando aire a 70 ºC a contracorriente en un secadero en el que entran 150 toneladas/da de producto (malta) a 17 ºC y con 45% de agua en peso, saliendo prcticamente a 70 ºC y con 4% de agua. El aire a 70 ºC se obtiene calentando aire desde las condiciones ambiente: 15 ºC, 93 kPa y 60 % de humedad. Se pide
a) Estimar la potencia trmica que requiere la evaporación del agua del producto, y compararla con la requerida para calentar el producto seco suponiendo que su capacidad trmica es de 2000 J/(kg×K).
b) Establecer el balance msico de H2O y el balance energtico del secadero despreciando el posible aporte de agua en el precalentamiento del aire.
c) Estimar la humedad absoluta de salida del aire del secadero, suponiendo que sale saturado, y calcular el gasto de aire.
d) Calcular el consumo energtico de combustible de 40 MJ/kg de poder calorfico, usando un quemador para calentar el aire.

Datos:

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

>	su1:="Aire":su2:="H2O":su3:="CCl2F2":dat:=[T1=(70+273)*K_,m[p]='(150e3/86400)*kg_/s_',T3=(17+273)*K_,y[w,3]=0.45,T4=(70+273)*K_,y[w,4]=0.04,T0=(15+273)*K_,p0=100e3*Pa_,phi0=0.6,PC=40e6*J_/kg_,c[ps]=2000*J_/(kg_*K_)];

Esquema:

>

>

Eqs. const.:

>

pvR12:=get_pv_data(su3):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:R12dat:=get_gas_data(su3),get_liq_data(su3):

a) Estimar la potencia trmica que requiere la evaporación del agua del producto, y compararla con la requerida para calentar el producto seco suponiendo que su capacidad trmica es de 2000 J/(kg×K).

>

m[p]=evalf(subs(dat,m[p]),3);m[ps]:=m[p]*(1-y[w,3]);m_[ps]:=subs(dat,m[ps]):'m[ps]'=evalf(%,2);m[w3]:=m[p]*y[w,3];m_[w3]:=subs(dat,m[w3]):'m[w3]'=evalf(%,2);eq4:='m[w4]/(m[ps]+m[w4])=y[w,4]';m_[w4]:=solve(eq4,m[w4]);m__[w4]:=subs(dat,m_[w4]):'m[w4]'=evalf(%,2);m[w]:='m[w3]-m[w4]';m_[w]:=m_[w3]-m__[w4]:'m[w]'=evalf(%,2);Q[eva]:=m[w]*h[lv];Q_[eva]:=subs(dat,Wdat,dat,m_[w]*h[lv0]):'Q[eva]'=evalf(%,2);Q[cal]:='m[ps]*c[ps]*(T4-T3)';Q_[cal]:=subs(dat,Q[cal]):'Q[cal]'=evalf(%,3);

i.e. en primera aproximacin basta considerar la evaporacin del agua.

.b) Establecer el balance msico de H2O y el balance energtico del secadero despreciando el posible aporte de agua en el precalentamiento del aire.

>

eqBMw:=m[a]*w1+m[p]*y[w,3]=m[a]*w2+m[p]*y[w,4];w1_:=evalf(subs(dat,w(phi0,T0,p0))):'w1'=evalf(%,2);eqBMw:=m[a]*(w2-w1_)=m_[w]:evalf(%,2);eqBE:='m[a]*h1+m[p]*y[w,3]*hw3+m[ps]*hps3=m[a]*h2+m[p]*y[w,4]*hw4+m[ps]*hps4';h1_:=subs(dat,Adat,Wdat,T=T1,dat,h(T,w1_)):'h1'=evalf(%,2);hw3_:=subs(dat,Wdat,c*(T3-T[tr])):'hw3'=evalf(%,2);hps3_:=subs(dat,c[ps]*(T3-T[tr])):'hps3'=evalf(%,2);h2_:=subs(dat,Adat,Wdat,T=T2,h(T,w2)):hw4_:=subs(dat,Wdat,c*(T4-T[tr])):'hw4'=evalf(%,2);hps4_:=subs(dat,c[ps]*(T4-T[tr])):'hps4'=evalf(%,2);eqBE_:=subs(h1=h1_,h2=h2_,hw3=hw3_,hps3=hps3_,hw4=hw4_,hps4=hps4_,dat,eqBE):'eqBE'=evalf(%,2);

c) Estimar la humedad absoluta de salida del aire del secadero, suponiendo que sale saturado, y calcular el gasto de aire.

Despreciando el efecto del producto seco, se tratara de una tpica humidificacin adiabtica.

Usando el diagrama h-w se llega a que la w de salida es de 0,024. Analticamente:

>	eqBEha:='h1=h2(T2,w2)';w2_:=subs(dat,w(1,T2,p0)):'w2'=evalf(%,2);eqBEha:=subs(w2=w2_,h1_=h2_):T2_:=fsolve(subs(SI0,eqBEha),T2=250..350)*K_:'T2'=evalf(%,4);w2__:=evalf(subs(T2=T2_,dat,w2_)):'w2'=evalf(%,2);m_[a]:=subs(w2=w2__,dat,solve(eqBMw,m[a])):'m[a]'=evalf(%,2);

Nótese que con esa aproximación la eqBE no vale para calcular m[a].

Si se tiene en cuenta el calentamiento del producto seco, puede hacerse iterando en el diagrama h-w a lo largo de phi=1. Analticamente:

>	T2_:=fsolve(subs(w2=w2_,SI0,solve(eqBMw,m[a])=solve(eqBE_,m[a])),T2=100..400)*K_:'T2'=evalf(%,4);m_[a]:=evalf(subs(w2=w2_,T2=T2_,solve(eqBMw,m[a]))):'m[a]'=evalf(%,3);

d) Consumo energtico de combustible de 40 MJ/kg de poder calorfico, usando un quemador para calentar el aire.

>	Q[a]:=m[a]*c[pa]*(T1-T0);Q_[a]:=subs(m[a]=m_[a],Adat,dat,Q[a]):'Q[a]'=evalf(%,2);m[F]:='Q[a]/PC';m_[F]:=subs(m[a]=m_[a],Adat,dat,m[F]):'m[F]'=evalf(%,2);m[F,aprox]:=Q[w]/PC;m_[F,aprox]:=subs(Q[w]=Q_[eva],dat,m[F,aprox]):'m[F,aprox]'=evalf(%,2);

>

La aproximación era pensar que sólo se gasta energa en evaporar el agua.

En cuanto al posible aporte de agua debido al quemador, depender de su contenido de hidrógeno; el mximo tendra lugar quemando gas natural, que generara mF*9/4=0,14 kg/s de agua que ya empieza a ser significativo frente a los 0,74 kg/s de agua a evacuar y necesitara algo ms de aire. En realidad se quema carbón (turba) que apenas contiene hidrógeno, o se usa una caldera para calentar agua y que sta caliente el aire.