ABSORCIÓN
ABSORCIÓN
ABSORCIÓN
DATOS PERSONALES:
INTEGRANTES
G.H: 01-Q
SEMESTRE: 2014- A
INDICE
I. COEFICIENTES GLOBALES DE TRANSFERENCIA DE MASA ................................................................ 3
COEFICIENTES LOCALES DE TRANSFERENCIA DE MASA .................................................................. 8
II. MODELAMIENTO DE LA TRANSFERENCIA DE OXÍGENO EN MEDIOS BIOLÓGICOS ....................... 10
CONCEPTOS BÁSICOS .................................................................................................................... 10
Características reológicas. ......................................................................................................... 13
Solubilidad de gases .................................................................................................................. 14
Área interfacial gas – líquido ..................................................................................................... 14
Coeficiente volumétrico global de transferencia de masa ....................................................... 16
COEFICIENTES DE TRANSFERENCIA DE MASA ............................................................................... 18
Modelado de la absorción de dióxido de azufre en el contacto de membranas .......................... 20
Modelo basado en coeficientes de transferencia de materia .................................................. 22
Modelos Matemáticos de Transferencia de Masa en Deshidratación Osmótica ......................... 32
Modelo de Crank ....................................................................................................................... 34
Modelo de Magee ..................................................................................................................... 37
III. PROCESOS DIFUSIONALES ABSORCION GASES
ABSORCION SIN REACCION QUIMICA O ABSORCION FISICA......................................................... 39
ABSORCION CON REACCION QUIMICA O ABSORCION QUIMICA ................................................. 41
ELIMINACION DE SO2 CONTENIDO EN GASES DE ESCAPE DE UNA COMBUSTIÓN ..................... 43
ABSORCION DE NO2 EN DISOLUCION DE NaHS y NaOH ............................................................... 43
COEFICIENTES GLOBALES DE TRANSFERENCIA DE MASA
-(hl / hG)
P
PAG
PA = m.CA+b, ecuación de
la curva de equilibrio.
PAi M -(hl / hG), es constante en
un punto, no es a la largo
de una columna
PA*
Donde:
Luego:
𝑃𝐴𝐺 − 𝑃𝐴𝑖 𝑘𝑙
=−
𝐶𝐴𝑙 − 𝐶𝐴𝑖 𝑘𝐺
Y de:
𝑃𝐴𝐺 − 𝑃𝐴𝑖 𝑘𝑙
=−
𝐶𝐴𝑙 − 𝐶𝐴𝑖 𝑘𝐺
Se obtiene:
𝑘𝑙
𝑃𝐴𝐺 − 𝑃𝐴𝑖 = (𝐶 − 𝐶𝐴𝑙 )
𝑘𝐺 𝐴𝑖
Entonces:
𝑃𝐴𝑖 = 𝑚. 𝐶𝐴𝑖 + 𝑏
𝑃𝐴 ∗ = 𝑚. 𝐶𝐴𝑙 + 𝑏
Luego:
𝑘𝑙
𝑃𝐴𝐺 − 𝑃𝐴 ∗ = . (𝐶𝐴𝑖 − 𝐶𝐴𝑙 ) + 𝑚. (𝐶𝐴𝑖 − 𝐶𝐴𝑙 )
𝑘𝐺
𝑘𝑙
𝑃𝐴𝐺 − 𝑃𝐴 ∗ = ( + 𝑚) . (𝐶𝐴𝑖 − 𝐶𝐴𝑙 )
𝑘𝐺
Reemplazando en:
𝑘𝑙
𝑁𝐴 = 𝐾𝐺 (𝑃𝐴𝐺 − 𝑃𝐴 ∗ ) = 𝐾𝐺 . ( + 𝑚) . (𝐶𝐴𝑖 − 𝐶𝐴𝑙 ) … … . . (𝛼)
𝑘𝐺
Y como:
𝑁𝐴 = 𝑘𝑙 (𝐶𝐴𝑖 − 𝐶𝐴𝑙 ) … … … . (𝛽)
𝑘𝑙
𝑘𝐺 . ( + 𝑚) . (𝐶𝐴𝑖 − 𝐶𝐴𝑙 ) = 𝑘𝑙 (𝐶𝐴𝑖 − 𝐶𝐴𝑙 )
𝑘𝐺
𝑘𝑙
𝐾𝐺 . ( + 𝑚) . (𝐶𝐴𝑖 − 𝐶𝐴𝑙 ) = 𝑘𝑙 (𝐶𝐴𝑖 − 𝐶𝐴𝑙 )
𝑘𝐺
Entonces:
𝑘𝑙 𝑘𝑙
+𝑚 =
𝑘𝐺 𝐾𝐺
Expresión del coeficiente global de transferencia en función de los coeficientes
locales de transferencia:
1 1 𝑚
= +
𝐾𝐺 𝑘𝐺 𝑘𝐿
1 𝑘𝐿 + 𝑚𝑘𝐺
=
𝐾𝐺 𝑘𝐺 . 𝑘𝐿
Análogamente:
𝑁𝐴 = 𝐾𝐿 ( 𝐶𝐴∗ − 𝐶𝐴𝐿 )
De:
𝑃𝐴𝐺 − 𝑃𝐴𝑖 𝑘𝐿
=
𝐶𝐴𝑖 − 𝐶𝐴𝐿 𝑘𝐺
𝑘𝐺
( 𝐶𝐴𝑖 − 𝐶𝐴𝐿 ) = (𝑃 − 𝑃𝐴𝑖 )
𝑘𝐿 𝐴𝐺
De:
𝑃𝐴 = 𝑚 𝐶𝐴 + 𝑏
𝑃𝐴𝐺 = 𝑚 𝐶𝐴∗ + 𝑏
𝑃𝐴𝑖 = 𝑚 𝐶𝐴𝑖 + 𝑏
1
( 𝐶𝐴∗ − 𝐶𝐴𝑖 ) = (𝑃 − 𝑃𝐴𝑖 )
𝑚 𝐴𝐺
Luego:
𝐾𝐿 𝐾𝐿 𝑘𝐺
𝑁𝐴 = 𝐾𝐿 ( 𝐶𝐴∗ − 𝐶𝐴𝐿 ) = (𝑃𝐴𝐺 − 𝑃𝐴𝑖 ) + (𝑃𝐴𝐺 − 𝑃𝐴𝑖 )
𝑚 𝑘𝐿
𝐾𝐿 𝐾𝐿 𝑘𝐺
𝑁𝐴 = + (𝑃𝐴𝐺 − 𝑃𝐴𝑖 ) = 𝐾𝐺 (𝑃𝐴𝐺 − 𝑃𝐴𝑖 )
𝑚 𝑘𝐿
De donde:
𝐾𝐿 𝐾𝐿 𝑘𝐺
+ = 𝑘𝐺 … ( 𝐷𝑖𝑣𝑖𝑑𝑖𝑒𝑛𝑑𝑜 𝑝𝑜𝑟 ( 𝐾𝐿 ∗ 𝑘𝐺 )
𝑚 𝑘𝐿
𝐾𝐿 𝐾𝐿 𝑘𝐺 𝑘𝐺
+ =
𝑚 ∗ 𝐾𝐿 𝑘𝐺 𝑘𝐿 ∗ 𝐾𝐿 𝑘𝐺 𝐾𝐿 𝑘𝐺
1 1 1
= +
𝐾𝐿 𝑘𝐺 𝑚 𝑘𝐿
1 𝑘𝐿 + 𝑚 𝑘𝐺 1 𝑘𝐿 + 𝑚 𝑘𝐺
= =
𝐾𝐿 𝑘𝐿 𝑚 𝑘𝐺 𝑚 𝑘𝐿 𝑘𝐺
Haciendo:
𝑘𝐿 + 𝑚 𝑘𝐺 𝑘𝐿 𝑚 𝑘𝐺 1 𝑚
= + = +
𝑘𝐿 𝑘𝐺 𝑘𝐿 𝑘𝐺 𝑘𝐿 𝑘𝐺 𝑘𝐺 𝑘𝐿
1 𝑚 1
+ =
𝑘𝐺 𝑘𝐿 𝐾𝐺
Entonces:
1 1
=
𝐾𝐿 𝑚 𝐾𝐺
𝐾𝐿 = 𝑚 𝐾𝐺
El desarrollo es general, pueden ser dos líquidos irrisibles en lugar de una fase
liquida y una gaseosa; el desarrollo se hace para un punta en la columna (relación
puntual)
Para gas muy soluble “m” pequeño y en este caso 𝑘𝐺 ~ 𝑘𝐺
𝑁𝐴 𝑁𝐴
𝑁𝐴 − 𝑌𝐴𝐿 𝑁𝐴 − 𝑋𝐴𝐿
∑ 𝑁𝑖 ∑ 𝑁𝑖
𝑁𝐴 = ∗ 𝐹 ∗ 𝐿𝑛 ( )= ∗ 𝐹 ∗ 𝐿𝑛 ( )
∑ 𝑁𝑖 𝐺 𝑁𝐴 ∑ 𝑁𝑖 𝐿 𝑁𝐴
− 𝑌𝐴𝑖 − 𝑋𝐴𝑖
∑ 𝑁𝑖 ∑ 𝑁𝑖
Dónde:
∑ 𝑁𝑖 = 𝑁𝐴 + 𝑁𝐵 + 𝑁𝐶+⋯
La ecuación:
𝐹𝐿
𝑁𝐴 𝑁𝐴 𝐹𝐺
− 𝑌𝐴𝐿 − 𝑋𝐴𝐿
∑ 𝑁𝑖 ∑ 𝑁𝑖
=( )
𝑁𝐴 𝑁𝐴
− 𝑌𝐴𝑖 − 𝑋𝐴𝑖
∑ 𝑁𝑖 ∑ 𝑁𝑖
𝑁𝐴 𝑁𝐴
𝑁𝐴 − 𝑌𝐴∗ 𝑁𝐴 − 𝑋𝐴𝐿
∑ 𝑁𝑖 ∑ 𝑁𝑖
𝑁𝐴 = ∗ 𝐹𝑂𝐺 ∗ 𝐿𝑛 ( )= ∗ 𝐹𝑂𝐿 ∗ 𝐿𝑛 ( )
∑ 𝑁𝑖 𝑁𝐴 ∑ 𝑁 𝑖 𝑁𝐴
− 𝑌𝐴𝐺 − 𝑋𝐴∗
∑ 𝑁𝑖 ∑ 𝑁𝑖
𝑁𝐴 𝑁
(𝑁 𝐴 )
𝑁𝐴 𝑁
(𝑁 𝐴 )
(𝑁 ) − 𝑌𝐴∗
𝐴 ∗𝐹
∑ 𝑁𝑖 𝑂𝐺
𝐴 ∗𝐹
𝐺 ∑ 𝑁 𝑖
𝐴 ∗𝐹
𝐿
𝑒 = 𝑒 ∑ 𝑁𝑖 + 𝑚′ ∗ ∗ (1 − 𝑒 ∑ 𝑁𝑖 )
𝑁𝐴
− 𝑌𝐴𝐺
∑ 𝑁𝑖
𝑁𝐴 𝑁
(𝑁 𝐴 )
𝑁𝐴 𝑁
(𝑁 𝐴 )
1 ∑ 𝑁𝑖 − 𝑌𝐴
(𝑁 ) ∗
𝐴 ∗𝐹 𝐴 ∗𝐹 𝐴 ∗𝐹
𝑒 ∑ 𝑁𝑖 𝑂𝐿 = 𝑒 ∑ 𝑁𝑖
𝐺
+ ′′ ∗ ∗ (1 − 𝑒 ∑ 𝑁𝑖
𝐺
𝑚 𝑁𝐴
− 𝑌𝐴𝐺
∑ 𝑁𝑖
−𝑁𝐴 𝑁𝐴 𝑁𝐴
( 𝑁𝐴 ) −𝑌𝐴𝐺 ( 𝑁𝐴 )
( )𝐹
∑ 𝑁𝑖 𝑜𝐿
1 ∑ 𝑁𝑖 ( )𝐹
∑ 𝑁𝑖 𝐺
𝑒 = ( 𝑁𝐴 ) (1-𝑒 )
𝑚´´ −𝑌𝐴𝐿
∑ 𝑁𝑖
𝑁𝐴
∑ 𝑁 = 𝑁𝐴 ; =1
∑ 𝑁𝑖
𝑁𝐴 𝑁𝐴 −𝑁𝐴
1−𝑋𝐴𝐿 ( )
𝑒 𝐹𝑂𝐺
=𝑒 𝐹𝐺 + m´ ( 1−𝑌 )(1-𝑒 𝐹𝐿
𝐴𝐺
𝑁𝐴 𝑁𝐴 −𝑁𝐴
1 1−𝑌
𝑒 𝐹𝑂𝐿 = ( 𝐴𝐺) (1 − 𝑒 𝐹𝐺 ) + 𝑒 𝐹𝐿
𝑚´´ 1−𝑋 𝐴𝐿
B) Contradifusión equimolar
1 1 𝑚´
= +
𝐹𝑂𝐺 𝐹𝐺 𝐹𝐿
1 1 1
= +
𝐹𝑂𝐿 𝑚´´𝐹𝐺 𝐹𝐿
MODELAMIENTO DE LA TRANSFERENCIA DE OXÍGENO EN MEDIOS
BIOLÓGICOS
CONCEPTOS BÁSICOS
Para especies poco solubles en agua como oxígeno o hidrocarburos, las dos
concentraciones en equilibrio en la internase CGi y CLi de la película gaseosa y
líquida respectivamente, pueden ser relacionadas en forma lineal a través de una
Ley de partición:
ningún otro lugar del sistema. De esta forma pueden escribirse dos ecuaciones
iguales para la velocidad de transferencia de oxígeno:
Donde:
CA = solubilidad de A en la fase líquida [mol A / m3]
P = presión total en la fase gaseosa [atm]
PA = presión parcial de A en la fase gaseosa [atm]
H = constante de la ley de Henry [atm. m3 / mol A]
Ya= fracción molar de A en la fase gaseosa
En general, el área interfacial por unidad de volumen del líquido (a) puede
calcularse como:
Donde:
por:
Donde:
En sistemas aire - agua éste diámetro es generalmente del orden de 0,45 cm. La
viscosidad aparente µapp puede ser la viscosidad de los fluidos puros
Newtonianos (µapp ~ µL).
Modelado de la absorción de dióxido de azufre en el contacto de membranas
(1)
(2)
(4)
Que en el caso de sistemas con bajas concentraciones (sistemas diluidos),
(5)
Donde k es el coeficiente de transferencia de materia.
La elección entre ambos modelos representa un compromiso entre la profundidad
del conocimiento deseado y la base experimental. Trabajar con coeficientes de
difusión da lugar a la obtención delmodelo más general y fundamental, pero en
muchos casos, unmodelo fenomenológico y sencillo como el basado en
coeficientes detransferencia de materia puede ser igualmente útil para describir
elsistema.
En términos generales, en el modelo de difusión la variabledependiente
(concentración) puede cambiar con todas las variablesindependientes (posición,
tiempo), mientras que el modelo detransferencia de materia presenta valores
medios (e.g.concentración media logarítmica) (Cussler, 1997).
∗ ∗
Las concentraciones en la interfase del componente que difunde (𝐶𝑆𝑂2,𝑔 𝑦 𝐶𝑆𝑂2,𝑙 )
generalmente difieren pero se asume que estánrelacionadas mediante las leyes
del equilibrio termodinámico (Perryy Green, 1997). De esta forma, se considera
que el equilibriotermodinámico se alcanza en la interface gas-líquido cuando el
gasy el líquido entran en contacto. Si se asume que en la interface secumple la ley
∗ ∗
de Henry, las concentraciones 𝐶𝑆𝑂2,𝑔 𝑦 𝐶𝑆𝑂2,𝑙 ) estánrelacionadas mediante la
siguiente ecuación:
(9)
(11)
(11a)
(12a)
o expresado como resistencias,
(12b)
Para una membrana porosa, la transferencia de materia a través de la membrana
se calcula como el coeficiente de difusión efectiva, Deff, del gas a través de los
poros de la membrana dividido por el espesor de la membrana (Cussler, 1997;
Gabelman y Hwang, 1999):
(13a)
(13b)
Las difusividades del dióxido de azufre en el gas ocluido en los poros de la
membrana, DSO2,g, y en la fase líquida, DSO2,l,
El coeficiente de transferencia de materia referido a la fase gas fluyendo por el
exterior y paralelo a la fibra puede ser estimado a partir de la siguiente correlación,
con una desviación en torno al ±10% (Cussler, 1997: Yang et al, 2006):
(14)
(15)
Donde dh es el diámetro hidráulico; DSO2,b es el coeficiente de difusión en el
seno de la fase gas; L es la longitud de la fibra; y es la viscosidad cinemática.
En el caso del líquido, dado que fluye con flujo laminar por el interior de la fibra, la
ecuación de Lévêque permite estimar el coeficiente de transferencia de materia kl,
con una desviación de ±1% (Cussler, 1997; Yang et al, 2006):
(16)
Donde di es el diámetro interior de la fibra y DSO2,l es el coeficiente de difusión en
la fase líquida.
Resultados y discusión
El coeficiente global de transferencia de materia promedio, calculado
experimentalmente en función de la concentración media logarítmica de dióxido de
azufre en la fase gas entre la entrada y la salida del contactor de membranas,
presenta un valor de Koverall = (1.10±0.11) x10-5 m\s-1, empleando
N,Ndimetilanilina como líquido de absorción.
Este coeficiente está relacionado con los coeficientes individuales de transferencia
de materia según la Ecuación (11) o (12), dependiendo de si se considera que los
poros de la membrana están llenos de gas o de líquido, respectivamente. Los
coeficientes referidos a la fase gas, kg, la fase líquida, kl, y la membrana, kmg y
kml, son estimados a partir de las Ecuaciones (13) a (16) y de los valores de los
parámetros listados en la Tabla 1. Los valores obtenidos se muestran en la
Tabla1.
Tabla 1. Parámetros utilizados en las Ecuaciones (13) a (16).
Coeficiente de difusión de SO2 en N,N-dimetilanilina, estimado según el Apéndice
A.
a. Coeficiente de difusión de SO2: Knudsen, bulk y global en aire, estimados
según el Apéndice A.
b. Viscosidad cinemática de la mezcla SO2-aire, estimada utilizando
correlaciones para viscosidad dinámica de Poling (2001) (Poling y
Prausnitz, 2001).
c. Estimación según el proceso industrial para una corriente gaseosa con
5%vol. De SO2, según Apéndice B. El dato para SO2 puro se muestra entre
paréntesis (van Dam et al, 1997)
(17a)
(17b)
Donde W es la fracción de mojado, definida como la relación entre la
longitud de poro mojado por el líquido (x) y el espesor total de la membrana
(δ = 5*10-4 m).
(18a)
Las difusividades efectivas para cada líquido se indican en la Tabla 4. Los valores
obtenidos tienen un orden de magnitud típico de sistemas donde el líquido se
encuentra ocluido en una matriz sólida (10-8 – 10-12 m2\s-1) (Perry y Green,
1997), lo cual refuerza el sentido de las conclusiones previamente obtenidas.
Introducción
𝑗 𝑗 ∞
𝑀0 − 𝑀𝑡 8 2
𝜋 2 𝐹𝑜
𝑗 𝑗
= 1−∑ exp (−(2𝑛 + 1) ) … (1)
𝑀0 − 𝑀∞ (2𝑛 + 1)2 𝜋 2 4
𝑛=0
𝑗
𝐷𝑒
Donde el número de Fourier (Fo) está dado por: 𝐹𝑜 = 𝑡
𝑙2
𝑗 𝑗 ∞
𝑀0 − 𝑀𝑡 𝑛
𝑗 𝑗
= 2(𝐹𝑜)0.5 (𝜋 −0.5 + 2 ∑(−1)𝑛 𝑖𝑒𝑟𝑓𝑐 ) … (2)
𝑀0 − 𝑀∞ 𝐹𝑜
𝑛=1
𝑗 𝑗
𝑀0 − 𝑀𝑡 8 𝜋 2 𝐹𝑜
𝑗 𝑗
=1− exp (− ) … (3)
𝑀0 − 𝑀∞ 𝜋2 4
𝑗 𝑗
𝑀0 − 𝑀𝑡 𝐹𝑜 0.5
𝑗 𝑗
= 2 ( ) … (4)
𝑀0 − 𝑀∞ 𝜋
Este modelo fue propuesto por Hawkes y Flink pero varios autores lo atribuyen a
Magee, quien hizo algunas modificaciones.
𝑀0 − 𝑀𝑡
= 𝑘𝑡 0.5 + 𝑘0
𝑀𝑡0