>	restart:#"m15_p21"

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

Se quema acetileno isobáricamente con 120% de aire teórico. Se pide:

a) Aire teórico.

b) Poder calorífico.

c) Temperatura adiabática.

Datos:

>	su1:="Aire":su2:="H2O":fuel:=C2H2:dat:=[A=1.2*A[0]];

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.

>	Comp[1];eq15_2;eq1:=Ateo(fuel);

i.e. 11,9 mol de aire por mol de acetileno.

b) Poder calorífico.

Por ser con exceso de aire, será igual al teórico:

>	eq:=eq_fit(fuel+a*O2=b*CO2+c*H2O);PCI_:=PCI(eq);PCS_:=PCS(eq);PCSm_:=subs(dat,%/rhs(Mf(fuel)));

i.e., 1,3 MJ/mol (50 MJ/kg); nótese que son casi iguales porque el acetileno tiene poco hidrógeno y genera poca agua.

c) Temperatura adiabática.

>	i:='i':eq15_7_2;eqMIX(subs(dat,A[0]=rhs(eq1),a*fuel+a*A*(c21*O2+c79*N2))=[2,3,4,6]);sol1_:=solve(subs(dat,{eqNX,eqBC,eqBH,eqBO,eqBN}),{a,x[O2],x[CO2],x[H2O],x[N2]});subs(sol1_,cpComp_,dat,eqTa);

i.e. unos 2560 K.

>