> restart:#"m09_p14"

Calcular la fuerza electromotriz (f.e.m.) máxima de una pila electroquímica de H2 y O2 a 25 °C.

Datos:

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

> dat:=[T=(273+25)*K_,F=96500*J_/(mol_*V_)];

`:=`(dat, [T = `+`(`*`(298, `*`(K_))), F = `+`(`/`(`*`(96500, `*`(J_)), `*`(mol_, `*`(V_))))])

Esquema:

> `:=`(Sistemas, [VC])

Eqs. const.:

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

a) fem

> eq:=H2O=H2+(1/2)*O2;eqFEM:=epsilon=-gr/(z*F);eqEQ(eq):z:=2;hr_:=hgs_r25(eq)[1]:'hr'=evalf(%,2);gr_:=hgs_r25(eq)[2]:'gr'=evalf(%,2);subs(gr=gr_,dat,dat,eqFEM):evalf(%,2);

`:=`(eq, H2O = `+`(H2, `*`(`/`(1, 2), `*`(O2))))

`:=`(eqFEM, epsilon = `+`(`-`(`/`(`*`(gr), `*`(z, `*`(F))))))

`:=`(z, 2)

hr = `+`(`/`(`*`(0.29e6, `*`(J_)), `*`(mol_)))
gr = `+`(`/`(`*`(0.24e6, `*`(J_)), `*`(mol_)))

epsilon = `+`(`-`(`*`(1.2, `*`(V_))))

>