INSTITUTO TECNOLOGICO DE ORIZABA

PROCESOS DE SEPARACION

Examen 2 (1a. Parte)

TEMA:
MTODOS RIGUROSOS PARA EL CLCULO DE
OPERACIONES DE SEPARACIN DE MEZCLAS
MULTICOMPONENTES

Mtodo de la suma de caudales (SR). Aplicacin a absorcin y

desorcin

INSTRUCTOR:
DR. JUAN MANUEL MENDEZ CONTRERAS

PRESENTA:
I.Q ADRIANA SANOVAL RAMOS
I.Q LIZBETH CONTRERAS ROMERO
I.Q ALFONSO FLORES MEZA
I.Q IGNACIO SANCHEZ BAZAN
NOVIEMBRE 2004

1
 INSTITUTO TECNOLOGICO DE ORIZABA
PROCESOS DE SEPARACION

MTODOS RIGUROSOS PARA EL CLCULO DE OPERACIONES DE

SEPARACIN DE MEZCLAS MULTICOMPONENTES

Mtodo de la suma de caudales (SR). Aplicacin a absorcin y desorcin

INTRODUCCIN

Hasta 1950, los clculos correspondientes a los problemas de rectificacin tenan que
realizarse a mano. Aunque se dispona de mtodos rigurosos, eran difciles de aplicar,
incluso en el caso de columnas muy pequeas. En esta poca, los mtodos aproximados
constituan la principal herramienta de diseo y los rigurosos slo se aplicaban para
columnas muy pequeas o para comprobar el diseo final. La aparicin de los
ordenadores invirti la situacin: los mtodos rigurosos, cuya aplicacin hasta el
momento requera varios das de trabajo, o incluso semanas, suponan ahora muy poco
tiempo de clculo con el ordenador y por tanto no tena sentido aceptar como buenas las
inexactitudes e imprecisiones de los mtodos aproximados. En la actualidad, los mtodos
rigurosos han pasado a ser la principal herramienta de diseo, quedando como nico
papel de los mtodos aproximados el proporcionar una estimacin inicial para los
clculos rigurosos y eliminar las opciones menos adecuadas.

El diseo final del equipo de etapa mltiple requiere una determinacin rigurosa de la
temperatura, la presin, los caudales y composiciones de las corrientes y las velocidades
de transferencia de calor para cada etapa. Dada la elevada no linealidad de las ecuaciones
que describen el proceso, es necesario resolver el problema general por el uso de
procedimientos iterativos y los diversos mtodos difieren en la seleccin del conjunto de
variables independientes.

Histricamente, los primeros intentos para la resolucin del sistema de ecuaciones

formado por los balances de materia, las relaciones de equilibrio, los sumatorios de
fracciones molares y los balances de entalpa (sistema de ecuaciones MESH) fueron los

2
 INSTITUTO TECNOLOGICO DE ORIZABA
PROCESOS DE SEPARACION
mtodos clsicos de clculo etapa a etapa y ecuacin a ecuacin de Lewis-Matheson
(1932) y Thiele-Geddes (1933), aplicables a columnas convencionales (se introduce un
nico alimento y se extraen dos productos, el destilado por la cabeza y el residuo por la
cola). Sin embargo, la mayor parte de los mtodos modernos son mtodos componente a
componente, haciendo referencia este tipo de designacin a la forma en que se agrupan
y se resuelven los sistemas de ecuaciones.

En la actualidad se utilizan mtodos que resuelven las ecuaciones MESH a partir de un

conjunto de ecuaciones lineales simultneas para las N etapas, agrupadas para cada
componente. Estos mtodos se pueden clasificar dentro de cuatro grandes grupos:

- Mtodos de punto de burbuja (BP), que calculan la temperatura de los pisos por
resolucin de la ecuacin del punto de burbuja.

- Mtodos de suma de caudales (SR), que usan los balances de energa para obtener la
temperatura de los pisos.

- Mtodos de Newton 2N, que calculan las temperaturas y el caudal total al mismo
tiempo, pero las composiciones en un paso separado.

Estas tres clases de mtodos separan las ecuaciones MESH en grupos, que resuelven a
travs de una serie de pasos.

- Mtodo de Newton global o de correccin simultnea (SC), que resuelve todas las
ecuaciones MESH simultneamente. Recientemente se han desarrollado otros mtodos
(mtodos inside-out, mtodos de relajacin, mtodos de homotopa, mtodos de no
equilibrio, etc.) que intentan solucionar a los problemas de convergencia que presentan
los mtodos anteriores, as como dar cobertura a la resolucin de problemas de
separacin de mezclas cada vez ms complejas.

3
 INSTITUTO TECNOLOGICO DE ORIZABA
PROCESOS DE SEPARACION
En general, los mtodos modernos estn ya programados, son rpidos y requieren un
espacio mnimo de memoria en un ordenador, aunque normalmente estn limitados a un
nmero relativamente reducido de posibilidades para la especificacin de variables.

Algoritmo de la matriz tridiagonal

Frecuentemente, los sistemas de ecuaciones que resultan al aplicar los mtodos
componente a componente suelen dar lugar a una matriz de coeficientes en forma de
matriz tridiagonal, y se resuelven mediante un mtodo de eliminacin progresiva que
recibe el nombre de algoritmo de Thomas.

Las ecuaciones MESH (M = balance de materia, E = equilibrio, S = sumatorio de

fracciones molares o msicas y H = balance de entalpa) para una etapa j pueden
escribirse:

1. M- balances de materia para cada componente (c ecuaciones para cada etapa):

M ij = L j 1 + Vi , j +1 + FjZ ij (L j + U j )X ij (V j + W j )Yij = 0 (1)
2. E- relaciones de equilibrio entre fases para cada componente (c ecuaciones por etapa):

Eij = Yij K ij X ij = 0 (2)

3. S - Sumatorios de las fracciones molares (una para cada etapa):
C (3)
(S Y ) j = Yij 1 = 0
i =1
C (4)
(S X ) j = X ij 1 = 0
i =1

4. H- Balance de energa (uno para cada etapa):

H j = L j 1 HL j 1 + V j +1 HV j +1 + FjHF j (L j + U j )HL j (V j + W j )HV j Qj = 0 (5)
Si se modifican las ecuaciones MESH de forma que en las ecuaciones M se sustituyen los
yij por Kijxij y los Lj por su valor en funcin de Wj, Uj y Vj dado por el balance global
de materia, se puede obtener la siguiente ecuacin para cada componente y etapa, donde
se ha omitido el subndice i para los trminos B, C y D:
A j X i , j 1 + B j X i , j + C j X i , j +1 = D j (6)

4
 INSTITUTO TECNOLOGICO DE ORIZABA
PROCESOS DE SEPARACION
donde, si se tiene en cuenta que Un balance global entre la etapa j-1 y la cabeza de la
columna proporciona:

(7)

(8)

(9)

siendo el entorno seleccionado para realizar los balances aquel en el que entran las
corrientes

Por tanto, en la ecuacin (6), los valores de los coeficientes son:

(10)

(11)

(12)

(13)

Si las ecuaciones M modificadas (6) se agrupan por componentes, pueden escribirse

como una serie de c sistemas de ecuaciones, uno para cada componente, en los que las
matrices de coeficientes son matrices tridiagonales y donde la variable de salida para
cada sistema es la composicin xi para ese componente en toda la cascada en
contracorriente de N etapas:

5
 INSTITUTO TECNOLOGICO DE ORIZABA
PROCESOS DE SEPARACION

(14)

Los coeficientes Bj y Cj dependen exclusivamente de las variables de tanteo Tj y Vj,

eliminndose de esta manera en el clculo de los coeficientes la dependencia de la
composicin.

La resolucin del conjunto de ecuaciones linealizadas (14) se simplifica mucho con la

aplicacin del algoritmo de Thomas, que es un mtodo Gaussiano de eliminacin en el
que se procede inicialmente a una eliminacin progresiva comenzando en la etapa 1 y
operando hasta alcanzar la etapa N, aislando finalmente xiN. Se obtienen as otros valores
de xij comenzando con xi, N-1 mediante una sustitucin por el camino inverso.

Las ecuaciones utilizadas en el algoritmo de Thomas son las siguientes: para la etapa 1, la
ecuacin (6) es de la forma B1xi1 + C1xi2 = D1, y puede resolverse para xi1 en funcin
de xi2 para dar:

(15)

Haciendo el cambio de variable:

(16)
(17)

6
 INSTITUTO TECNOLOGICO DE ORIZABA
PROCESOS DE SEPARACION
Por tanto, los coeficientes de la matriz se transforman en B11, C1p1 y D1q1,
donde quiere decir reemplazado por. Slo se necesita conservar los valores de p1 y
q1.
Para la etapa 2:

(18)

tomando

(19)

resulta
(20)
De esta forma, A20, B21 y C2p2 y D2q2. Slo se necesita conservar los
valores de p2 y q2. De forma general, se puede definir:

(21)

por tanto
(22)
Donde Aj0, Bj1 y Cjpj y Djqj. Slo se necesita conservar los valores de pj y qj.
De esta forma, partiendo de la etapa 1, se van calculando de forma progresiva los valores
de pj y qj en el siguiente orden: p1, q1, p2, q2,..., pN-1, qN-1, qN. En la etapa N, la
ecuacin (22) aisla xiN como:

(23)
se calculan entonces, a partir de (22), valores sucesivos de xi sustituyendo a lo largo del
clculo inverso en la forma:
(24)

El algoritmo de Thomas es de gran eficacia y es superior a las rutinas alternativas de

inversin de matrices.

7
 INSTITUTO TECNOLOGICO DE ORIZABA
PROCESOS DE SEPARACION

MTODO COMPONENTE A COMPONENTE

Mtodo de la suma de caudales (SR). Aplicacin a absorcin y Desorcin

Los componentes qumicos presentes en la mayora de los absorbedores y desorbedores

cubren un intervalo relativamente amplio de volatilidades, por lo que la utilizacin del
mtodo BP para la resolucin de las ecuaciones MESH suele fallar ya que, como se ha
dicho, el clculo de la temperatura de las etapas mediante el clculo del punto de burbuja
es demasiado sensible a la composicin de la fase lquida y el balance energtico de la
etapa es mucho ms sensible a la temperatura de las etapas que a los caudales interetapas.
En estos casos se puede aplicar con resultados satisfactorios el mtodo de la suma de
caudales (SR), desarrollado por Burninham y Otto, junto con la formulacin de la matriz
tridiagonal para resolver las ecuaciones M modificadas.

8
 INSTITUTO TECNOLOGICO DE ORIZABA
PROCESOS DE SEPARACION
La figura muestra el algoritmo para el mtodo SR de Burningham-Otto, donde las
especificaciones del problema son las condiciones y etapa de localizacin de todas las
alimentaciones, presin en cada etapa, caudales totales de todas las corrientes laterales,
caudales de transferencia de calor desde o hacia cualquier etapa, y nmero total de etapas.

Para iniciar los clculos se supone un conjunto inicial de variables de tanteo Tj y Vj. Para
la mayora de los problemas es suficiente con suponer un conjunto de valores Vj basados
en la consideracin de caudal molar constante interetapas, operando desde la cola del
absorbedor utilizando los caudales de las alimentaciones y de las corrientes laterales de
fase vapor especificadas.

Generalmente puede obtenerse un conjunto adecuado de valores Tj iniciales a partir de

los valores supuestos para las etapas de cabeza y cola, suponiendo una variacin lineal
con el nmero de etapa.

Los valores de x ij se obtienen resolviendo (14) mediante un algoritmo de Thomas. Los

valores as obtenidos se utilizan directamente para calcular nuevos valores de Lj mediante
la ecuacin de suma de caudales:
c
L(jk +1) = L(jk ) X ij (25)
i =1

donde los valores de L(kj ) . Los correspondientes valores de Vj(k+1) se determinan a partir
de:
N
V j = L j 1 LN + (Fm Wm U m ) (26)
m= j

Se calculan a continuacin los valores normalizados de xij a partir de.

(X ) ij normalizado = c
X ij
(27)
X
i =1
ij

Los valores de yij correspondientes se calculan por el equilibrio, segn (2). Finalmente,
el perfil de temperatura para la siguiente iteracin se obtiene por resolucin del balance
de entalpa en cada etapa.

9
 INSTITUTO TECNOLOGICO DE ORIZABA
PROCESOS DE SEPARACION

Wang y Henke sugieren el siguiente criterio basado exclusivamente en los sucesivos

valores de T:

= (T j( k ) T j( k 1) ) 0.01N
N
2
(28)
j =1

EJEMPLO DEMOSTRATIVO:
Realcese una iteracin segn el mtodo SR para la columna de destilacin
propuesta.

10