Principios de Transferencia de Masa (Difusion Molecular)
Principios de Transferencia de Masa (Difusion Molecular)
Principios de Transferencia de Masa (Difusion Molecular)
FRANCISCO DE MIRANDA
AREA DE TECNOLOGIA
PROGRAMA DE INGENIERIA QUÍMICA
OPERACIONES UNITARIAS III
I. Transporte de Masa.
Introducción
El transporte de masa ocurre cuando un constituyente (de un sistema de varios
componentes) viaja desde una región de alta concentración a otra de baja
concentración.
La mayoría de los problemas típicos de Ingeniería Química caen en el campo
del transporte de masa, por consiguiente, la característica que distingue a un
ingeniero químico es su habilidad para diseñar y operar equipo en el cual se
preparan reactivos, o donde se da lugar a reacciones químicas y son realizadas
separaciones de los productos resultantes.
El proceso fundamental de transferencia, esto es, transferencia de masa ,
interviene en los procesos de destilación, absorción, secado, extracción líquido-
líquido, adsorción y procesos de membrana. Cuando se transfiere masa de una
fase a otra o a través de una sola fase el mecanismo básico es el mismo, ya sea
que se trate de gases, líquidos o sólidos. Esto también se demostró para la
transferencia de calor en la cual el transporte de calor por conducción obedece la
ley de Fourier en gases, líquidos y sólidos.
La ecuación general de transporte molecular como La transferencia de calor,
momento lineal y de masa se caracteriza o fundamenta por el mismo tipo general
de ecuación
fuerza impulsora
Velocidad de un Proceso de Transferencia =
resistencia
Esta ecuación puede escribirse como sigue: Para la difusión molecular de las
propiedades de momento lineal de calor y de masa:
A (2)
B B
B B B
B B
B B B
A
(1)
Ci = Ni / V
Concentración Másica (ρi)
ρi = Ci x Mi = mi / v
La concentración total de masa o densidad, ρ, es la masa total de la mezcla
contenida en la unidad de volumen, esto es:
n
ρ= Σ ρi
i=1
nt m n j m
C= = ∑ = ∑ ρi
V j =1 V J =1
m m m
ρ = mt = ∑ j = ∑ ρ i
V i =1 V i =1
mi ρ i m
ϖi = = Además ∑ϖ i = 1
mt ρ j =1
ci
Xi = ( Líquidos y Sólidos)
c
En el caso de gases se establece:
ni C i
yi = =
nt C
Para un gas ideal :
pi
yi = Pi = Presión parcial del componente
pt
PT = Presión Total
En el caso de líquidos no electrolitos (Solución como moléculas) se representa
por:
ci
xi =
c
Unidades de Velocidad
En un sistema de componentes múltiples, las diferentes especies se moverán
de manera normal a diferentes velocidades; por lo tanto, para evaluar la
velocidad de la mezcla de gases, se necesitan promediar la velocidad de cada
una de las especies presente.
m c i*v i m
.V = ∑ = ∑ y i*
i =1
m
i =1 c
∑Ci
i
i =1
Para un Gas
n ci * y i
V =∑
i =1 C
De acuerdo a la ley de Fick, una especie puede tener una velocidad relativa a
la masa o una velocidad molar media, solamente si existen los gradientes en
la concentración.
Flujo Molar
El flujo molar indica la variación del gradiente de concentración a través de
Fick fue quien primero postuló una relación empírica para este flujo molar y,
d cA
J AB = − D AB *
dz
Z = (m) ó (pie)
• Flujo Másico
dϖ A
= Es el gradiente de concentración en
dZ
función de la fracción de masa
dy A
J AZ = c A (v Az − V z ) = − cD AB
dz
después de reordenar esta ecuación, queda:
dy A
c A v Az = − cD AB + c AV z
dz
Vz =
1
(c A v Az + c B v Bz )
c
ó
c A V z = y A (c A v Az + c B v Bz )
Al sustituir esta ecuación en la expresión anterior, se obtiene una nueva
ecuación;
+ y A (c A v Az + c B v Bz )
dy A
c A v Az = cD AB (*6)
dz
como las velocidades componentes v Az , v Bz son velocidades relativa al eje
Cualquiera de estas cantidades puede ser una parte importante del fluido
molar total , NA .Cuando la ecuación (*7) se usa para describir una difusión
Las cuatros ecuaciones que definen los flujos JA, jA, NA y nA son enunciados
equivalentes de la ecuación de Fick de la rapidez o difusión.
∇y A J A = − cD AB ∇y A
JA Constante c
∇c A J A = − D AB ∇c A
La tabla ante indicada presenta un resumen de las formas equivalentes de la
ecuación de Fick de la rapidez.
4- Aplicación Ley de Fick para la difusión molecular.
Para arrancar aplicando la ley de Fick tómese primero en cuenta la difusión
molecular cuando la totalidad del fluido esta inmóvil, es decir, estacionario. La
difusión de las moléculas se debe a un gradiente de concentración.
Recordemos
que la ecuación general de la ley de Fick puede escribirse como sigue para
una mezcla binaria de A y B
dX A
J Az = cD AB (1),
dz
⎛ Kgmol A + B ⎞
donde: c= Es la concentración total de A y B en ⎜⎜ 3
⎟
⎟
⎝ m ⎠
cdx A = d (cx A ) = dc A
p A2
=0.2 atm. Calcule en unidades SI y CGS el flujo específico de He en estado
D AB (c A1 − c A 2 )
J Az = (3)
z 2 − z1
con base en la Ley de los gases ideales, pAV = nART, y
5 - Coeficiente de Difusión
La proporcionalidad de la ecuación de Fick, DAB, se conoce con el nombre de
coeficiente de difusión. Sus dimensiones fundamentales, que pueden
obtenerse a partir de la siguiente ecuación:
⎛ ⎞
⎜ ⎟
J A, z ⎛ M ⎞⎜ 1 ⎟ L 2
D AB = =⎜ ⎟ =
dc A ⎝ L 2 t ⎠⎜ M 1 ⎟ t
⎜ 3*L⎟
dz ⎝L ⎠
son idénticas a las dimensiones fundamentales de las otras propiedades de
transferencia : la viscosidad cinemática, v y la difusividad térmica ,∝ o su
razón equivalente ,K / ρ cp. La difusividad de la masa se ha dado en cm2/seg,
las unidades SI son m2/seg; o sea un factor 10-4 veces menor. En el sistema
inglés, pie2/hr, son las unidades utilizadas. Para hacer convecciones de uno
a otro de estos sistemas se utilizan las siguientes relaciones:
DAB (Cm2/seg) = 104 , DAB (Pie2/hr) = 3.87
DAB (m2/seg) DAB (cm2/seg)
NA = Kc (CL1 – CLi)
-
ING. CARMEN BROCK Página 15 de 35
TRANSFERENCIA DE MASA
DIFUSIÓN MOLECULAR
I.6- DIFUSIÓN MOLECULAR EN GASES
• Contradifusión Equimolar En Gases
En la figura.1. se muestra un diagrama para dos gases, A y B, a presión total
constante P, en dos cámaras grandes, conectadas por un tubo que sirve para que
se verifique la difusión molecular en estado estacionario. Una agitación en ambas
cámaras mantiene uniformes sus concentraciones. La presión parcial PA1> PA2 y
PB2>PB1. Las moléculas de A se difunden hacia la derecha y las de B hacia la
izquierda. Puesto que la presión total P es constante en todo el sistema, los moles
netos de A que se difunden hacia la derecha deben ser iguales a los moles netos de
B, que lo hacen hacia la izquierda. Si no fuera así, la presión total no se mantendría
constante. Esto significa que,
JAZ = - JBZ
pA1 pA2
1 2
pB1/ P JA* pB2/ P
*
JB
P
PA1
PB2
PB, ó P
PA2
PB1
vA
vAd vM
NA + NB
vM = (12)
c
Sustituyendo la ecuación (12) en la (10)
cA (
NA = JA + N A + N B) (13)
C
dx A c A (
N A = − cD AB + N A + N B) (14)
dz c
D AB dp A p A
NA = − + NA (17)
RT dz P
Reordenando e integrando,
⎛ pA ⎞ D dP A + c A
N A ⎜⎜1 − ⎟⎟ = − AB (18)
⎝ P⎠ RT dz c
z2 p
= − D AB A 2 dp A
NA ∫ dz ∫ (19)
RT p pA
A1 1 −
z1
P
D AB P − P A1
NA = Ln (20)
RT ( z 2 − z1) P − P A2
La ecuación (20) es la expresión final adecuada para calcular el flujo de A.
Sin embargo, con frecuencia se escribe también de otra forma. Primero se
define la media logarítmica de B inerte. Puesto que:
P= pA1+ pB1 = pA2 + pB2, pB1= P – pA1
pB2= P - pA2,
pB2 − pB2 p A1 − p A2
p BM = = (21)
⎛p ⎞ ⎡ (P − p A 2 ) ⎤
Ln⎜⎜ B 2 ⎟⎟ Ln ⎢ ⎥
⎝ p B1 ⎠ ⎣ (P − P A1) ⎦
NA=
D AB P ( p − p A2 ) (22)
RT (z 2 − z1) p BM A1
Ejemplo 3. Difusión de Agua a través de aire en reposo que no se
difunde
El agua en el fondo de un tubo metálico estrecho se mantiene a temperatura
constante de 293 ° K. La presión total del aire (que se supone seco) es 1.01325x
105 Pa (1.0 atm) y la temperatura es 293° K (20°C). El agua se evapora y se difunde
a través del aire en el tubo y la trayectoria de difusión z2-z1 tiene 0.1542 m (0.5 pie)
ING. CARMEN BROCK Página 21 de 35
TRANSFERENCIA DE MASA
DIFUSIÓN MOLECULAR
de longitud. Calcule la velocidad de evaporación en estado estacionario en lbmol /
h.pie2 y en kgmol / s.m2. La difusividad del vapor de agua a 293 °K y 1atm de
presión es 0.250x10-4 m2/s. Suponga que el sistema es isotérmico. Utilice unidades
SI y del sistema inglés.
AIRE(B)
PA2
z1 – z2
z
PA1
zf
AGUA (A)
p B 2 − p B1 100 − 0.9769
p Bm = = = 0.988atm = 1.001x105 p a
⎛p ⎞ ⎛ 1.00 ⎞
Ln⎜⎜ B 2 ⎟⎟ Ln⎜ ⎟
⎝ p B1 ⎠ ⎝ 0.9769 ⎠
Puesto que pB1 es cercano a pB2 puede emplearse la media lineal (pB1+pB2)/2 para
un valor muy similar a pBM.
Sustituyendo en la ecuación (22) con z2-z1 = 0.5 pie (0.1524m)
NA =
D AB
( p A1 − p A2 ) = (0.969)(1.0)(0.0231 − 0)
RT ( z 2 − z1) p bm (0.730)(528)(0.5)(0.988)
Kgmol
N A = 1.595 x10
−7
seg. m2
ρA z F
D AB P ( p A1 − p A 2 ) t F
∫ zdz = ∫ dt (25)
M A z0 RTp BM 0
Al despejar tF,
ING. CARMEN BROCK Página 23 de 35
TRANSFERENCIA DE MASA
DIFUSIÓN MOLECULAR
ρ A (z 2F − z 02 ) RTp BM
tF = (26)
2M A D AB P ( p A1 − p A 2 )
El método que se muestra en el ejemplo 4 se ha usado para determinar
experimentalmente la difusividad DAB. En este experimento, la longitud de
trayectoria inicial z0 se mide en t=0 y también la final zF en tF. Así, la ecuación
(26) se usa para calcular DAB.
• .- Difusión a través de un área de corte transversal variable.
NA
NA = (27)
A
Donde NA es kg mol de A que se difunde por segundo o kgmol / s. En
NA
NA = (28)
4πr
2
dr r1
r1 r NA r2
pA2
pA1 NA
z
1 2
(a) (b)
NA dp A
NA = = − D AB (29)
4πr
2
RT ⎛ pA ⎞
⎜1 − P ⎟dr
⎝ ⎠
Nótese que ∂r ha reemplazado a dz. Reordenando e integrando entre r1 y un
punto r2 a gran distancia,
p
N A r 2 dr D AB A 2 dp A
∫ 2 =− ∫ (30)
4π r 1 r RT p ⎛ p A1 ⎞
A1 ⎜1 −
⎝ P ⎟⎠
N A ⎛⎜ 1 − 1 ⎞⎟ = D AB P Ln P − p A2 (31)
4π ⎝ r1 r 2 ⎠ RT P − p A1
2 D AB (
N A1 = c A1 − c A 2 ) (33)
DI
Esta ecuación también se puede usar par líquidos, donde DAB es la difusividad
de A en el líquido.
Ejemplo 5. Evaporación de una esfera de naftaleno
Una esfera de naftaleno con radio de 2.0 mm está suspendida en gran
volumen de aire estacionario a 318 °K y 1.01325 x 105 Pa (1atm). Se puede
suponer que la temperatura superficial del naftaleno es 318 °K y su presión de
vapor a esta temperatura es 0.555 mm de Hg. El valor de DAB del naftaleno en
aire a 318 °K es 6.92 x 10-6 m2/s. Calcule la rapidez de evaporación del
naftaleno en la superficie.
Solución: El diagrama de flujo es similar al de la figura 3 (a) . DAB=6.92X10-6
m2/s , pA1 = (0.555/760)(1.01325 x 10 5
) = 74.0 Pa, pA2 = 0, r1 =
2/1000 m, R = 8314 m3. Pa / Kgmol. K, pB1 = P - pA1 = 1.0135 x 105 – 74.0
=1.01251 x 10 5 Pa, pB2= P-pA2 = 1.01325 X 10 5 – 0. Puesto que los valores
de pB1 y pB2 son muy similares,
N A = − D AB dp A
NA = − (35)
πr 2
RT ⎛ pA ⎞
⎜1 − P ⎟dz
⎝ ⎠
en esta geometría, el radio variable r se puede relacionar con la posición z de
la trayectoria como sigue :
⎛ r 2 − r1 ⎞
r = ⎜⎜ ⎟⎟ z + r1 (36)
⎝ z 2 − z1 ⎠
0.0018583 (T )
3/2 1/ 2
⎛ 1 1 ⎞
D AB = ⎜ ⎟
PΓ AB Ω D AB
2 ⎜ M +M ⎟
⎝ A B⎠
donde: T = temperatura del sistema °K
P = presión total del sistema (Atm.)
σi, j= Diámetro de colisión
ΩDi,j = Función adimensional de temperatura y del campo
de potencial intermolecular ( integral de colisión)
= F(KT / ε i-j)
K = constante de Boltzman ( 1.38x10-16 ergios /°K)
Γ A, B =
1
(Γ A + Γ B )
2
ε AB = ε A . ε B
K K K
calcularse .
Ahora bién, si se tiene disponible un valor de Difusividad experimental a una
temperatura T1. , entonces la Difusividad del mismo sistema a una temperatura
T2 se pueden estimar
⎛T2⎞
3/ 2 (Ω D AB )T
(Ω D AB )T
1
D AB ,T 2 = D AB ,T 1 ⎜⎜ ⎟
⎟
⎝ 1⎠
T 2
⎝ M ⎠ prom
donde cprom. Es la concentración total promedio de A + B en Kg mol /m3
M1 = es el peso molecular promedio de la solución en el punto 1 kgmasa / kgmol
ρ1= densidad promedio de la solución en el punto 1 Kg. / m3
La ecuación (1) usa el valor promedio de DAB, que puede variar con la
concentración, y el valor promedio de c, que también puede variar con la
concentración. La ecuación (2) se usa un promedio lineal de c. El caso de
contradifusión equimolar en la ecuación (1) es muy poco frecuente.
D AB c prom. ( − )
NA= ( − ) x x
z 2 z1 x BM A1 A2 (3)
x B 2 − x B1
x BM = ⎛ ⎞
ln⎜⎜ x B 2 ⎟
donde ⎝ x B1⎟⎠ (4)
Nótese que xA1 +xB1= xA2 +xB2 =1.0 . En soluciones diluidas, xBM es cercano a 1.0 y
c es esencialmente constante. Entonces, la ecuación (3) se simplifica a
D AB (c A1 − c A2)
NA=
z 2 − z1 (5)
ρ1 ρ2
+ 972.8 988.1
+
M 1 M2 20.07 18.75 = 3
c prom. = 2
=
2
50 . 6 .Kg .mol / m
para calcular xBM con la ecuación (4) se puede emplear la media lineal, pues xB1 y
xB2 son valores cercano entre sí,
diferencia de los gases, las difusividad DAB no es igual que DBA para líquidos.
líquida. Todos los valores son aplicables a soluciones diluidas del soluto que
Stokes- Einstein, se obtuvo para una molécula esférica muy grande (A) que se
molar.
−16
=
9.96 x10 T
D AB
µ V 1A/ 3
(8)
dode: DAB = Difusividad en m2 / seg.
ING. CARMEN BROCK Página 33 de 35
TRANSFERENCIA DE MASA
DIFUSIÓN MOLECULAR
T = temperatura en Kelvin
µ = viscosidad de la solución Pa. Seg o Kg / mol.seg
VA = volumen molar del soluto a su punto de ebullición
Normalen m3 / Kg mol.
La ecuación (8) no es válida para solutos de volumen molar pequeño. Se ha
intentado obtener otras deducciones teóricas, pero las fórmulas obtenidas no
predicen las difusividades con precisión razonable. Debido a esto, se han
desarrollado diversas expresiones semiteóricas como la correlación de WILKE-
CHANG , puede usarse para la mayoría de los propósitos generales cuando el
soluto (A) está diluido con respecto al disolvente (B).
D AB = 1.173x10
−16
(ϕM B) 1/ 2 T
µ BV 0A.6 (4)
REFERENCIA BIBLIOGRAFICA
∑ (∆τ )
j =n
ψ
j
j =1
= j =n
A
∑R
j =1
j