Diferencia entre revisiones de «Reacciones de autocatalisis Grupo 9A»
(→Introducción) |
(→Ecuación diferencial) |
||
| Línea 88: | Línea 88: | ||
== Ecuación diferencial == | == Ecuación diferencial == | ||
| − | + | Procedemos a resolver el PVI mediante metodos numericos estudiados en la asignatura interpretando los resultados relacionados con el proceso quimico: | |
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
=== Método de Euler === | === Método de Euler === | ||
Revisión del 19:04 4 mar 2015
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Reacciones con autocatálisis. Grupo A2 |
| Asignatura | Ecuaciones Diferenciales |
| Curso | Curso 2014-15 |
| Autores | David Carmona Rodriguez,Alejandro Muñoz Cotter, Daniel Alonso Palop, Luis Bermeosolo Echeverria |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Contenido
1 Introducción y planteamiento
La autocatálisis es el proceso mediante el cual un compuesto químico induce y controla una reacción química sobre sí mismo. Los compuestos autocatalíticos no son catalizadores en sentido estricto ya que su estructura química resulta alterada durante el proceso.
Consideramos una solución bien mezclada a temperatura y volumen constantes. En esta solución tiene lugar una reacción química en la que en el momento inicial se encuentran dos reactivos A y B. A medida que avanza el tiempo se forma el producto 2B, teniendo en cuenta que la presencia de B hace de efecto catalítico en la reacción y satisfaciendo 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.
Reacción bimolecular: [math] A + B \rightarrow _{k1} 2B [/math]
Para tiempo [math](t=0)[/math] nombramos e identificamos las variables:
A: concentracion inicial [math]x_0[/math] (mol/l)
B:concentracion inicial [math]y_0[/math] (mol/l)
Sucede la reaccion[math] \rightarrow [/math] El tiempo comienza [math](t\gt0)[/math]
A: concentracion en funcion de t. [math]x(t)[/math] (mol /l)
B:concentracion en función de t. [math]y(t)[/math] (mol/l)
Como el volumen se mantiene constante [math](V=cte)[/math]
[math]V=volumen[/math]
[math]M_A(t)[/math] = masa del compuesto A
[math]M_B(t)[/math] = masa del compuesto B
Según la ley de concentración de la masa: [math]M_A(t)+M_B(t)=k[/math] siendo [math]k=''cte''[/math] en el proceso, si dividimos por V:
[math]\frac{M_A(t)}{V} + \frac{M_B(t)}{V} =\frac{k}{V}[/math] (renombramos [math]k^*=\frac{K}{V}[/math])
Sustituimos por nuestros términos:
[math]x(t)+y(t) = k^*[/math]
Derivando [math]x'(t)=+y'(t)=0[/math] ec.(1)
Segun la ley de acción de masas:
[math]\mbox{Velocidad de reacción = (cte)·(Cantidad de reactivo A)·(cantidad de reactivo de B)}[/math]
[math]y'(t)=k1*x(t)*y(t)[/math] ec.(2)
si integramos la ec.(1) [math]x(t)+y(t)=k^*[/math] (con [math]k^*=\frac{k}{V}[/math]) despejamos [math]x(t)=k^*-y(t)[/math]
y sustituyendo en ec.(2) ya tenemos planteado el P.V.I con las condiciones iniciales dadas en el enunciado :
[math] \left \{
\begin{matrix}
y'(t) = k_{1}*(cte - y(t))*y(t) \\
y(0)=0.01
\end{matrix}
\right . [/math]
Ahora procedemos a comprobar si tiene solucion y ademas es unica mediante la aplicacion del teorema de la existencia y unicidad:
Existe solución para el PVI planteado si existe una "Bola" de radio [math] r [/math] alrededor del punto de estudio [math] (t_{0},y_{0}) [/math] en la que la función sea continua. Esta solución será además única cuando también se cumpla que [math] \displaystyle{ \delta f \over \delta y} [/math] sea también contínua en [math] D \cap B((t_{0},y_{0}),r) [/math].
Observamos que nuestra funcion: [math] f(t,y(t))=f= k_{1}*(cte - y(t))*y(t) [/math] Y: [math] \displaystyle{ \delta f \over \delta y} = k_{1}(cte - 2y) [/math]
Ambas funciones son de clase [math]C^\infty[/math] podemos afirmar entonces, que existe solución única.
2 Ecuación diferencial
Procedemos a resolver el PVI mediante metodos numericos estudiados en la asignatura interpretando los resultados relacionados con el proceso quimico:
2.1 Método de Euler
Resolvemos la ecuacion mediante el programa de EULER
2.2 Método del Trapecio
Resolvemos la ecuacion mediante el metodo del trapecio.
2.3 Método de Runge-Kutta
Resolvemos mediante Rounge Kutta
