>	restart:#"m15_p13"

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

Calcular la temperatura de combustión adiabática del octano con 400% de aire teórico.

Datos:

>	su1:="Aire":su2:="H2O":fuel:=C8H18:dat:=[A=4*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) Temperatura de combustión adiabátic

>

eqEST:=eq_fit(fuel+a*O2=b*CO2+c*H2O);eq15_2;eqA0:=Ateo(fuel);A_:=subs(dat,A[0]=rhs(eqA0),A);PCI_:=PCI(eqEST);eqMX:=eqMIX(a*fuel+b*(c21*O2+c79*N2)=[2,3,4,6]);eqDat:=b/a=A_;sol1:=subs(dat,solve({eqNX,eqBC,eqBH,eqBO,eqBN,eqDat},{a,b,x[Comp[2]],x[Comp[3]],x[Comp[4]],x[Comp[6]]}));i:='i':subs(sol1,cpComp,dat,eqTa);

Podemos ver cómo varía la Ta con el aire.

>	n:=9:for j from 0 to n do eqDat:=b/a=(1+3*(j/n))*rhs(eqA0);sol1:=solve(subs(dat,{eqNX,eqBC,eqBH,eqBO,eqBN,eqDat}),{a,b,x[N2],x[O2],x[CO2],x[H2O]}):i:='i':Ta||j:=subs(sol1,cpComp,dat,rhs(eqTa)):od:plot([[[1,0],[1,4000]],[seq([(1+3*(j/n)),Ta||j/K_],j=0..n)]],lambda=0..4,Ta=0..3000,color=black);

i.e. la temperatura adiabática disminuye desde el máximo (con el aire teórico) de 2200 K y es de 880 K para lambda=A/A0=4.

>