Diferencia entre revisiones de «Reacciones complejas GRUPO 1A»

De MateWiki
Saltar a: navegación, buscar
Línea 1: Línea 1:
{{ TrabajoED | Reacciones complejas grupo 1A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Pablo Medina Higueras 1326;Jesús Caballero Pozo 1392;Jaime Delage Ramírez 1425 }}
+
{{ TrabajoED | Reacciones complejas Grupo 1A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Pablo Medina Higueras  
 +
Jesús Caballero Pozo  
  
 +
Jaime Delage Ramírez  }}
  
==Enunciado==
+
 
 +
==Introducción==
  
 
Se considera una reacción química irreversible en una solución bien mezclada. Supondremos que la reacción ocurre para un volumen y temperatura constantes. Al inicio se encuentran dos reactivos A y B que van formando un producto C en lo que se conoce como una reacción bimolecular, es decir, una molécula de A y una de B producen una de C,
 
Se considera una reacción química irreversible en una solución bien mezclada. Supondremos que la reacción ocurre para un volumen y temperatura constantes. Al inicio se encuentran dos reactivos A y B que van formando un producto C en lo que se conoce como una reacción bimolecular, es decir, una molécula de A y una de B producen una de C,
 
                                                   <big>'''A + B → C'''</big>
 
                                                   <big>'''A + B → C'''</big>
Supondremos que se satisface la ley de acción de masas que establece que la velocidad de reacción es proporcional al producto de las concentraciones de los reactivos. Se pide:
+
Supondremos que se satisface la ley de acción de masas que establece que la velocidad de reacción es proporcional al producto de las concentraciones de los reactivos.  
 
+
1. Comprobar que la concentración del producto C a lo largo del tiempo puede obtenerse a partir de la ecuación diferencial
+
+
                                  <math>y'(t)=k_{1}</math>(<math>a_{0}</math> −<math>y(t)</math>)(<math>b_{0}</math> −<math>y(t)</math>), <math>t > 0</math>,
+
 
+
para algunas constantes <math>a_{0}</math>, <math>b_{0}</math> y <math>k_{1}</math>. Interpretar dichas constantes y enunciar un problema de
+
valor inicial adecuado para calular la concentración de C a lo largo del tiempo. ¿Tiene solución única?
+
 
+
2. ¿Cómo cambiaría la ecuación diferencial si el proceso fuera reversible?
+
 
+
3. Suponiendo <math>a_{0}</math> = 3 mol/l, <math>b_{0}</math> = 1 mol/l y <math>k_{1}</math> = 1 mol/s resolver el PVI por el método de Euler, eligiendo un paso h = 0.1, en los primeros 2 segundos. Dibujar en una misma gráfica las concentraciones de los dos reactivos y el producto. Interpretar las gráficas.
+
 
+
4. Resolver numéricamente el sistema para un tiempo suficientemente grande como para establecer una aproximación de los límites de las concentraciones cuando t → ∞. ¿Es lógico lo que se observa?
+
 
+
5. Resolver el PVI numéricamente usando el método del Trapecio y el método de Runge-Kutta de cuarto orden con el mismo paso de tiempo.
+
  
6. Supongamos ahora que se produce una reacción consecutiva de la forma
 
                                    <big>'''A + B →<math>k_{1}</math> C →<math>k_{2}</math> D'''</big>
 
  
Plantear ahora un sistema de dos ecuaciones diferenciales para las concentraciones de C y D, basado en la ley de acción de masas. Obtener el PVI asociado.
+
==Interpretación de las constantes y PVI para calcular la concentración de C==
  
7. Si <math>k_{2}</math> = 5, lo que supone que la segunda reacción es mucho más rápida que la primera, resolver numéricamente el problema en los primeros 10 segundos con los métodos de Euler y Runge-Kutta. Dibujar las gráficas de las concentraciones en una misma figura e interpretar los resultados. ¿Se puede usar el paso h = 0.3?
+
En primer lugar, interpretaremos las constantes y variables de la siguiente ecuación diferencial,esta ecuación diferencial será la de la ley de masas:
 +
                                          <big>'''y'(t) = k1(a0 − y(t))(b0 − y(t)), siendo t>0'''</big>  
 +
donde:
 +
  y'(t) = Representa la velocidad de la reacción química, es decir, la velocidad con que se produce C.
 +
  k1 = Constante de proporcionalidad.
 +
  a0 = Concentración inicial de A.
 +
  y(t) = Evolución de la concentración de  C a lo largo del tiempo.
 +
  b0 = Concentración inicial de B.
  
8. Si <math>k_{2}</math> = 1/5, lo que supone que la segunda reacción es mucho más lenta que la primera, resolver numéricamente el problema en los primeros 10 segundos con los métodos de Euler y Runge-Kutta. Dibujar las gráficas de las concentraciones en una misma figura e interpretar los resultados.
+
Es conveniente mencionar que la concentración de A respecto al tiempo será de a0 - y(t), y para B será de b0 - y(t), tal y como nos indica el enunciado, tomaremos el tiempo como  mayor que 0: t>0.
  
 +
Además, sabemos que para t = 0, la concentración de C será nula, por lo que: y(t = 0) = 0.
  
 +
Con todo lo deducido anteriormente, llegamos a la conclusión que estamos ante un Problema de Valor Inicial (PVI) o de Cauchy:
 +
                                        <big>'''y'(t) = k1(a0 − y(t))(b0 − y(t)), siendo t>0'''</big>
 +
                                                            <big>'''y(0) = 0'''</big>
  
 +
Ahora procederemos a ver si nuestro problema tiene solución única o no, a través del Teorema de existencia y Unicidad. Observamos que la función f(t,y) = y'(t) es continua en el intervalo (I = (0,∞) ∩ B(0,0),r>0), por lo que admite al menos una solución. Por otro lado, observamos la derivada parcial:
 +
                                                            <big>'''∂f/∂y=k(2y-a-b)'''</big>
 +
y vemos que no hay problemas de continuidad, porque la derivada parcial nos resulta un polinomio, así que podremos afirmar que tiene una única solución.
 +
 
  
  

Revisión del 20:43 27 abr 2016

Trabajo realizado por estudiantes
Título Reacciones complejas Grupo 1A
Asignatura Ecuaciones Diferenciales
Curso Curso 2015-16
Autores Pablo Medina Higueras

Jesús Caballero Pozo

Jaime Delage Ramírez

Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura


1 Introducción

Se considera una reacción química irreversible en una solución bien mezclada. Supondremos que la reacción ocurre para un volumen y temperatura constantes. Al inicio se encuentran dos reactivos A y B que van formando un producto C en lo que se conoce como una reacción bimolecular, es decir, una molécula de A y una de B producen una de C,

                                                  A + B → C

Supondremos que se satisface la ley de acción de masas que establece que la velocidad de reacción es proporcional al producto de las concentraciones de los reactivos.


2 Interpretación de las constantes y PVI para calcular la concentración de C

En primer lugar, interpretaremos las constantes y variables de la siguiente ecuación diferencial,esta ecuación diferencial será la de la ley de masas:

                                         y'(t) = k1(a0 − y(t))(b0 − y(t)), siendo t>0 

donde:

  y'(t) = Representa la velocidad de la reacción química, es decir, la velocidad con  que se produce C.
  k1 = Constante de proporcionalidad.  
  a0 = Concentración inicial de A.
  y(t) = Evolución de la concentración de  C a lo largo del tiempo.
  b0 = Concentración inicial de B.

Es conveniente mencionar que la concentración de A respecto al tiempo será de a0 - y(t), y para B será de b0 - y(t), tal y como nos indica el enunciado, tomaremos el tiempo como mayor que 0: t>0.

Además, sabemos que para t = 0, la concentración de C será nula, por lo que: y(t = 0) = 0.

Con todo lo deducido anteriormente, llegamos a la conclusión que estamos ante un Problema de Valor Inicial (PVI) o de Cauchy:

                                       y'(t) = k1(a0 − y(t))(b0 − y(t)), siendo t>0 
                                                            y(0) = 0 

Ahora procederemos a ver si nuestro problema tiene solución única o no, a través del Teorema de existencia y Unicidad. Observamos que la función f(t,y) = y'(t) es continua en el intervalo (I = (0,∞) ∩ B(0,0),r>0), por lo que admite al menos una solución. Por otro lado, observamos la derivada parcial:

                                                           ∂f/∂y=k(2y-a-b) 

y vemos que no hay problemas de continuidad, porque la derivada parcial nos resulta un polinomio, así que podremos afirmar que tiene una única solución.