> restart:#"m16_p23"

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

[Ateo, Mf, PCI, PCS, eqEQ, eqMIX, eq_fit, get_hgs_data, hgs_r25, nulist, seqEBE]

En un combustor adiabático se quema metano con un 110% de aire. Se pide:
a) Composición y temperatura de salida suponiendo combustión completa.
b) Estimar la fracción de monóxido de carbono y de monóxido de nitrógeno a la salida suponiendo que se alcanza el equilibrio químico.
c) Hacer sendos esquemas del combustor, incluyendo un perfil típico de llama, suponiendo que la alimentación fuese premezclada o separada.

Datos:

> fuel:=CH4;dat:=[lambda=1.1];

CH4
[lambda = 1.1]

Eqs. const.:

> dat:=op(dat),op(subs(g=g0,[Const])),SI2,SI1:

a) Composición y temperatura de salida suponiendo combustión completa

> eq:=eqMIX(a*fuel+a*A*(c21*O2+c79*N2)=[2,3,4,6]);evalf(Ateo(fuel),2);A0_:=rhs(Ateo(fuel)):A=A[0]*lambda;eqDat:=A=subs(dat,A0_*lambda):evalf(%,3);eqNX;eqBC;eqBH;eqBO;eqBN;sol1_:=solve(subs(eqDat,dat,{eqNX,eqBC,eqBH,eqBO,eqBN}),{a,x[O2],x[N2],x[CO2],x[H2O]}):evalf(%,2);i:='i':eq15_5;PCI_:=subs(sol1_,PCI(eq)/a):'PCI'=evalf(%,3);i:='i':eq15_7_2;Ta_:=rhs(subs(sol1_,cpComp_,dat,eq15_7_3));

`+`(`*`(a, `*`(CH4)), `*`(a, `*`(A, `*`(`+`(`*`(c21, `*`(O2)), `*`(c79, `*`(N2))))))) = `+`(`*`(x[O2], `*`(O2)), `*`(x[N2], `*`(N2)), `*`(x[CO2], `*`(CO2)), `*`(x[H2O], `*`(H2O)))
A[0] = 9.6
A = `*`(A[0], `*`(lambda))
A = 10.5
1 = `+`(x[O2], x[N2], x[CO2], x[H2O])
0 = `+`(x[CO2], `-`(a))
0 = `+`(`*`(2, `*`(x[H2O])), `-`(`*`(4, `*`(a))))
0 = `+`(`-`(`*`(2, `*`(a, `*`(A, `*`(c21))))), `*`(2, `*`(x[O2])), `*`(2, `*`(x[CO2])), x[H2O])
0 = `+`(`-`(`*`(2, `*`(a, `*`(A, `*`(c79))))), `*`(2, `*`(x[N2])))
{a = 0.87e-1, x[CO2] = 0.87e-1, x[H2O] = .17, x[N2] = .72, x[O2] = 0.18e-1}
PC = `+`(`-`(Sum(`*`(nu[i], `*`(h[i])), i = 1 .. C)))
PCI = `+`(`/`(`*`(0.802e6, `*`(J_)), `*`(mol_)))
Ta = `+`(T25, `/`(`*`(a, `*`(PCI)), `*`(Sum(`*`(x[Com[i]], `*`(c[p, i])), i = 1 .. CP))))
`+`(`*`(2136., `*`(K_)))

b) Estimar la fracción de monóxido de carbono y de monóxido de nitrógeno a la salida suponiendo que se alcanza el equilibrio químico.

Serán muy pequeñas y la temperatura no disminuirá apreciablemente.

> subs(Comp='Comp',eq9_20);eq9_25;eqEQCO_:=evalf(subs(p=p0,T=Ta_,eqEQ(CO+(1/2)*O2=CO2)));xCO_:=solve(subs(sol1_,eqEQCO_),x[CO]):'x[CO]'=evalf(%,2);eqEQNO_:=evalf(subs(p=p0,T=Ta_,eqEQ(N2+O2=2*NO)));xNOeq_:=solve(subs(sol1_,eqEQNO_),x[NO])[1]:'x[NO]'=evalf(%,2);

Product(`^`(x[Comp[i]], nu[i]), i = 1 .. C) = `*`(`^`(`/`(`*`(p0), `*`(p)), Sum(nu[i], i = 1 .. C)), `*`(K(T)))
lnK = `/`(`*`(`+`(`-`(grR), `*`(hrR, `*`(`+`(1, `-`(`/`(`*`(T25), `*`(T)))))))), `*`(R[u], `*`(T25)))
`/`(`*`(x[CO2]), `*`(`^`(x[O2], `/`(1, 2)), `*`(x[CO]))) = 249.6
x[CO] = 0.26e-2
`/`(`*`(`^`(x[NO], 2)), `*`(x[O2], `*`(x[N2]))) = 0.7595e-3
x[NO] = 0.31e-2

c) Hacer sendos esquemas del combustor, incluyendo un perfil típico de llama, suponiendo que la alimentación fuese premezclada o separada.

No se han representado variantes típicas, cómo las bocas múltiples en los combustores de premezcla (para aumentar la estabilidad de la llama), ni la recirculación de aire en los combustores de gases no premezclados (para reducir el tamaño de la llama).

Para una potencia térmica dada (i.e. para caudales dados), el área de inyección vendría fijada por el rango de velocidades estables para llama de premezcla, en el primer caso, y por la longitud deseada de llama en el segundo combustor, que sería de mayor tamaño en cualquier caso.

Image

Fig. 1. Esquemas del combustor de premezcla y de no-premezcla.

>