Reacciones con autocatálisis, Grupo C26
1 INTRODUCCION
Este trabajo estudia las concentraciones de reactivos y productos a lo largo del tiempo en una reaccíon química bimolecular irreversible, con un volumen y una temperatura constantes y con la particularidad de que uno de los componentes (B) hace de efecto catalítico y actuará tanto de reactivo como de producto,es decir, partiendo de un contenido de A y B iniciales, reaccionaran para producir B que a su vez participara de nuevo como reactivo, por lo que es de la siguiente forma:
A +B →2B
Supondremos que en nuestra reacción se satisface la Ley de Accion de masas, que establece que la velocidad de reaccion es proporcional al producto de la concentracion de los reactivos, y la ley de de conservacion de masa ,por la que definimos que la masa de reactivos mas productos es constante en el tiempo Para la resolución del problema se usan los metodos de cálculo de Euler, del Trapecio , de Runge-Kutta
2 INTERPRETACION DEL PROBLEMA Y EXPRESION NUMÉRICA
comenzamos nombrando las variables y expresando el problema numéricamente: contenido de A= x(t) [mol/l] contenido de B= y(t) [mol/l]
Ley de Acción de Masas: por definición es:
y'(t)=K*x(t)*y(t) (I)
y'(t)= velocidad de creación del producto [mol/(l*s)] k = constante de proporcionalidad mol/s
Ley de conservación de masas
masa reactivos + masa productos = cte
x(t)+y(t) =(C) (II)
derivamos esta expresión y se obtiene:
x'(t)+y'(t)=0 x'(t)=-y'(t) (III)
La velocidad de desaparición de los reactivos es igual a la de aparicion de los productos(que es un concepto muy intuitivo) Por lo tanto la formula (III) se explica a partir derivar respecto al tiempo de La ley de conservación de masas. ahora sustitumos el contenido de A , x(t) de (II) en (I)
y'(t)=k*y(t)*(C-y(t)) → y'(t)=k*C*y(t)-K*(y(t)) (IV)
además conocemos las concentraciones iniciales(en t=0) y el valor de la k. A(0)=1mol/l B(0)=0,01mol/l k=mol/s
por lo que sustituyendo estos valores en (II) para deducir la (C) y en (IV) planteamos el problema de valor inicial
{y'(t)= 1,01*y(t)-y²(t)
{ y(0)=0,01 (t,y) ∈ [0,10]x R'
| |