INGENIERÍA DE LAS REACCIONES
QUÍMICAS II
Introducción a
Polymath
POLYMATH UNA HERRAMIENTA
PARA ANÁLISIS NUMÉRICO
• Se utiliza Polymath para…
Interpolación Ecuaciones Diferenciales
Derivación Ecuaciones no Lineales
Integración Ecuaciones Lineales y Sistemas
• ¿Que otros programas se pueden utilizar?
MATLAB, OCTAVE, EXCEL, R, etc…
• ¿Porque usar Polymath?
Extremadamente sencillo de usar.
POLYMATH
SISTEMAS DE ECUACIONES
LINEALES
INTERPOLACIÓN Y REGRESIÓN
NUMÉRICA
• Esta opción de Polymath muestra una tabla de datos que se utiliza para la entrada,
manipulación y almacenamiento de datos numéricos.
• Los datos se almacenan en forma de columnas donde cada columna está asociada con un
nombre (variable) y se puede tratar por separado.
• Se puede realizar una regresión (ajustando una línea recta, varias curvas y ecuaciones a los
datos usando técnicas de regresión lineal, lineal, polinomial y no lineal), analizados
(interpolados, diferenciados, integrados junto con varios cálculos estadísticos) y trazado.
INTERPOLACIÓN Y REGRESIÓN
NUMÉRICA
• La Tabla de Datos Polymath para la entrada y manipulación de datos se muestra arriba en la tabla de arriba. La fila
superior de celdas contiene los nombres de columna (los nombres predeterminados son C01 a C30). El resto son
celdas de datos, que tienen 300 filas por 30 columnas (valores predeterminados).
• La dirección de la celda activa actual se muestra en la esquina superior izquierda de la tabla de datos. Para
seleccionar las opciones de análisis numérico con los datos de las columnas de debe seleccionar los datos de
interés.
• Los nombres de columna se pueden cambiar usando el cambio "Nombre de columna" en el menú desplegable
"Columna" que se muestra a la izquierda. Un nombre debe comenzar con una letra, y puede contener letras y
dígitos. Consulte la sección Variables y expresiones para obtener más reglas de nomenclatura de columnas.
• Al ingresar los datos, es importante colocar los datos en todas las filas. Dejar filas vacías (en particular la primera
fila) puede impedir que POLYMATH identifique las columnas en regresión y / o análisis.
INTERPOLACIÓN Y REGRESIÓN
NUMÉRICA
• El contenido de una columna se puede definir como función de columnas previamente
definidas. Para definir tal columna, seleccione "Definición de columna" en el menú desplegable
"Columna". Las columnas deben ser previamente nombradas antes de que puedan ser usadas
en una expresión.
• También puede establecer o cambiar el nombre de columna en el lado izquierdo del signo "="
e ingresar la definición de columna en forma de una expresión donde la expresión puede
contener constantes y nombres de columnas previamente definidas.
Valor de la Celda o
Ecuación de Definición
Nombre de la Columna
INTERPOLACIÓN Y REGRESIÓN
NUMÉRIC A
Exporta a Excel los Resultados
Sitúa los Resultados del modelo de
salida en la tabla Daa.
Muestre un informe que incluye el modelo
de regresión los valores numéricos ,los Se visualiza el grafico de la curva
intervalos de confianza de los parámetros, calculada (o puntos) y los puntos de
etc. datos.
Se visualiza un gráfico residual que muestra
la desviación entre los datos y los valores
calculados de la variable dependiente.
INTERPOLACIÓN Y REGRESIÓN
NUMÉRIC A – REGRESIÓN NO
LINEAL
Luego debemos seleccionar una
rutina de análisis de regresión
apropiada.
Aquí se ingresa el modelo que
desea de regresión. Proporcionar conjeturas iniciales
para los parámetros de su modelo,
en este caso A, B y C. (La solución
que proporciona Polymath puede
ser muy sensible a las suposiciones
de valores iniciales, por lo que si
La primera solución de regresión
no es muy buena, es posible que
desee cambiar las conjeturas
iniciales y volver a ejecutar el
regresión).
INTERPOLACIÓN Y REGRESIÓN
NUMÉRICA
Y=f(X)
X Y
1 0.00422922
5 0.003944 • Escribir los valores deseados en las columnas.
10 0.0038219
• Clic en ‘Regression’
20 0.00369617
40 0.00356189
• Definir la Variable dependiente (Y) e
60 0.0034656 independiente (X).
100 0.00334169 •Clic en “Report”
200 0.00317108 •Selecciona el grado de interpolación
400 0.00299625 •Clic en la flecha morada
760 0.00283086 -3
x 10
5
4.5
Y(X=50) = ?
y
3.5
2.5
0 500 1000
x
INTEGRACIÓN NUMÉRICA
Y=f(X)
X Y
1 0.00422922
5 0.003944 -3
x 10
10 0.0038219 5
20 0.00369617
4.5
40 0.00356189
60 0.0034656 4
100 0.00334169
y
200 0.00317108 3.5
400 0.00299625 3
760 0.00283086
2.5
0 500 1000
x
Área bajo la curva= =?
DIFERENCIACIÓN NUMÉRICA
Y=f(X)
X Y
-3
1 0.00422922 5
x 10
5 0.003944
10 0.0038219 4.5
20 0.00369617
4
40 0.00356189
y
60 0.0034656 3.5
100 0.00334169
200 0.00317108 3
400 0.00299625 2.5
760 0.00283086 0 500 1000
x
=?
PROBLEMA
La reacción de cloruro de trifenilmetilo (A) y metanol (B), se
llevo a cabo en una solución de benceno y piridina a 25°C. La
piridina reacciona con HCl que después precipita como
clorhidrato de piridina, de tal modo que la reacción es
irreversible. Los datos de concentración – tiempo se obtuvieron
en un reactor intermitente. La concentración inicial de metanol
era de 0.5 mol/dm3. En un conjunto de experimentos, se
encontró que el orden de reacción respecto al metanol es de
primer orden.
Determinar el orden de reacción y la constante de reacción.
(𝐶6 𝐻5 )3 𝐶𝐶𝑙 + 𝐶𝐻3 𝑂𝐻 → (𝐶6 𝐻5 )3 𝐶𝑂𝐶𝐻3 + 𝐻𝐶𝑙
𝐴+𝐵 → 𝐶+𝐷
T, min 0 50 100 150 200 250 300
Ca, (mol/dm3) x 50 38 30.6 25.6 22.2 19.5 17.4
103
PLANTEAMIENTO
𝐴+𝐵 → 𝐶+𝐷
𝛽
−𝑟𝐴 = 𝑘𝐶𝐴𝛼 𝐶𝐵
𝐶𝐵 = 𝐶𝐵𝑜
𝑘 ′ = 𝑘𝐶𝐵
𝜶
−𝒓𝑨 = 𝒌′𝑪𝑨
METODO INTEGRAL
−𝒓𝑨 = 𝒌′𝑪𝜶𝑨
• Orden Cero: 𝐶𝐴 = 𝐶𝐴0 − 𝑘𝑡
• Orden Uno: ln 𝐶𝐴 = ln 𝐶𝐴0 − 𝑘𝑡
1 1
• Orden Dos: = + 𝑘𝑡
𝐶𝐴 𝐶𝐴0
RESULTADOS
1
= 20.11753 + 0.124794𝑡
𝐶𝐴
𝑑𝑚3
𝑘′ = 0.124794
𝑚𝑜𝑙 𝑚𝑖𝑛
𝑘′
𝑘=
𝐶𝐵
2
0.12479 𝑑𝑚3 1
𝑘= = 0.2495 .
0.5 𝑚𝑜𝑙 𝑚𝑖𝑛
−𝒓𝑨 = 𝟎. 𝟐𝟒𝟗𝟓𝑪𝟐𝑨 𝑪𝑩
METODO DIFERENCIAL
• Ajustas el polinomio: 𝐶𝐴 = 𝑓(𝑡)
• Linealizas la ecuación:
−𝒓𝑨 = 𝒌′𝑪𝜶𝑨
𝑑𝐶𝐴
ln − = ln 𝑘´ + 𝛼 ln 𝐶𝐴
𝑑𝑡
• Encuentras los parámetros: 𝑘 𝑦 𝛼
RESULTADOS
𝑑𝐶𝐴
ln − = −1.925 +2.049 ln 𝐶𝐴
𝑑𝑡
𝑑𝑚3
ln 𝑘′ = −1.925 ′
𝑘 = 0.1459
𝑚𝑜𝑙 𝑚𝑖𝑛
2
𝑘′ 0.1458 𝑑𝑚3 1
𝑘= 𝑘= = 0.292 .
𝐶𝐵 0.5 𝑚𝑜𝑙 𝑚𝑖𝑛
−𝒓𝑨 = 𝟎. 𝟐𝟗𝟐𝑪𝟐.𝟎𝟒𝟖𝟓
𝑨 𝑪𝑩
REGRESION NO LINEAL
−𝒓𝑨 = 𝒌′𝑪𝜶𝑨 • Elegir una variable conocida, en
este caso el tiempo
𝒅𝑪𝑨 • Desarrollar el modelo
− = 𝒌′𝑪𝜶𝑨
𝒅𝒕 • Debemos suponer k y α
• Iterar hasta llegar al valor
mínimo de varianza
𝒅𝑪𝑨 ′ • Es importante conocer los
− 𝜶 = 𝒌 𝒅𝒕
𝑪𝑨 valores iniciales para iterar
𝐶𝐴 𝑡
𝑑𝐶𝐴
න − 𝛼 = න 𝑘′𝑑𝑡
𝐶𝐴𝑜 𝐶𝐴 0 𝟏
𝒕= (𝟎. 𝟎𝟓)𝟏−𝜶 −𝑪𝟏−𝜶
𝑨
𝒌 𝟏−𝜶
1 1−𝛼
𝑡= 𝐶𝐴𝑜 − 𝐶𝐴1−𝛼
𝑘 1−𝛼
RESULTADOS
′
𝑑𝑚3
𝑘 = 0.1467193
𝑚𝑜𝑙 𝑚𝑖𝑛
2
0.1467 𝑑𝑚3 1
𝑘= = 0.293 .
0.5 𝑚𝑜𝑙 𝑚𝑖𝑛
−𝒓𝑨 = 𝟎. 𝟐𝟗𝟑𝑪𝟐.𝟎𝟒𝟒𝟕
𝑨 𝑪𝑩
EJERCICIO 1
Calcule la Energía de activación para la descomposición de
cloruro de benceno diazonio para dar cloro-benceno y
nitrógeno:
k, s-1 T, K 𝐸𝑎
0.00043 313 −
𝑘 = 𝑘0 𝑒 𝑅𝑇
0.00103 319
0.0018 323 𝐸𝑎 1
ln 𝑘 = ln 𝐾𝑜 −
0.00355 328 𝑅 𝑇
0.00717 333
EJERCICIO 2
t, min Ca, mol/L
0 1
Para una reacción a volumen constante A→B se ha 1 0.9
determinado los siguientes datos experimentales que 4 0.77
se muestran en la siguiente tabla. Determinar el 16 0.54
orden de reacción y la constante de reacción. 28 0.39
50 0.26
80 0.19
100 0.16
160 0.1