>	restart:#"m16_p21"

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

Representar gráficamente la hform, gform y hreac de los hidrocarburos saturados en función del número de carbonos, comentando la tendencia.

Datos:

Incluímos el H2 como si fuera el C=0, pero no es verdad.

Los consideramos todos en estado gaseoso, i.e. a presiones bajas.

>	dat:=op(subs(g=g0,[Const])),SI2,SI1:F:=["H2(g)","CH4(g)","C2H6(g)","C3H8(g)","C4H10(g)","C5H12(g)","C6H14(g)","C7H16(g)","C8H18(g)"];hgs0:=[get_hgs_data(F[1])];for i from 1 to 8 do cat(hgs,i):=[get_hgs_data(F[i+1])];od;

a) Representar gráficamente la hform.

>	plot(subs(SI0,[[0,hgs0[1]],[1,hgs1[1]],[2,hgs2[1]],[3,hgs3[1]],[4,hgs4[1]],[5,hgs5[1]],[6,hgs6[1]],[7,hgs7[1]],[8,hgs8[1]]]));

i.e. es casi lineal (quitando el H2).

b) Representar gráficamente la gform.

>	plot(subs(SI0,[[0,hgs0[2]],[1,hgs1[2]],[2,hgs2[2]],[3,hgs3[2]],[4,hgs4[2]],[5,hgs5[2]],[6,hgs6[2]],[7,hgs7[2]],[8,hgs8[2]]]));

Si no fuese por el heptano sería casi lineal (debe estar mal ese dato). Efectivamente, estaba mal en la fuente de datos que se estaba usando, y se corrigió. Para eso valen este tipo de ejercicios, para detectar anormalidades.

c) Representar gráficamente la hcomb.

>

eq:=eq_fit(H2+b*O2=c*CO2+d*H2O);hgsr0:=[hgs_r25(eq)];eq:=eq_fit(CH4+b*O2=c*CO2+d*H2O);hgsr1:=[hgs_r25(eq)];eq:=eq_fit(C2H4+b*O2=c*CO2+d*H2O);hgsr2:=[hgs_r25(eq)];eq:=eq_fit(C3H8+b*O2=c*CO2+d*H2O);hgsr3:=[hgs_r25(eq)];eq:=eq_fit(C4H10+b*O2=c*CO2+d*H2O);hgsr4:=[hgs_r25(eq)];eq:=eq_fit(C5H12+b*O2=c*CO2+d*H2O);hgsr5:=[hgs_r25(eq)];eq:=eq_fit(C6H14+b*O2=c*CO2+d*H2O);hgsr6:=[hgs_r25(eq)];eq:=eq_fit(C7H16+b*O2=c*CO2+d*H2O);hgsr7:=[hgs_r25(eq)];eq:=eq_fit(C8H18+b*O2=c*CO2+d*H2O);hgsr8:=[hgs_r25(eq)];plot(subs(SI0,[[0,hgsr0[1]],[1,hgsr1[1]],[2,hgsr2[1]],[3,hgsr3[1]],[4,hgsr4[1]],[5,hgsr5[1]],[6,hgsr6[1]],[7,hgsr7[1]],[8,hgsr8[1]]]));

La razón de que la hr sea tan lineal con n es que la hr de CnH2n+2 es igual a n*hCO2+(n+1)*hH2O-hfuel y la hfuel apenas cuenta si n es grande.

>