> restart:#"m09_p23"

La ecuación de Gibbs-Duhem enseña que en todo proceso a T y p constantes Snidmi=0. Se pide:
a) Poner la ecuación anterior en términos del grado de avance de la reacción.
b) Desarrollar el resultado anterior hasta verificar el aserto del enunciado.

Datos:

> read`../therm_chem.m`:with(therm_chem):with(therm_proc):i:='i':eqGD:=0=S*dT-V*dp+Sum(n[o]*dmu[i],i=1..C);

`:=`(eqGD, 0 = `+`(`*`(S, `*`(dT)), `-`(`*`(V, `*`(dp))), Sum(`*`(n[o], `*`(dmu[i])), i = 1 .. C)))

Esquema:

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

Eqs. const.:

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

a) Poner la ecuación anterior en términos del grado de avance de la reacción.

> i:='i':n[i]:=n[i,0]+xi*nu[i];nt:='sum(n[j,0]+xi*nu[j],j=1..C)';x[i]:='n[i]/nt';mu[i]:='mu[i,0]+R[u]*T*ln(x[i])';eqGD_Tp:='sum(n[i]*Diff(mu[i],xi),i=1..C)';eqGD_Tp_:=expand(value(%));'eqGD_Tp_'=expand(subs(C=9,%));

`:=`(n[i], `+`(n[i, 0], `*`(xi, `*`(nu[i]))))

`:=`(nt, sum(`+`(n[j, 0], `*`(xi, `*`(nu[j]))), j = 1 .. C))

`:=`(x[i], `/`(`*`(n[i]), `*`(nt)))

`:=`(mu[i], `+`(mu[i, 0], `*`(R[u], `*`(T, `*`(ln(x[i]))))))

`:=`(eqGD_Tp, sum(`*`(n[i], `*`(Diff(mu[i], xi))), i = 1 .. C))
`:=`(eqGD_Tp_, `+`(`-`(`/`(`*`(R[u], `*`(T, `*`(sum(`+`(`-`(`*`(nu[i], `*`(sum(`+`(n[j, 0], `*`(xi, `*`(nu[j]))), j = 1 .. C)))), `*`(sum(nu[j], j = 1 .. C), `*`(n[i, 0])), `*`(sum(nu[j], j = 1 .. C),...
eqGD_Tp_ = 0

>