539.sagews

# calculo con ecuacion 3 y 5
a = 1 # reactivo 1 (CH3OH) metanol
b = 1 # reactivo 2 (CH3COOCH2CH3) acetato de etilo
c = 1 # producto 1 (CH3CH2OH) etanol
d = 1 # producto 2 (CH3COOCH3) acetato de metilo
dn = c + d - a - b
jsmath( str(a) + "A \ + \ " + str(b) + "B \ <===> \ " + str(c) + "C \ + \ " + str(d) + "D")

\newcommand{\Bold}[1]{\mathbf{#1}}1A \ + \ 1B \ [removed] \ 1C \ + \ 1D

# entropias para cada especie
Sa = 239.9  #126.8 # entropia del reactivo 1
Sb = 362.75 #257.7 # entropia del reactivo 2
Sc = 281.6  #160.7 # entropia del producto 1
Sd = 324.4  #207.8 # entropia del producto 2
Srxn = Sc + Sd - Sa - Sb; jsmath("\Delta S \ = \ " + str(float(Srxn)) + "\ J/molK")

\newcommand{\Bold}[1]{\mathbf{#1}}\Delta S \ = \ 3.35\ J/molK

# entalpias para cada especie
Ha = -201.0 #-239.2 # entalpia del reactivo 1
Hb = -443.6 #-479.3 # entalpia del reactivo 2
Hc = -234.8 #-277.6 # entalpia del producto 1
Hd = -413.3 #-445.9 # entalpia del producto 2
Hrxn = Hc + Hd - Ha - Hb; jsmath("\Delta H \ = \ " + str(float(Hrxn)) + "\ kJ / mol")

\newcommand{\Bold}[1]{\mathbf{#1}}\Delta H \ = \ -3.5\ kJ / mol

# constantes y paquetes a utilizar
import fpformat
dec = 4 # cantidad de decimales a utilizar
T = 70 + 273.15 # temperatura alterable
Ts = 298.15 # temperatura estandar
R = 8.314472 # constante de gases ideales
Cr = 1 # concentracion de reactivo
dH = Hrxn * 1000 # entalpia de reaccion estandar
dG = dH - Ts * Srxn # energia libre de gibbs de reaccion estandar
jsmath("\Delta G \ = \ " + str(float(dG)) + "\ J / mol")

\newcommand{\Bold}[1]{\mathbf{#1}}\Delta G \ = \ -4498.8025\ J / mol

# calculo con ecuacion 1
Kp = e^-(dG/(R*Ts)); jsmath("K^{o}_{p}=" + str(fpformat.fix(float(Kp),dec)))

\newcommand{\Bold}[1]{\mathbf{#1}}K^{o}_{p}=6.1398

# calculo con ecuacion 4
K = Kp / (e^(1/Ts)/e^(1/T))^(dH/R); jsmath( "K_{p} =" + str(fpformat.fix(float(K),dec)))

\newcommand{\Bold}[1]{\mathbf{#1}}K_{p} =7.3887

# calculo con ecuacion 2
Ks = K / (R * Ts)^dn; jsmath( "K_{c}=" + str(fpformat.fix(float(Ks),dec)))

\newcommand{\Bold}[1]{\mathbf{#1}}K_{c}=7.3887

# calculo con ecuacion 6
C = solve(x^2 / (Cr - x)^2 == K, x, solution_dict = True)

\newcommand{\Bold}[1]{\mathbf{#1}}Posibles \ concentraciones \ de \ productos:

\newcommand{\Bold}[1]{\mathbf{#1}}x_{1} =0.7311\ mol / L

\newcommand{\Bold}[1]{\mathbf{#1}}x_{2} =1.5820\ mol / L

\newcommand{\Bold}[1]{\mathbf{#1}}Porcentaje \ de \ Rendmiento \ (\%R) \ = \ 73.1053\%