>	restart:#"m15_p42"

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

Una caldera ha de suministrar 10 MW a un sistema de calefacción por agua caliente, quemando gas natural con un 10% de exceso de aire, saliendo los gases a 50 ºC por encima de su temperatura de rocío. Se pide:
a) Composición esperada de los gases de escape.
b) Temperatura adiabática de combustión y temperatura de rocío.
c) Caudales necesarios de gases de alimentación.

Datos:

>	su1:="Aire":su2:="H2O":fuel:=CH4:dat:=[Q=10e6*W_,e=0.1,DTs=50*K_];

Eqs. const.:

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

a) Composición esperada de los gases de escape.

>	eq15_4;eq1_:=Ateo(fuel);A_:=solve(subs(A[0]=rhs(eq1_),dat,eq15_4),A);eq:=eqMIX(a*fuel+a*A_*(c21*O2+c79*N2)=[2,3,4,6]);sol1_:=solve(subs(A[0]=rhs(eq1_),dat,{eqNX,eqBC,eqBH,eqBO,eqBN}),{a,x[Comp[2]],x[Comp[3]],x[Comp[4]],x[Comp[6]]});

b) Temperatura adiabática de combustión y temperatura de rocío.

>	PCI_:=subs(sol1_,PCI(eq)/a);eq15_7_2;eqTa_:=subs(b=a*A_,sol1_,cpComp_,c[p,fuel]=0,c[p,Air]=c[p,N2],Tef_=T25,Tea_=Tea,dat,eqTa);eq8_2;Trocio_:=solve(subs(dat,eval(subs(x[v,sat]=x[H2O],sol1_,p[v]=pv,p=p0,dat,eq8_2))),T);Trocio=TKC(Trocio_);

Ta=2140 K y Trocío=330 K.

c) Caudales necesarios de gases de alimentación.

>	n:='n':eqBE:=Q=n[fuel]*PCI;n_[fuel]:=subs(sol1_,dat,Q/PCI_):n[fuel]=evalf(%);m_[fuel]:=n_[fuel]*rhs(Mf(fuel)):m[fuel]=evalf(%);n:='n':eqA:=n[a]=n[fuel]*A;eqA:=n[a]=n_[fuel]*A_;eqAm:=m[a]=subs(Adat,n_[fuel]*A_*M);

>