3.1.

- La transferencia de masa y sus aplicaciones

Las operaciones de transferencia de masa tienen gran importancia desde el punto de vista de las
operaciones y procesos de ingeniería química y otras ciencias afines, y, también en la vida
cotidiana; tal como enjabonarse y lavarse las manos implica la transferencia de masa de materiales
grasosos en la mano hacia el jabón debido a su solubilidad y esta mezcla hacia el agua debido a
su afinidad quimica; remojar un cereal en agua implica la transferencia del agua hacia el interior
del cereal y la transferencia de los componentes solubles del cereal (alcaloides en el caso de los
frijoles) hacia el agua, etc. Dentro del campo de la ingeniería química se utilizan estos principios
de transferencia de masa en operaciones como la destilación, absorción, adsorción,
humidificación, secado, etc. A nivel industrial casi la totalidad de las plantas utilizan las
operaciones indicadas con la finalidad de obtener diversos productos a partir de materias primas
con la inclusión de aditivos. Cabe mencionar que el proceso de transferencia de masa se produce
por dos mecanismos conocidos: La difusión molecular, que es un proceso espontaneo y que se da
por la diferencia de concentraciones y el transporte convectivo que se produce en presencia de
algún tipo de flujo.
3.2.- Conceptos fundamentales de transferencia de masa
3.2.1 La osmosis y la osmosis inversa
Cuando dos fluidos de distinta densidad se encuentran separados por una membrana
semipermeable, se produce una diferencia de presiones entre ambos; el fluido menos denso tiende
a pasar a través de la membrana hasta equilibrar la presion; este fenómeno se llama osmosis.
Esta tendencia natural, en la industria se fuerza a actuar en sentido inverso a lo que haría la
naturaleza; es decir forzar mediante la aplicación de una presion, que el solvente de la solución
concentrada pase al otro compartimento donde la solución es bastante mas diluida; este principio
se denomina osmosis inversa.
Figura 3.1: Esquema del proceso de osmosis y osmosis inversa
3.2.2. La diálisis
La diálisis es una operación en el cual tanto el solvente de una solución diluida (compartimiento
de la derecha) atraviesa una membrana permeable y desde el lado izquierda el soluto de bajo peso
molecular atraviesa la membrana hacia la solución diluida. Este procedimiento es aplicado
durante la depuración de la sangre de los enfermos con insuficiencia renal y recibe el nombre
especifico de hemodiálisis
Figura 3.2. Esquema de un proceso de dialisis

3.2.3-. Difusión molecular

Es el fenómeno por el cual las moléculas disueltas tienden a distribuirse uniformemente en el seno
de un fluido; fenómeno que también puede ocurrir a través de una membrana si es lo
suficientemente permeable, tal como se observa en la Figura 3.3

La difusión como fenómeno de transferencia de masa esta omnipresente en la naturaleza y es de

gran relevancia en la ingeniería, especialmente en aspectos relacionados, con los procesos
químicos, termodinámica, metalurgia, los aspectos ambientales, entre otros.
Uno de los fenómenos de transferencia mas importantes en las áreas de ingeniería, es la
transferencia de masa por difusión molecular, es decir, aquella cuya fuerza impulsora es la
diferencia de potencial químico; sin embargo, el calculo de esta propiedad es bastante dificultosa,
por tanto, se sustituye por la diferencia de concentraciones o su equivalente en diferencia de
presiones o composiciones molares.
Los principios de transferencia de masa por difusión se aplican al transporte de un solo
componente (por ejemplo la difusión de un gas A en el interior de un solido), a una mezcla binaria
(por ejemplo la difusión de cualquier gas en el aire o de un soluto en agua) o de mezclas
multicomponente.
En la figura 3.4 se observa la difusión de un soluto (colorante en agua), que consiste en el
transporte del soluto hacia las diferentes partes del agua, hasta lograr que la concentracion sea
uniforme en todo el recipiente; dependiendo de la naturaleza de los componentes el proceso de
difusión hasta lograr el equilibrio, puede tomar mucho tiempo.

3.2.4: Formas de expresión de la concentración

Existen diversas formas de expresar la concentración de los componentes de una mezcla;
representan en si, las diversas formas de contabilizar la cantidad de una especie (i) en una mezcla
de (n) componentes.
Como regla general, las diversas formas de concentración indican la cantidad de un componente
respecto a una mezcla; en la medida que las cantidades de los componentes y la mezcla son
escritas en diferentes unidades (masa, mol, volumen, etc.), aparecen las diferentes expresiones de
concentración que conocemos y algunas otras menos comunes, pero con aplicaciones muy
importantes en procesos particulares, tal como se presenta en la Tabla 3.1.
Tabla 3.1 .- Formas más comunes de expresar la concentración de un componente en una mezcla.

Definición de los tipos de concentración Formulas

Concentración másica (  i ): Esla cantidad de masa del m
i = i (i)
componente i por unidad de volumen ( v ) v
Concentración de la mezcla(  ): Es la suma de las n

contribuciones de las concentraciones másicas todos los

 = i
i =1
 (ii)
componentes de la mezcla
Fracción másica ( wi ) : Es la relación entre la 
wi = i (iii)
concentración másica del componente i la 
concentración másica de la mezcla (  )
 Suma de fracciones másicas: Esta suma debe ser n

equivalente a la unidad. w
i =1
i =1 (iv)

Concentración molar ( Ci ): Esla cantidad de moles de la ni

Ci = (v)
especie i por unidad de volumen ( v ) v
Concentración de la mezcla ( C ): Es la suma de las n

contribuciones de las concentraciones molares de todos C =  Ci (vi)

i =1
los componentes de la mezcla
Fracción molar ( xi ): Es la relación entre la Ci
xi = (vii)
concentración molar del componente i y la C
concentración molar de la mezcla. ( C )
Fracción molar para el caso de los gases ( yi ): Es la pi
yi = (viii)
relación entre la presion parcial del componente i ( pi ) P
y la presion de la mezcla, ( P )
Relación entre las concentraciones másica y molar i = Ci M i (ix)
La concentración másica total se obtiene mediante:
n

m m i
= T = i =1
(3.1)
v v
La concentración molar total se obtiene mediante:
n

n n i
C= T = i =1
(3.2)
v v
Para el caso de gases, la concentración total está dada por:

P
C= (3.3)
RT

3.2.5.- Relaciones entre las diversas formas de concentración y velocidades medias de

difusión.
Consideremos la difusión de un gas, cuyo transporte se observa a través de un plano imaginario,
tal como se muestra en la Figura 3.5.
Las moléculas de este gas se desplazan a velocidad de diferentes debido a que cada molécula
adquiere una determinada cantidad de energía. Denotamos con Vi a la velocidad de una molécula
especifica respecto a una coordenada fija o estacionaria. La densidad de flujo masico obtenemos
multiplicando la concentración música del componente i por su velocidad media: iVi . Del
mismo modo la densidad de flujo molar obtenemos multiplicando la concentración molar del
componente i por su velocidad media: CiVi . Utilizamos estas expresiones para obtener los
siguientes parámetros.

Velocidad media de masa: ( V ): Este parámetro se obtiene dividiendo la densidad de flujo de

masa total por la suma de las concentraciones másicas de todas las especies.
n

V i i
V= i =1
n
(3.4)

i =1
i

Velocidad media de molar: ( V * ): Este parámetro se obtiene dividiendo la densidad de flujo molar
total por la suma de las concentraciones molares de todas las especies.
n

CV i i
V =
* i =1
n
(3.5)
C
i =1
i

A partir de las ecuaciones anteriores, se puede inferir la diferencia que existe entre la velocidad
de difusión de i, ( Vi ) y la velocidad media de difusión ( V o V * )

Vi − V : es la velocidad de difusión de i respecto a la velocidad media de masa

Vi − V * : es la velocidad de difusión de i respecto a la velocidad media molar

Velocidad media de volumen ( V v ). Con frecuencia, probablemente la velocidad media de

volumen es la velocidad de referencia más importante para obtener datos experimentales, ya que
para un sistema fijo de volumen constante, esta velocidad es cero.
 n
V = Vi CiVi
v
(3.6)
i =1

Vi : Volumen molar parcial de i.

3.2.6. Densidades de flujo masico, molar y volumétrico para un componente
En la Tabla 3.2 se presentan diversas relaciones de densidad de flujo masico, molar y volumétrico
respecto a un punto fijo y relativo a las velocidades media de masa, molar y de volumen
Tabla 3.2 Relaciones de las diversas formas de expresar la densidad de flujo respecto a las formas
de expresión de las concentraciones y velocidades.
Densidades de flujo de De masa: ni = iVi (i)
masa ( ni ) y molar ( N i ) Molar: N i = CiVi ii
relativa a coordenadas
estacionarias
Densidades de flujo de De masa: ji = i (Vi − V ) Iii
masa ( ji ) y molar ( J i ) Molar: J i = Ci (Vi − V ) iv
relativa a la velocidad
media de masa
Densidades de flujo de
*
De masa: ji* = i (Vi − V * ) V
masa ( ji ) y molar ( J vi
i )
relativa a la velocidad
Molar: J i* = Ci (Vi − V * )
media molar
Densidades de flujo de
V V
De masa: jiV = i (Vi − V V ) vii
masa ( j ) y molar ( J ) viii
i
relativa a la velocidad
i
Molar: J iV = Ci (Vi − V V )
media de volumen

3.2.7.- Resumen de las diversas relaciones para sistemas binarios

En las siguientes tablas mostraremos las relaciones útiles cuando se requiere desarrollar
ecuaciones de transferencia de masa para sistemas binarios, las mismas que fueron deducidas de
las ecuaciones anteriores.
Tabla 3.3: Notación de las diversas formas de concentración para sistemas binarios

Defi Densidad de la solución (g/cm3):  =  A +  A (i)

Nicio Concentración de masa de A (g A/cm3 de solución):  A = C A M A
Nes (iii)
Basi Fracción de masa de A: wA =  A /  (iii)
cas Fracción molar de A: x A = C A / C (iv)
Concentración molar de la solución (molg/cm3): C = C A + CB
(v)
Concentración molar de A (molg A/cm3 de solución): C A =  A / M A (vi)
Peso molecular medio de la mezcla: M =  / C (vii)
Defi x A + xB = 1 viii
Nicio ix
M = x A M A + xB M B
Nes
adicio wA / M A x
xA =
nales wA / M A + wB / M B
xi
 dwA
dxA =
M A M B ( wA / M A + wB / M B ) 2

wA + wB = 1 Xii
wA wB 1
+ = Xiii
MA MB M
xA M A Xiv
wA =
x A M A + xB M B
M A M B dxA xv
dwA =
( xA M A + xB M B )2

Tabla 3.4: Notación de velocidades para sistemas binarios

Defi V A : Velocidad de la especie A relativa a coordenadas estacionarias

Nicio
VA − V : Velocidad de difusión de la especie A relativa a V
Nes
Basi VA − V * : Velocidad de difusión de la especie A relativa a V *
cas
V = (1/  )(  AVA +  BVB ) = wAVA + wBVB : Velocidad media de masa
V * = (1/ C )(C AVA + CBVB ) = xAVA + xBVB : Velocidad media molar
Relacio V − V * = wA (VA − V * ) + wB (VB − V * )
Nes adi
Ciona V * − V = xA (VA − V ) + xB (VB − V )
les

Problema 3.1
En dos compartimientos permeables denominados compartimiento I y II, inicialmente se tienen
las cantidades de A y B que se muestran en la Figura 1. Encuentre las cantidades relativas de cada
uno de los componentes A y B en los compartimientos I y II cuando se ha llegado al equilibrio.
Solución
Se observa en una etapa inicial que las masas de A y B en cada compartimiento son diferentes
(100 Kg de A y 200 kg de B); debido a sus pesos moleculares, resulta que el numero de moles en
cada compartimiento y para cada componente es el mismo (2 molKg)
Figura P_3.1: Esquema del proceso de transferencia de masa en dos etapas: Etapa inicial y etapa
de equilibrio

Proceso inicial Proceso en equilibrio

mA = 100 KgA mB = 200 KgB mA = 50 KgA mB = 100 KgB
M A = 50 Kg / molKg M B = 100 Kg / molKg mB = 100 KgB mA = 50 KgA
nA = 2 molKg nB = 2 molKg nA = 1 molKg nB = 1 molKg
xA = 1 xB = 1 nB = 1 molKg nA = 1 molKg
x A = 0.5 xB = 0.5
xB = 0.5 x A = 0.5
Compartimiento I Compartimiento II Compartimiento I Compartimiento II
Cuando se llega al equilibrio el número total de moles en cada compartimiento sigue siendo el
mismo (1 mol Kg de A y 1 mol Kg de B = 2 mol Kg); sin embargo, para que se produzca el
equilibrio hubo un desplazamiento neto de materia de un compartimento a otro (en este caso 50
Kg del compartimiento II al I), lo cual implica que V * = 0 y V  0 .

Problema 3.2
Supongamos que en la Figura 3.2 se tiene una mezcla binaria compuesta de A y B componentes,
para lo cual se dispone de los siguientes datos: x A = 1 / 6 , V * = 12 , VA − V * = 3 , M A = 5M B .
Tomando en consideración estos datos, encuentre:

a) V A , b) VB , c) VB − V * , d) V , e) VA − V , f) VB − V

Figura P_3.2: Distribución de concentración de un sistema de difusión binaria, en un determinado

instante.

Solución

a) A partir de la ecuacion (iii) de la Tabla 3.4 se tiene:

VA − V * = 3 , V * = 12 , con lo cual se obtiene VA = 15 (velocidad de A respecto a

coordenadas estacionarias)

b) A partir de la ecuacion (v) de la Tabla 3.4 se tiene: V = x AVA + xBVB . Ahora

*

1 5 2
reemplazando datos: 12 = 15 +  VB ; de donde VB = 11
6 6 5

c) Dado que ya se disponen de los valores, lo que queda es reemplazar los datos obtenidos
2 3
previamente, con lo cual, VB − V = 11 − 12 = −
*

5 5
 d) A partir de la ecuacion (vi), se tiene: V − V * = wA (VA − V * ) + wB (VB − V * )

Ahora deberemos obtener las fracciones en peso a partir de las fracciones molares y de
los pesos moleculares. A partir de la ecuacion (x) de la Tabla 3.3 se tiene:
xA M A xA (1/ 6) 1
wA = = = =
xA M A + xB M B xA + xB M B / M A (1/ 6) + (5 / 6)(1/ 5) 2
Reemplazando se tiene

1 1 3 1
V − 12 = (15 − 12) + (− ) luego: V = 13
2 2 5 5

1 9 4
e) VA − V = 15 − 13 = = 1 Ojo:
5 5 5
A partir de la Tabla 3.2, se puede obtener las densidades de flujo para sistemas binarios
a partir del conocimiento de los diversos parámetros de concentración y velocidad
Tabla 3.5 Resumen de las densidades de flujo para sistemas binarios

Sistema de referencia Flux molar Flux de masa

Coordenadas fijas N A = C AVA i nA =  AVA ii
Velocidad media molar J A* = C A (VA − V * ) iii j *A =  A (VA − V * ) iv
Velocidad media de masa J A = C A (VA − V ) v j A =  A (VA − V ) vi
Velocidad media de volumen J = C A (VA − V ) vii
V
A
V
j =  A (VA − V V ) viii
V
A

En la Tabla 3.6 se observa más ecuaciones que relacionan las concentraciones de las especies A
y B expresadas en sus diversas formas y son útiles para la aplicación de las ecuaciones de
transferencia de masa.

Tabla 3.6. Otras relaciones de concentración para mezclas binarias

Velocidad de la especie A (cm/s)
VA VA − V VA − V *
nA =  AVA jA =  A (VA − V ) j *A =  A (VA − V * )
Densidad de flujo de molar de la especie A (molg/cm2s)
N A = C AVA j A = C A (VA − V ) J A* = C A (VA − V * )
Suma de las densidades de flujo de masa (g/cm2s)
nA + nB = V j A + jB = 0 j *A + jB* =  (V − V * )
Suma de las densidades de flujo de molar (molg/cm2s)
N A + N B = CV * J A + J B = C (V * − V ) J A* + J B* = 0
 Densidad de flujo en función a n A y nB
N A = nA / M A j A = n A − w A ( n A + nB ) j*A = nA − xA (nA +
MA
nB )
MB
Densidad de flujo en función a NA y
nA = N A M A J A = N A − wA ( N A +
MB
NB ) J A* = N A − xA ( N A + N B )
MA
Densidad de flujo en función a JA y V
nA = j A +  AV J A = jA / M A j *A = ( M / M B ) j A
de flujo en función a J A* y V*
N A = J A* + C AV * J A = ( M B / M ) J A* j *A = J A M A
Nota: La fórmula de color azul es una de las más utilizadas
Problema 3.3:
Relaciones de densidades de flujo (I)

Basándose en la definición de Vi (Volumen molar de i) y considerando V C i i = 1 pruebe que

V i
V
Ci = 0

Solución

Partimos de la ecucion (iii) de la Tabla 3.2

J iV = Ci (Vi − V V ) (1)

Multiplicando ambos miembros de (1) por Vi

Vi J iV = Ci Vi (Vi − V V ) (2)

Aplicando sumatoria a ambos lados de la igualdad

n n

 JiV Vi =  Ci Vi (Vi − V V )
i =1 i =1
(3)

De la ecuacion (3.6) se tiene

n
V v = Vi CiVi (4)
i =1

Reemplazando (4) en (3)

n n n

 JiV Vi = V v − V v  Ci Vi = V v (1 −  Ci Vi )
i =1 i =1 i =1
(5)

Dado que V C i i =1
Reemplazando en (5) se obtiene
n

J
i =1
V
i Vi == V v (1 − 1) = 0 (6)

Problema 3.4

Relaciones de densidad de flujo (II)

n
Demuestre que se cumple la relación:  i =1
J i* = 0

Solución

A partir de la ecuacion (vi) de la Tabla 3.2 se tiene:

J i* = Ci (Vi − V * ) (1)

Aplicando sumatoria a ambos miembros de la ecuacion (1)

n n n n

 Ji* =  Ci (Vi − V * ) =  CiVi − CiV *

i =1 i =1 i =1 i =1
(2)

La ecuacion (2) puede ser modificada del siguiente modo:

n n

n CV C i i i n

J *
i = i =1
n
i =1
−  CiV * (3)
i =1
C
i =1
i
i =1

A partir de (3) se obtiene

n n n

J
i =1
*
i =V *
C − C V
i =1
i
i =1
i
*
=0 (4)

Problema 3.5

Relaciones de densidad de flujo (III)

Demuestre que se cumple la siguiente igualdad

 2 DAB
J A* = − wA
CM A M B

Solución

La fracción molar del componente A viene dada por:

wA / M A
xA = (1)
xA / M A + xB / M B

Diferenciando la ecuacion (1)

 dwA
dxA = (2)
M A M B ( wA / M A + wB / M B ) 2
Aplicando operadores a cada miembro de la ecuacion (7)

wA
xA = (3)
M A M B ( wA / M A + wB / M B ) 2
Utilizando la ley de Fick

J A* = − DABC A = −CDABx A (4)

Nota: Para esta demostración, se requiere el conocimiento de la ley de Fick dada por la ecuacion
(4), que por cierto aun no se ha definido; sin embargo, a manera de adelanto en base a la
ecuacion (8.141) ( = − ) del anexo, la ley general del transporte difusivo adecuarse para
la demostración de este ejemplo.

Reemplazando (3) en (4)

CDAB wA CDAB wA

J A* = − =− (5)
M A M B ( wA / M A + wB / M B ) 2
 1  1 2
M AM B ( A  + B )
 MA  MB

CDAB wA  2CDAB wA

J A* = − =− (6)
C C
M A M B ( A + B )2 M A M BC 2
 
Finalmente se obtiene

 2 DAB
J A* = − wA (7)
CM A M B

Problemas propuestos

Demuestre las siguientes identidades

C 2 M A M B DAB
1.- j A = − xA


2.- N A = −CDABxA + xA ( N A + N B )

3.- nA = −  DABwA + wA (nA + nB )

MA
4.- J A = VB
V
jA
MB

 xB 
5.- J VA + J BV =  N A − N B  (CVB − 1)
 xA 
 n
xi
6.- 
i =1
J iV
Ci
=0

CDAB
7.- C (VA + VB ) = − xA
xA xB

C A VA
8.- J VA = 1 −
M A wB

Ojo todo lo que es volumen molar con color rosado se debe cambiar a una letra minúscula

3.3.- Ecuacion general del transporte de masa molecular

La ecuacion del transporte de masa por difusión molecular, conocida como la velocidad de
transferencia de masa expresada como la densidad de flujo de molar o flux molar, viene dada
por:

(3.7)

El coeficiente de transferencia de masa es un factor que multiplica a la fuerza impulsora, y,

además, depende de las propiedades físicas (p. ej. la difusividad) de la geometría (longitud,
radio, etc.) y las condiciones de operación (presion, temperatura, etc.)

La difusividad mide la movilidad de las moléculas y puede ser obtenida experimentalmente o

usando correlaciones para el caso de los gases y líquidos.

La fuerza impulsora es el gradiente de concentración molar o su equivalente expresado en

gradiente de presiones parciales o gradiente de composiciones molares .

3.4.- La primera ley de Fick para difusión molecular

El transporte de masa molecular se define como la transferencia de moleculas individuales por
medio de desplazamientos desordenados de las mismas. Se supone que cada molécula se
desplaza en línea recta y cambia su dirección al chocar con otras moleculas, cuyo movimiento
es al azar, por lo que se conoce también como trayectoria aleatoria, tal como se observa en la
Figura 3.6
En la figura 3.6 (a), se observa que las moléculas se desplazan en forma desordenada de la
región 1 a la región 2 chocando elásticamente ya sea con las moléculas del mismo tipo o con las
de otro tipo. Si hay un mayor número de moléculas del tipo A en la región 1, entonces habrá
mayor número de moléculas desplazándose de la región 1 a la región 2; por tanto, habrá un
movimiento neto de las moléculas del tipo A de izquierda a derecha, lo cual permite concluir
que el proceso de transferencia de masa se da desde una región de mayor concentración a una
región de menor concentración.

En la Figura 3.6 (b) y (c) se observa el mismo hecho del fenómeno de transferencia de masa
descrito de mayor a menor concentración, hasta alcanzar el equilibrio

Ahora consideremos el esquema de la Figura 3.7, donde un soluto A se esta difundiendo en un

material solido o un material liquido (comprendido entre dos placas planas y paralelas), en la
dirección z

En la parte superior se mantiene su concentración molar a C A = C Ah  0 , mientras que en la

parte inferior es mantenida a CA = 0 (en régimen estacionario); cuando el sistema es no
reaccionante, se puede observar el desarrollo del perfil lineal de concentración del soluto A, tal
como se muestra en la Figura 3.7.
El flujo molar de A, ( n A ) transferido por unidad de area desde el plano superior, es proporcional
al gradiente de concentracion.

nA C Ah
 (3.8)
A h
nA C A
 (3.9)
A z
El segundo termino de (2) se puede expresar en forma diferencial, y para convertir la constante
de proporcionalidad en una igualdad, se introduce una constante de proporcionalidad, que en
este caso se denomina coeficiente de difusión, D ; el primer termino se denomina densidad de
*
flujo molar, ( J Az ), denominado densidad de flujo molar o flux molar relativo a la velocidad
media molar V * y para efectos prácticos, el carácter (*) será obviado, salvo sea necesario su uso
y por evitar errores de interpretación de las respectivas densidades de flujo.

La ecuacion (2) escrita en forma diferencial, conduce a:

dC A
J Az = − D (3.10)
dz
La ecuacion (9), se conoce como la primera ley de Fick

J Az : Densidad de flujo molar (moles/supxtiempo); se ha omitido (*) en el superíndice.

C A : concentración molar del componente A (mol/vol)

D : Coeficiente de difusión del componente A (sup/tiempo)
z : Dirección del movimiento molecular del fluido (long)

Nota: El signo negativo aparece debido a que el transporte de masa se da siempre en el sentido
decreciente.

La ecuacion (3.10) puede ser escrita en términos de otras propiedades, tales como

dy A
J Az = −CD (3.11)
dz
D dp A
J Az = − (3.12)
RT dz
dwA
J Az = −  D (3.13)
dz
La nomenclatura de las ecuaciones (3.11), (3.12) y (3.13) son las mismas que se ha descrito en
los ítems anteriores.

Se puede definir una ecuacion general para la ley de Fick, la cual fue derivada por Hooyman y
otros con la finalidad de transformar las ecuaciones de densidad de flujo de un sistema de
referencia a otro, la misma que está dada por:
 1− X A 1
J VA = − DAB  xA (3.14)
1 − xA V

V : Volumen molar de la mezcla

x A : Fracción molar

X A : Factor de peso relacionado con la velocidad de referencia apropiado

J VA : Densidad de flujo molar relativo a la velocidad media volumétrica.

Usando la ecuacion anterior, las distintas formas de la Ley de Fick pueden ser deducidas
fácilmente. Los factores de peso respectivo son:

Tabla 3.7 Factores de peso para la ecuacion (14)

Ley de Fick Velocidad de referencia X A (Factor de peso)

jA V wA
J A* V* xA
J VA VV C A VA

Problema 3.6. Relación de densidades de flujo (IV)

Usando la ecuacion (3.14) encuentre las siguientes densidades de flujo

a) La densidad de flujo molar respecto a la velocidad media molar

b) La densidad de flujo de masa respecto a la velocidad media de masa.

Solución

a) De la ecuacion (3.14) se tiene

1− X A 1
J A* = − DAB  xA (1)
1 − xA V

De la Tabla ccc X A = xA y C=
1
V

Reemplazando en (1)

1 − xA
J A* = − DAB  CxA = −CDABxA (2)
1 − xA
La ecuacion (2) representa una de las formas de la ley de Fick para transporte difusional

b) De la ecuacion (3.14) se tiene

1− X A 1
J A = − DAB  xA (3)
1 − xA V
De la Tabla ccc X A = wA y C= 1
V

Reemplazando en (3)

1 − wA 1
J A = − DAB  xA (4)
1 − xA V
Por otro lado del problema ddddd

dwA dwA
dxA = = (5)
M A M B ( wA / M A + wB / M B ) 2
M A M B (1/ M ) 2
Reemplazando (5) en (4)

wB CM 2
J A = − DAB  wA (6)
xB M A M B

Por otro lado: wB = ( M B / M ) xB y  = CM , jA = J A M A

jA = −  DABwA (7)

3.5 Mecanismos de transferencia de masa

Los mecanismos de transferencia de masa conocidos son el difusivo y el convectivo

El transporte de masa difusivo se produce únicamente por diferencia de concentraciones

expresado en términos de concentración o fracción molares o presiones parciales, la misma que
es aplicable a un solo tipo de molécula o a una mezcla. Por otro lado, estudios realizados en un
mismo componente estando a las mismas condiciones (presion, temperatura, concentración),
se observó que las moléculas de este mismo componente podían difundirse, sin que se produzca
transferencia de materia; esto ocurre porque el verdadero potencial de transferencia de materia
es el potencial químico, y debido a su dificultad en el cálculo, para efectos prácticos se considera
la diferencia de concentración expresada en su diferentes formas como el potencial de
transferencia de masa. En el caso de la difusión de un componente en si mismo, el coeficiente
de difusión se denomina coeficiente de autodifusion.

Son ejemplos de transporte de masa por difusión molecular: Transporte de masa de un fluido
en el interior de un sólido poroso, el transporte de masa de un soluto en un solvente sin
movimiento o en flujo laminar; el transporte de masa en una película estancada, etc.

El transporte de masa convectivo se da debido al movimiento global del fluido o la acción de un

mecanismo externo que propicie el movimiento del fluido conteniendo el componente
especifico que se esta estudiando. Estos mecanismos pueden ser una corriente de fluido, una
fuente de calor o un sistema de agitación que propicie la turbulencia. El mecanismo de
transferencia de materia por convección se produce en la interfase fluido-fluido o solido fluido;
este mecanismo de transferencia de masa en cualquier caso, es el más importante. Son ejemplos
de transporte de masa convectivo, la transferencia de masa desde el interior de una partícula
solida hacia una corriente de fluido, el transporte de masa de un soluto en un flujo turbulento,
el transporte de masa de un soluto en un recipiente mediante una agitación mecánica.

Algunos ejemplos de transporte de masa difusivo y convectivo son mostrados en la figura xxx

3.6.- Estimación de los coeficientes de difusión de gases a baja presion

3.6.1: Método de Chapman y Enskog

Chapman y Enskog trabajando en forma independiente y utilizando la teoría cinetico molecular

de los gases para relacionar las propiedades de los gases con las fuerzas que actúan entre las
moléculas, obtuvieron la ecuacion para el calculo del coeficiente de difusion de un gas A en otro
gas B.

1/2
1.858 10−27 T 3/2  1 1  (3.15)
DAB =  + 
P AB
2
D  MA MB 

 AB = ( A +  B ) / 2 (3.16)

 AB  
1/2

= A  B  (3.17)
k  k k 

DAB : Coeficiente de difusión, m2/s

T : Temperatura absoluta, K
P : Presion absoluta, atm.

 AB : Diámetro promedio de colisión, m.

D : Integral de colisión, que depende de  AB / k y de T

 : energía potencial mínima donde las fuerzas de atracción y repusion entre moléculas se
encuentran en equilibrio, J
k : Constante de Boltzman, J/K

 AB / k : Parametro energético, K

Los parámetros  / k y D son conocidos como parámetros de Lennard-Jones y el diámetro de

colisión, se encuentran en tablas; sin embargo, la integral de colisión puede ser calculada
mediante la ecuacion (3.18)

Para gases no polares

1.06036 0.1930 1.03587 1.76474

D = + 0.47635T * + 1.52956T * + 3.89411T * (3.18)
T *0.1561 e e e

T * = T / ( AB / k ) (3.19)

Para gases polares

D = D + 0.19 AB
2
/T* (3.20)

 AB = ( A   B )1/2 (3.21)

1/3
 1.585  veb 
 = 2  (3.22)
 1 + 1.3 

 p2
 = 1.94 10 3
(3.23)
vebTeb

( / k ) = 1.18(1 + 1.3 2 )Teb (3.24)

Teb : Temperatura de ebullición, K

 p : Momento dipolar, Debyes

veb : Volumen molar en estado condensado a su temperatura de ebullición, cm3/molg

 : Parámetro adimensional
 / k : parámetro energético, K
Cuando los diámetros de colisión y las constantes energéticas no se encuentran disponibles,
dichos valores se pueden estimar a partir de:

 = v1/3
C
 = 11.8veb1/3  / k = 0.75Tc  / k = 1.21Teb

vc : Volumen molar crítico, m3/kmol

veb : Volumen molar crítico, m3/kmol

Teb : Temperatura de ebullición normal, K

Tc : Temperatura critica del gas , K

−10
 : Diámetro de colisión, 10 m

Los coeficientes de difusion calculados a determinadas condiciones de temperatura y presion

pueden ser usados para predecir su valor a otras condiciones mediante la relación:

k (T1 )
3/2
P1  T2 
DAB ,T2 , P2 = DAB ,T1 , P1   (3.25)
P2  T1  k (T2 )

Problema 3.6

Coeficiente de difusión binaria de una mezcla a dos condiciones de temperatura

Estime el coeficiente de difusión del sistema hidrogeno-amoniaco a 25oC y 1 atm.; con este valor
obtenga el coeficiente de difusión de dicha mezcla a 85oC y 1 atm.

Solución

Identificamos a los componentes de la mezcla como, hidrogeno: A; amoniaco: B

De la Tabla xxxx obtenemos los siguientes parámetros.

 A = 2.877 10−10 m ( / k ) A = 59.7 K (1)

 B = 2.900 10−10 m ( / k ) B = 558.3K (2)

Ahora es posible calcular el parámetro energético de la mezcla

 AB  
1/2

= A  B  = (59.7  558.3)1/2 = 182.57 K (3)

k  k k 
La temperatura adimensional se obtiene mediante:

T1* = T1 ( K ) / ( AB / k ) = (25 + 273) /182.57 = 1.63225 (4)

Reemplazando en (3.18)

1.06036 0.1930 1.03587 1.76474

D = + + + = 1.158 (5)
T *1.63225 e1.63225T * e1.63225T * e1.63225T *
El diámetro de colisión se obtiene mediante un promedio aritmético, ecuacion (3.16)

 AB = ( A +  B ) / 2 = (2.877 + 2.900) / 2 10−10 m (6)

Estos datos obtenidos son reemplazados en (3.15)

 1/2
1.858 10−27 T 3/2  1 1  1.858 10 −27 2983/2
1/2
1 1 
DAB =  +  =  + 
P AB D
2
 MA MB 
−10 2
1 (2.8635 10 ) 1.158  2 17 

DAB = 7.52 10−5 m 2 / s (7.83 10−5 m 2 / s) (6)

Ahora, cuando la temperatura es 85oC

La temperatura adimensional se obtiene mediante:

T2* = T2 ( K ) / ( AB / k ) = (85 + 273) /182.57 = 1.961 (7)

La integral de colisión se obtiene mediante:

1.06036 0.1930 1.03587 1.76474

D = + 1.961T * + 1.961T * + 1.961T * = 1.082 (8)
T *1.961 e e e
Estos nuevos valores son reemplazados en la ecuacion (zzzz)

k (T1 )
3/2
 T2 
3/2
P1  1  358   1.158 
−5
DAB ,T2 , P2 = DAB ,T1 , P1   = 7.52  10     
P2  T1  k (T2 )  1  298   1.082 

DAB = 1.06 10−4 m2 / s (1.093 10−4 m2 / s) (9)

Nota: Los valores presentados entre paréntesis son los valores experimentales para cada una de
las temperaturas dadas.

3.6. 2 Métodos empíricos

Método de Fuller, Schettler y Giddings

Propusieron la siguiente ecuacion empírica para el calculo del coeficiente de difusion binaria de
una mezcla, de gases.

110−3  T 1.75 (1/ M A + 1/ M B )1/2

DAB = (3.26)
P ( vA )1/3 + ( vB )1/3 

DAB : Coeficiente de difusión, m2/s

T : Temperatura absoluta, K
P : Presion absoluta, atm.

 v : Es la suma de los volúmenes atómicos de difusion mostrados en la Tabla 8

Tabla 3.8.- Volúmenes atómicos y moleculares de difusion para la Ecuacion (3.26)

Volúmenes atómicos y estructuras parciales

C 16.5 (Cl) 19.5
H 1.98 (S) 17.0
 O 5.48 Anillo aromático -20.2
(N) 5.69 Anillo heterociclo -20.2
Volúmenes moleculares de difusion de moléculas sencillas v
H2 7.07 CO 18.9
O2 16.6 CO2 26.9
He 2.88 N 2O 35.9
N2 17.9 NH3 14.9
D2 6.70 H2O 12.7
Aire 20.1 (CCl2F2) 114.8
Ar 16.1 (SF6) 64.7
Kr 22.8 (Cl2) 37.7
(Xe) 37.9 (Br2) 67.2
- -
Ojo revisar por que algunas formulas estan entre paréntesis y también revisar el contenido de
la ultima fila.

Problema 3.6

Difusividad de una mezcla binaria por el método de Fuller et. al.

Estime el coeficiente de difusión de la mezcla metano-dióxido de carbono a OoC y 1 atm- El valor

experimental a estas condiciones es 0.153 cm 2/s.

Solución

Para el calculo correspondiente hacemos uso de la ecuacion (3.26)

1 10−3  T 1.75 (1/ M A + 1/ M B )1/2

DAB = (1)
P  ( v A )1/3 + ( vB )1/3 
2

Ahora con la ayuda de la Tabla 3.8 obtenemos los volúmenes moleculares del metano y del
dióxido de carbono

( v)CH 4 = C + 4  H = 16.5 + 4 1.98 = 24.42

( v)CO2 = 26.9

Reemplazando en (1)

1 10−3  T 1.75 (1/ M A + 1/ M B )1/2 110−3  273.151.75 (1/16 + 1/ 44)1/2

DAB = = = 0.154cm 2 / s
P ( v A ) + ( vB ) 
1/3 1/3 2
1  24.2 + 26.9 
1/3 1/3 2

3.7. Estimación de los coeficientes de difusión de líquidos

Método de Wilke-Chang

Al considerarse el proceso de difusión como el transporte de moléculas esféricas en soluciones

diluidas, teóricamente se llega a la ecuacion conocida como la ecuacion de Stokes-Einstein

R T kT
DAB = = B (3.26)
N 6 B rA 6 B rA
Basandose en esta expresión de Sokes-Einstein, Wilke y Chang propusieron la siguiente
modificación empírica de la ecuacion (3.26)

7.4 10−8 ( M B )1/2 T

0
DAB = (3.27)
 B vA0.6

DAB : Coeficiente de difusion binaria a cualquier concentracion de soluto y solvente m2/s

0
DAB : Coeficiente de difusion de una solución de un soluto A en un disolvente B, m2/s

B : Viscosidad del disolvente, cP.

rA : Radio de las moléculas, m

k B : Constante Boltzman, J/K

R : Constante univeral de los gases, J/molg K
N : Numero de Avogadro, moléculas/molg

v A : Volumen molar del soluto a su temperatura normal de ebullición, cm3/molg

M B : Masa molecular del disolvente, g/molg.

 : Factor de asociación del disolvente , cuyos valores son 2.6 para el agua, 1.9 para el metanol,
1.5 para el etanol y 1.0 para disolventes no asociados.

El volumen molar en su punto de ebullición normal puede ser calculado usando una ecuacion
de estado. Con este conocimiento, se puede hacer uso directamente de la ecuacion (3.27) o del
diagrama de la Figura 3.8

Figura 3.8: parámetros para el cálculo del coeficiente de difusión de mezclas liquidas diluidas en
función del volumen molar del soluto.
Método de Scheibel

Scheibel modifico el método de Wilke-Chang para el cálculo del coeficiente de difusión binaria
de mezclas liquidas, eliminando el factor de asociación,  , con lo cual obtuvo la siguiente
correlación:

kT
0
DAB = (3.28)
 B v1/3
A

  3v  2/3 
−8
k = 8.2  10 1 +  B   (3.29)
  v A  

Si el disolvente no es agua y v A  vB , el valor de k se puede suponer igual a 25.2 10−8 .Si el

disolvente es el benceno y v A  2vB el valor de k se puede suponer igual a 18.9 10−8 . Para
otros disolventes si vA  2.5vB , el valor de k se puede suponer igual a 17.5 10−8 . Este método
es recomendable sobre todo para mezclas de sustancias orgánicas.

Método de Hayduk-Laudie

La ecuacion propuesta por estos investigadores, es la siguiente:

DA0 , H 2O = 13.26 10−5  H−1.4

2O A
v −0.589 (3.30)

DA0 , H 2O : Difusividad del soluto A en el agua, cm2/s

 H 2O : Viscosidad del agua, cP.

v A : Volumen molar del soluto A, cm3/molg

Ojo en todos los casos falta escribir adecuadamente el volumen molar

Dependencia del coeficiente de difusión con la concentración

Vignes sugirió la siguiente correlación para expresar la variación del coeficiente de difusión en
función de la concentración

DAB = ( DAB ) (D )
0 xB 0 xA
BA
(3.31)

DAB : Es el coeficiente de difusión a cualquier concentración obtenida a partir de la

concentración a dilución infinita de A en B y de B en A. Esta ecuacion no puede ser utilizada para
mezclas con una asociación fuerte.

Problema 3.7

Difusividad del ácido acético en agua

Estime el coeficiente de difusión de ácido acético diluido en agua a 12.5oC cuyo valor
experimental es 0.95x10-5 cm2/s. Cual seria dicho valor del coeficiente de difusión si la
concentración fuese del 20% de ácido acético?

Datos

Densidad del ácido acético a su temperatura normal de ebullición: 0.973 g/cm3

Densidad del agua a su temperatura normal de ebullición: 0.958 g/cm3

Viscosidad del ácido acético a 12.5 oC: 1.4 cp.

Viscosidad del agua a 12.5 oC: 1.2 cp.

Solución

La nomenclatura para usar es la siguiente. Ácido acético (A), agua (B)

a) Hacemos uso de la ecuacion de Wilke-Chang

7.4 10−8 ( M B )1/2 T

0
DAB = (1)
 B vA0.6
Ahora calculamos el volumen molar de A

1 1
vA = MA =  60 = 64.034 cm / mo lg
A 0.937

El factor de asociación del solvente (agua) es 2.6

Reemplazando datos en (1)

7.4 10−8 (2.6 18)1/2 (12.5 + 273)

0
DAB = = 9.93 10−6 cm2 / s
1.2  64.034 0.6

Respecto al valor experimental, hay una desviación del 9.1%

b) Ahora utilizamos el método de Scheibel.

 kT
0
DAB = (2)
 B v1/3
A

  3v  2/3 
−8
k = 8.2  10 1 +  B   (3)
  v A  

Ahora calculamos el volumen molar de B

1 1
vB =  MB = 18 = 18.79 cm / mo lg
B 0.958
Reemplazando en (2)

  3 18.79 2/3 
−8 −8
k = 8.2 10 1 +    = 15.73 10
  64.034  
Ahora reemplazando en (1)

kT 15.73 10−8 (12.5 + 273)

0
DAB = = = 0.935 cm2 / s
 B vA
1/3
1.2  64.034 2/3

La desviación respecto al valor experimental es de 2.8%

c) Método Hayduk-Laudie

Hacemos uso de la ecuacion (3.30)

DA0 , H 2O = 13.26 10−5  H−1.4

2O
 vA−0.589 (4)

Reemplazando datos en la ecuacion (4)

DA0 , H 2O = 13.26 10−51.2−1.4  64.034−0.589 = 0.887  10−5 cm2 / s

Ahora, la desviación respecto al valor experimental es -2.5%

El coeficiente de difusión a cualquier concentración se obtiene mediante la ecuacion (3.31)

DAB = ( DAB ) (D )
0 xB 0 xA
BA
(5)

El coeficiente de difusión de B en A calculamos por el método de Scheibel

  3v 2/3    3  64.034 2/3 

−8
8.2 10 1 +  A
  T 8.2 10 1 + 
−8
  (12.5 + 273)

  vB  
 = 
  18.79  
0
DBA =
 A vB
1/3
1.4 18.791/3

0
DBA = 3.59 10−5 cm2 / s
Finalmente reemplazando en (5)
DAB = ( DAB ) (D ) = ( 0.9375 10−5 ) ( 3.59 10 )
xB xA 0.8 −5 0.2
0 0
BA = 1.224 10−5 cm2 / s