> restart:#"m15_p03"

Determinar el aire teórico necesario y el poder calorífico en la combustión del octano:

a) En base molar.

b) En base másica.

Datos:

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

> su1:="Aire":su2:="H2O":su3:="C8H18":dat:=[]:

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

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:Fdat:=get_liq_data(su3):fuel:=convert(su3,symbol):

Estequiometría

> eq1:=C8H18+(25/2)*O2=8*CO2+9*H2O;

`+`(C8H18, `*`(`/`(25, 2), `*`(O2))) = `+`(`*`(8, `*`(CO2)), `*`(9, `*`(H2O)))

a) En base molar.

> eq15_2;A0_:=subs(dat,rhs(Ateo(fuel)));i:='i':eq15_5;PCI_:=PCI(eq1);PCS_:=PCS(eq1);

A[0] = `/`(`*`(`+`(u, `*`(`/`(1, 4), `*`(v)), `-`(`*`(`/`(1, 2), `*`(w))), y)), `*`(c21))
59.5
PC = `+`(`-`(Sum(`*`(nu[i], `*`(h[i])), i = 1 .. C)))
`+`(`/`(`*`(0.507e7, `*`(J_)), `*`(mol_)))
`+`(`/`(`*`(0.547e7, `*`(J_)), `*`(mol_)))

b) En base másica.

> Mf_:=rhs(Mf(fuel));Ma_:=subs(Adat,M);eq3:=A0m=A[0]*Ma/Mf;A0m_:=A0_*Ma_/Mf_;PCIm_:=PCI_/Mf_;PCSm_:=PCS_/Mf_;

`+`(`/`(`*`(.114, `*`(kg_)), `*`(mol_)))
`+`(`/`(`*`(0.29e-1, `*`(kg_)), `*`(mol_)))
A0m = `/`(`*`(A[0], `*`(Ma)), `*`(Mf))
15.2
`+`(`/`(`*`(0.445e8, `*`(J_)), `*`(kg_)))
`+`(`/`(`*`(0.480e8, `*`(J_)), `*`(kg_)))

i.e., el n-octano tiene un poder calorífica superior de 5,5 MJ/mol=48 MJ/kg; este último valor es representativo de todos los hidrocarburos, mientras que el valor molar varía mucho de uno a otro.

>