Diferencia entre revisiones de «Reacciones con autocatálisis, Grupo C26»

Revisión del 20:24 6 mar 2015

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 →k2B

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'

UNICIDAD DE SOLUCIÓN

Aplicamos a continuación el Teorema de Existencia y Unicidad:

Sea f(t,y)=y′(t)⇒f(t,y)=1.01y(t)−y(t), y δf(t,y)/δx=1.01−2y(t).

y B[(t0,y0),r], r>0,en(t0,y0)=(0,0.01)

f(t,y) es continua en B∩D⇒∃ y(t)/y′(t)=1,01y(t)−y²(t)

y su derivada respecto de y →f'(t,y) es continua en la intersección con el dominio

Por lo que hemos demostrado que sí tiene solución y ésta es única.

3 RESOLUCION NUMÉRICA DEL PROBLEMA DE VALOR INICIAL

ahora procedemos a resolver el P.V.I dado anteriormente por loe métodos numéricos de Euler, del Trapecio y de Runge-kutta. Estudiamos el problema a lo largo de los 10 primeros segundos.

3.1 METODO DE EULER

%Método de Euler

%Datos iniciales
t0=0;
y0=0.01;
h=0.1;
t0=0;
tN=10;
 
%Cálculo de subintervalos equidistantes
N=(tN-t0)/h;
t=t0:h:tN;
 
%vector solución:
y=zeros(1,N+1);
y(1)=y0;
 
%aplicamos el metodo de Euler
for i=1:N;
    y(i+1)=y(i)+h*(1.01*y(i)-(y(i)).^2);
end;
%Tabla de Resultados
[t',y'];
x=1.01-y;
hold on ;
plot(t,y,'r')
plot(t,x,'b')
hold off;
legend('Sustancia B: y(t)','Sustancia A: x(t)','Location','best');
xlabel('Tiempo (s)')
ylabel('Concentración (mol/l)')
grid on

Con esto vemos que si bien al principio la veocidad de reaccion (o de creacion de producto)es lenta porque hay poco contenido en B y mucho en A, va progresando hasta su punto máximo cuando A Y B son iguales y luego vuelve a disminuir porque se va agotando A , y esto ocurre por la Ley de acción de Masas y'(x)=kx(t)y(t)

3.2 METODO DEL TRAPECIO ahora calculamos el P.V.I con el metodo del trapecio con el mismo paso del tiempo que anteriormente. %Método del Trapecio

%Datos iniciales
t0=0;
y0=0.01;
h=0.1;
t0=0;
tN=10;
 
%calculo de subintervalos equidistantes:
N=(tN-t0)/h;
t=t0:h:tN;
 
%vector solución:
y=zeros(1,N+1);
y(1)=y0;
%metodo del trapecio:
for i=1:N
y(i+1)=(((1-(h/2)*1.01)-sqrt(((1-(h/2)*1.01)^2)+4*(h/2)*((1+(h/2)*1.01)*y(i)-(h/2)*(y(i))^2)))/-h);
end;
%Sacamos la tabla de resultados
[t',y'];
x=1.01-y;
%Gráfico
hold on;
plot(t,y,'r');
plot(t,x,'b');
legend('Sustancia B: y(t)','Sustancia A: x(t)','Location','best');
xlabel('Tiempo (s)')
ylabel('Concentración (mol/l)')

grid on

3.3 METODO DE RUNGE KUTTA DE ORDEN CUATRO

El método de Runge Kutta es mas complicado por lo que hay que usar funciones ya en Matlab para que sea mas sencillo de expresarlo

%Método de Runge-Kutta
 
%Función f(t)=y'(t)
fy=inline('1.01*y-y^2','y');
 
%Datos del problema
y0=0.01;
t0=0;
tN=10;
h=0.1;
 
t=t0:h:tN;
y=zeros(1,length(t));
y(1)=y0;
 
for i=1:length(t)-1
    k1=fy(y(i));
    k2=fy(y(i)+k1/2*h);
    k3=fy(y(i)+k2/2*h);
    k4=fy(y(i)+k3*h);
    y(i+1)=y(i)+h/6*(k1+2*k2+2*k3+k4);
end
 
x=1.01-y;
 
hold on
plot(t,x,'b')
plot(t,y,'r')
legend('A: x(t)','B: y(t)')
xlabel('Time(s)')
ylabel('Concentración (mol/l)')
grid on

5. REACCIÓN DE LOTKA. INTERPRETACIÓN

Consideramos a partir de ahora la reacción consecutiva propuesta por Lotka en 1920:

A + X --> {k1} 2X
X + Y --> {k2} 2Y                      
Y --> {k3} B

Donde A, X, Y, D son sustancias distintas e interpretaremos como la concentración de dichas sustancias. La reacción consume A para producir B y la velocidad y mezcla de las reacciones están controladas por X e Y. Según la ley de la conservación de la masa: "En una reacción química ordinaria la masa permanece constante, es decir, la masa consumida de los reactivos es igual a la masa obtenida por los productos". En nuestro caso, al no tratarse de una reacción nuclear, podemos aplicar esta ley de la siguiente manera: la suma de las concentraciones de las sustancias implicadas debe de permanencer constante.

A + X + Y + B = cte

Sin embargo, debemos recurrir también a la ley de acción de masas: "En una reacción química reversible en equilibrio a una temperatura constante, una relación determinada de concentraciones de reactivos y productos, tienen un valor constante". En este punto, es crítico ver que las sustancias X e Y no solo se consumen, sino que las primeras dos reaccionas también las producen. Por ello mismo, restaremos la porción producida a la porción inicial o consumida.

A partir de esto podemos calcular el sistema de ecuaciones que corresponde a nuestro problema:

X' = k1*A*X - k2*X*Y
Y' = k2*X*Y - k3*Y                      
B' = k3*Y

Si derivamos la expresión obtenida a partir de la ley de conservación de la masa, obtenemos:

A' + X' + Y' + B' = 0

Expresión en función de X', Y' y B', las cuales ya tenemos, y podemos por tanto sustituir en nuestro problema:

A' + k1*A*X - k2*X*Y + k2*X*Y - k3*Y + k3*Y = 0
A' + k1*A*X = 0                      
A' = -k1*A*X

Donde k1 es la velocidad de reacción.

Por ello, en los próximos apartados deberemos resolver el siguiente problema de valor inicial con los datos que se proporcionan, que son:

k1 = k2 = 2*k3 = 0.1
Y' = 0.1*X*Y - 0.1/2*Y
B' = 0.1/2*Y                     
A' = -0.1*A*X

con valores iniciales:

A(0) = 5; X(0)=0.0005; Y(0)=0.00001; B(0)=0