Procesos Químicos

CARACTERÍSTICAS DEL MODELADO DE PROCESOS QUÍMICOS

Desde el punto de vista de los sistemas que se estudian, los modelos dinámicos se clasifican en estocásticos y determinísticos; los modelos estocásticos se basan en probabilidad de eventos que conlleva a variables aleatorias, por otro lado tenemos a los modelos determinísticos que se basa en la existencia de la correlación entre variables donde se aplican leyes físicas, químicas, biológicas, físico-químicas etc, esta basado en leyes de conservación de energía así como en reacciones nucleares.

Dentro del planteamiento con ecuaciones algebraicas se tiene:

Balance de masa

Entrada = Salida

(en estado estacionario)

• Con acumulación:

dm = me - ms

• Sin acumulación:

me = ms

Ee = Ee

      m = de masa

      e = de energía

      en este caso, es un modelo estático

Recapitulando un poco, la simulación puede clasificarse en determinística o estocástica:

Determinística: Las ecuaciones dependen de parámetros y variables conocidas con certeza; no existen leyes de probabilidad.

Estocástica: Ciertas variables estarán sujetas a incertidumbre.

El modelo estocástico esta sujeto a la ley de probabilidad de resultados por lo que no entra en los procesos químicos.

Por otro lado recordamos la simulación de eventos discretos, donde existen variables de interés que no tienen un comportamiento continuo como el diseño de plantas Batch multiproducto.

•  SECUENCIAL - MODULAR
  
  Los procesos con reciclos den ser descompuestos así como los cálculos se deben realizar unidad por unidad.
  El modelo de operación unitaria es un elemento básico que es construido a partir de balances de masa energía y momentum; se obtiene un conjunto de ecuaciones algebraicas no- lineales (debe ser compatible y determinado):
  
  donde:u = variable de entrada/salida
  x = estado interno de la variable (temperatura, presión, concentración, etc).
  d = variable dependiente de la geometría (volumen, área de intercambio de calor etc)
  p = propiedades físicas (entalpías, valores de k etc).
  
• 
  
  ORIENTADA A ECUACIONES
  Todas las ecuaciones algebraicas no lineales se integran a un único conjunto y se resuelven simultáneamente.
  El reciclo no es dificultad para el modelo.
  Se refiere mas a modelos dinámicos.
  La ecuación que lo representa:
  
  Para un estado estacionario, el término derivativo es igual a 0
  Los resultados se muestran en tablas o gráficas.
  
• MÓDULOS SIMULTÁNEOS
  Es una mezcla de los anteriores.
  Muestra un enfoque a los simuladores comerciales futuros.

¿POR QUÉ LOS PROCESOS QUÍMICOS SE CONSIDERAN DETERMINÍSTICOS, CONTINUOS, ESTÁTICOS Y/O DINÁMICOS?

Porque las variables endógenas y exógenas no se pueden tomar como datos al azar en los modelos determinísticos. Por lo que las relaciones de estas variables son exactas y la probabilidad queda fuera de ellos.

Porque si se refiere a un modelo continuo las variables dependen del tiempo y en la mayoría de los procesos químicos existirá una variable por mínima que sea que dependerá del tiempo.

Porque ya sea en modelos estáticos (no dependen del tiempo) o dinámicos (se toma en cuenta el tiempo) existirá en un proceso.

¿SE PUDIERAN CONSIDERAR MODELADOS ESTOCÁSTICOS DE LOS PROCESOS QUÍMICOS?

No se puede considerar debido a que en algún elemento cuando no se conoce con anticipación se esta incorporando a la incertidumbre. En un proceso ningún dato se puede tomar al azar debido a que todo esta debidamente controlado.

Un proceso estocástico tiene numerosos conceptos de probabilidad y estadística que sirven para realizar grandes series de muestreos.

CONCLUSIONES
En los procesos químicos difícilmente se darán los procesos estocásticos debido a la incertidumbre que maneja, en un proceso a nivel Industria todo debe estar regularizado y monitorizado porque cualquier falla implica pérdidas para la empresa.

⨠ REFERENCIAS

• Apuntes de clase Scenna Nicolás José. (1999).
• Modelado, simulación y optimización de procesos químicos. Argentina: ISBN

Pruebas de bondad de ajuste

PRUEBA DE BONDAD DE AJUSTE JI CUADRADA

Se utiliza para encontrar la distribución de probabilidad de una serie de datos. La metodología se basa en lo siguiente:

1. Se colocan los n datos históricos en una tabla de frecuencias de 
   

   

   intervalos. Se obtiene la frecuencia observada en cada intervalo i(FOi). Se calcula la media y la varianza de los datos.
2. Se propone una distribución de probabilidad de acuerdo con la forma de la tabla de frecuencias obtenida en el paso 1.
3. Con la distribución propuesta, se calcula la frecuencia esperada para cada uno de los intervalos (FEi) mediante la integración de la distribución propuesta y su posterior multiplicación por el número total de datos.
4. Se calcula el estimador:
   

   
5. Si el estimador C es menos o igual al valor correspondiente (jicuadrada) con m-k-1 grados de libertad (k = número de parámetros estimados de la distribución) y a un nivel de confiabilidad de 1-alfa, entonces no se puede rechazar la hipótesis de que la información histórica sigue la distribución propuesta en el punto 2.

Ejemplo:

Mediante la prueba x^2 determine el tipo de distribución de probabilidad que sigue la demanda de automóviles a un nivel del 95%, si a través del tiempo se ha registrado el comportamiento consignado en la figura 1.

Obtenga la tabla de frecuencias de la figura 1 considerando 7 intervalos y cuantificando la frecuencia para cada uno de ellos:

Figura 1. Comportamiento de la demanda

La distribución de probabilidad esperada que se propone, observando los datos de la columna FO, es una distribución uniforme entre  y  automóviles por día, o sea:

Integrando la función f(x):

Donde
LI = límite inferior de cada intervalo.
LS = límite superior de cada intervalo.

Sustituyendo los valores de los límites para obtener F(x) y multiplicándolos por el total de datos, se obtiene FE para cada intervalo.

Calculando el estadístico C con los datos de FEi y FOi, se obtiene

C = 4.092

El valor C = 4.092, comparado con el valor de la tabla de

indica que no podemos rechazar que los datos anteriores se comportan de acuerdo a una distribución uniforme entre 0 y 13 automóviles demandados por día con un nivel de confianza del 95%. Entonces,

PRUEBA DE AJUSTE DE BONDAD KOLMOGOROV-SMIRNOV

Esta prueba encuentra el tipo de distribución de probabilidad de una serie de datos, siendo más eficiente que la de (jicuadrada) ya que trabaja con la distribución de probabilidad acumulada. La metodología se basa en lo siguiente:

1. Se colocan n datos históricos en una tabla de frecuencias con 
   

   

   intervalos. Para cada intervalo se tendrá la frecuencia observada i(FOi). Se calcula la media y la varianza de los datos.
2. Se divide la frecuencia observada de cada intervalo por el número total de datos. A este resultado para obtener la probabilidad observada i(POi).
3. Se calcula la probabilidad acumulada observada de cada intervalo (PAOi) del paso 2.
4. Se propone una distribución de probabilidad de acuerdo con la forma de la tabla de frecuencias obtenida en 1.
5. Con la distribución propuesta se calcula la probabilidad esperada para cada uno de los intervalos (PEi) mediante la integración de la distribución propuesta.
6. Se calcula la probabilidad acumulada esperada (PAEi) para cada intervalo de clase.
7. Se calcula el valor absoluto entre PAOi y PEOi para cada intervalo y se selecciona la máxima diferencia, llamándola DM.
8. El estimador DM se compara con un valor límite con n datos y a un nivel de confiabilidad de 1-alfa. Si el estimador DM es menor o igual al valor límite, entonces no se puede rechazar que la información histórica sigue la distribución propuesta en el paso 4.

Ejemplo:

Mediante la prueba de Kolmogorov-Smirnov determine el tipo de distribución de probabilidad que siguen los datos del ejemplo anterior, con un nivel de confianza del 95%.

Obtenga la tabla de frecuencias, considerando 7 intervalos:

La distribución de probabilidad esperada que se propone, según los datos de la columna FO, es una distribución uniforme entre 0 y 13 automóviles por día, o sea:

Integrando la función f(x):

Donde:
LS = límite superior de cada intervalo.

Evaluando la ecuación anterior, se obtiene la tabla siguiente:

Al obtener la diferencia término a término entre PEA y POA, se tiene:

EL valor DM es igual a la máxima diferencia, o sea, 0.0694, que comparando contra el valor  , indica que los datos anteriores siguen una distribución uniforme entre 0 y 13 automóviles demandados por día, a un nivel de confianza del 95%. Por lo tanto,

⨠ REFERENCIAS

• Reza M. y García, E. (1996). Simulación y análisis de modelos estocásticos. (pp. 78-80). México: McGraw-Hill.