>	restart:#"m15_p18"

>	read`../therm_chem.m`:with(therm_chem);with(therm_proc):

Se quema butano con un 80% de exceso aire ambiente que está a 20 °C y tiene una humedad del 75%. Se pide:

a) Aire teórico suponiendo aire seco.

b) Temperatura de rocío de los gases de escape habiendo supuesto aire seco.

c) Aire teórico teniendo en cuenta la humedad.

d) Temperatura de rocío real de los gases de escape.

e) Poder calorífico inferior.

f) Temperatura de combustión adiabática

Datos:

>	su1:="Aire":su2:="H2O":fuel:=C4H10:dat:=[e=0.8,T0=(20+273)*K_,phi0=0.75];

Eqs. const.:

>	Adat:=get_gas_data(su1):Wdat:=get_gas_data(su2),get_liq_data(su2):get_pv_data(su2):dat:=op(dat),op(subs(g=g0,[Const])),Adat,SI2,SI1:

a) Aire teórico suponiendo aire seco.

>	eq1:=eq15_2;eqA0:=Ateo(fuel);

b) Temperatura de rocío de los gases de escape habiendo supuesto aire seco.

>

eq:=eqMIX(a*fuel+b*(c21*O2+c79*N2)=[2,3,4,6]);eqDat_:=b/a=subs(A[0]=rhs(eqA0),dat,(1+e)*A[0]):evalf(%,2);sol1_:=solve(subs(dat,dat,{eqNX,eqBC,eqBH,eqBO,eqBN,eqDat_}),{a,b,x[N2],x[O2],x[CO2],x[H2O]});eq8_2;Trocio_:=evalf(subs(dat,solve(eval(subs(x[v,sat]=x[H2O],p[v]=pv,sol1_,p=p0,dat,eq8_2)),T)));Trocio=TKC(Trocio_);

c) Aire teórico teniendo en cuenta la humedad

>	x[v,0]:=solve(subs(eq8_2,eq8_3),x[v]);xv0_:=evalf(subs(p=p0,T=T0,phi=phi0,dat,eval(subs(p[v]=pv,%)))):'x[v,0]'=evalf(%);A0h:=A0/(1-xv0);A0h:=rhs(eqA0)/(1-xv0_);

d) Temperatura de rocío real de los gases de escape.

>	eq:=eqMIX(a*fuel+b*(c21*O2+c79*N2+xv0_*H2O)=[2,3,4,6]);eqDat_:=b/a=subs(dat,(1+e)*A0h);sol1_:=solve(subs(dat,dat,{eqNX,eqBC,eqBH,eqBO,eqBN,eqDat_}),{a,b,x[N2],x[O2],x[CO2],x[H2O]});eq8_2;Trocio_:=evalf(subs(dat,solve(eval(subs(x[v,sat]=x[H2O],p[v]=pv,sol1_,p=p0,dat,eq8_2)),T)));Trocio=TKC(Trocio_);

e) Poder calorífico inferior

>	PCI_:=subs(sol1_,PCI(eq)/a);

f) Temperatura de combustión adiabática

>	i:='i':eq15_7_2;subs(cpComp,sol1_,dat,eqTa);

i.e. unos 1420 K.

>