LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

COEFICIENTE DE DIFUSIVIDAD

Angel Ariztizabal Pedro Ballen

Facultad de Ingeniera Departamento de Ingeniera Qumica y Ambiental.

1. OBJETIVOS
1.1 General

Determinar el coeficiente de difusin de un vapor a travs de fase gaseosa.

1.2 Especficos

Comparar el resultado experimental contra los resultados obtenidos mediante correlaciones

semiempiricas y con los hallados en la literatura para el coeficiente de difusin del sistema
acetona-aire.

2. MARCO TEORICO
Difusin molecular: es el viaje de uno o ms componentes a travs de otros ocasionados por un gradiente
de concentraciones o de potencial qumico cuando se ponen en contacto dos fases inmiscibles, que se
encuentran estancadas o en rgimen laminar. La rapidez con la cual se transfiere un componente en una
mezcla depender del gradiente de concentracin existente en un punto y en una direccin dados. Su
movimiento est descrito por el flux, el cual est relacionado con la difusividad por medio de la Primera
Ley de Fick para un sistema isobrico e isotrmico.
Primera Ley de Fick: las leyes de transferencia de masa, muestran la relacin entre el flujo de sustancia
que se difunde y el gradiente de concentracin responsable de dicha transferencia. La relacin bsica para
difusin molecular define el flux molar relativo a la velocidad molar promedio, el cual se designa por J A.
Una relacin emprica para este flux molar, postulada por Fick, define la difusin del componente A en un
sistema isobrico o isotrmico, as:

J A =D AB

dCA
dz

(1)

Para difusin en la direccin z, donde J A es el flux molar en la direccin z relativa a la velocidad molar
promedio, dCA/dz es el gradiente de concentracin de A en la direccin z y DAB es el coeficiente de
difusin molecular que se define como el factor de proporcionalidad del componente A en el componente
B de la ecuacin de Fick; el signo negativo hace nfasis en el hecho que la difusin ocurre en el sentido de
la disminucin de la concentracin.
Una relacin de flujo ms general, que no es restringida a sistemas isobricos e isotrmicos, para la Ley de
Fick, es la siguiente:

J Az=c D AB

dy A
dz

(2)

LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

Donde c es la concentracin global de la mezcla, mientras que y A es la fraccin molar de A.

Difusin de A a travs de B estancado
Estado estable. El coeficiente de difusin para un sistema gaseoso, puede ser medido experimentalmente
en una celda de difusin. Consta de un tubo angosto parcialmente lleno con lquido puro A, el cual se
mantiene a temperatura y presin constante por medio de un bao de agua. Un gas B se hace fluir a travs
del extremo abierto del tubo, este debe tener una solubilidad despreciable en el lquido A al tiempo que
debe ser inerte qumicamente a l. El componente A se vaporiza y difunde dentro de la fase gaseosa. La
velocidad de vaporizacin de A, puede ser expresada matemticamente en trminos del flujo msico o
molar.
Estado seudo-estable. En muchas operaciones de transferencia de masa, uno de los lmites puede
moverse con el tiempo, tal es el caso real de la celda de difusin de Arnold. Si la longitud de la trayectoria
de difusin vara muy poco en un tiempo relativamente largo, un modelo de difusin seudo-estable se
puede usar.

Figura 1. Esquema de difusin en estado seudo estable

Realizando un balance de materia a travs de un volumen diferencial a lo largo del tubo, teniendo a A t
como rea transversal al flujo se obtiene:

N Az A t ]z + z N Az A t ]z=0

(3)

Dividiendo esta expresin por el volumen AtZ y evaluando el lmite cuando Z tiende a cero, se obtiene
la ecuacin diferencial siguiente:

d N Az
=0
dz

(4)

Esta relacin refleja un flujo molar constante de A en la fase gaseosa desde z 1 hasta z2. Ahora,
considerando solamente el plano que pasa por z 1 y la restriccin que el gas es insoluble en el lquido A, se
deduce que el flujo NBZ es igual a cero a lo largo de toda la columna, de tal manera que el componente B
se comporta como un gas estancado. El flujo molar constante de A se puede describir por la siguiente
ecuacin:

N Az=D AB c

d yA C A
+
( N Az + N Bz )
dz
c

(5)

LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

Reemplazando

CA
= y A , se llega a:
c

N Az=

D AB c d y A
1 y A dz

(6)

Esta ecuacin se puede integrar en las condiciones lmites:

Z = Z1 yA = yA1;

Z = Z2 yA = yA2

Integrando se obtiene:

N Az=D AB

1 y A 2
c
ln
z 2z 1
1 y A 1

(7)

La media logartmica de la concentracin del componente B se define como:

y A 1 y A 2

y BML =

ln

1 y A 2
1 y A 1

(8)

Si se reemplaza (8) en (7) se obtiene:

N Az=

c D AB y A 1 y A 2
z 2z 1 y BML

(9)

La ecuacin (9) describe el flujo de masa en la pelcula de gas estancado. Colocando el flux de A como
funcin del volumen y la concentracin de A en el lquido se tiene:

N Az=

C AL dV
A t dt

(10)

Sabiendo que V = At * h:

dV
dh
=A t
dt
dt

(11)

y que dh = -dz, se obtiene:

N Az=C AL

dz
dt

(12)

Igualando (12) y (9) e integrando se obtiene:

z 2t z 2
1 D AB P ( y A 1 y A 2)
=
t
2
C AL RT
y BML

(13)

LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

Donde: Zt = Z2 Z1,t

Zt0 = Z2 Z1,t =0

Esta ecuacin corresponde a una lnea recta, donde la pendiente es

1 D AB P y A 1 y A 2
C AL RT
y BML

(14)

A partir de esta ecuacin y teniendo los datos de altura para diferentes tiempos, se puede hallar la
difusividad de un sistema binario empleando la Celda de Arnold.
Ecuaciones semiempricas para el clculo de coeficientes de difusividad
Cuando no se cuenta con datos experimentales de los coeficientes de difusividad para distintos sistemas,
se puede recurrir a expresiones basadas en su mayora en la teora cintica de los gases, para el clculo de
la difusividad.
Ecuacin de Hirschfelder- Bird- Spotz
DAB

104 1.084 0.249

1/ M A 1/ M B T 3 / 2
PT

rAB

1/ M A 1/ M B

f kT / AB

(15)

Siendo:
DAB = difusividad, m2 / s
T = temperatura absoluta, K.
MA, MB = peso molecular de A y B, kg / kmol.
PT = presin absoluta, N / m2.
rAB = dimetro de colisin efectivo, nm.
= (rA + rB) / 2.
AB = constante de fuerza de Lennard-Jones.
=(A B)1/2
k = constante de Boltzmann.
f (kT/AB) = D = integral de colisin.
Los valores de r y son caractersticos para cada molcula, estos se encuentran tabulados para un nmero
considerable de compuestos o pueden ser calculados partiendo de contribucin de tomos y la temperatura
normal de ebullicin, respectivamente.

Ecuacin de Wilke Lee

DAB

Siendo:

3.03 0.98 / M AB1 / 2

PT M AB1 / 2 rAB

10 T
3

f kT / AB

3/ 2

(16)

LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

DAB = difusividad, cm2 / s

T = temperatura absoluta, K.
MA, MB = peso molecular de A y B, kg / kmol.
MAB = 2[(1 / MA) + (1 / MB)]-1
PT = presin absoluta, N / m2.
rAB = dimetro de colisin efectivo, nm.
= (rA + rB) / 2.
AB = constante de fuerza de Lennard-Jones.
=(A B)1/2
k = constante de Boltzmann.
f (kT/AB) = D = integral de colisin.
De igual manera que el anterior caso, los valores de r y son caractersticos para cada molcula.
Ecuacin de Fuller, Shettler y Giddings

DAB

0.00143 T 1.75

PT M AB1/ 2

1/ 3
A

1/ 3
v B

(17)
Siendo:
DAB = difusividad, cm2 / s
T = temperatura absoluta, K.
MA, MB = peso molecular de A y B, kg / kmol.
MAB = 2[(1 / MA) + (1 / MB)]-1
PT = presin absoluta, N / m2.
v = se calcula para cada compuesto sumando los volmenes atmicos de difusin
tabulados en la referencia.

Coeficientes de difusin para sistemas gaseosos binarios a bajas presiones: prediccin a partir de la
teora
El desarrollo hecho para la construccin de esta ecuacin se basa en la teora cintica de los gases, para
hacer un estudio detallado de como ocurre el proceso difusivo en fase gaseosa, teniendo en cuenta la
trayectoria libre media de las molculas. En este tambin consideran los parmetros de la integral de
colisin y la constante de fuerza de Lennard Jones para la mezcla.
DAB

0.001858 T 3 / 2 M A M B / M A M B

2
PT rAB f kT / AB

Siendo
DAB = difusividad, cm2 / s
T = temperatura absoluta, K.
MA, MB = peso molecular de A y B, kg / kmol.
PT = presin absoluta, atm.
rAB = dimetro de colisin efectivo, A.

1/ 2

(18)

LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

= (rA + rB) / 2.
AB = constante de fuerza de Lennard-Jones.
=(A B)1/2
k = constante de Boltzmann.
f (kT/AB) = D = integral de colisin

Las estimaciones de DAB, basadas en la anterior ecuacin se conocen como valores tericos, aunque estn
implicados en el clculo empirismos evidentes. Adems, sobre la base de est ecuacin, es que se han
planteado diferentes modelos de clculo de coeficientes de difusividad. Esta expresin tambin se conoce
como la ecuacin de Chapman Enskog.
3. LISTA DE MATERIALES
Materiales

Cantidad

Equipo Celda de Arnold

Teodolito

Manmetros

Termmetro

Cronometro

Tabla 1. Lista de materiales

4. PROCEDIMIENTO EXPERIMENTAL
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.

Lavar el capilar primero con agua destilada.

Secar el capilar perfectamente, mediante aire a presin o por calentamiento.
Con una jeringa adicionar sustancia (acetona), en el capilar gota a gota.
Colocar el tapn, luego colocar el tubo de difusin en la cuba centrado y en posicin vertical.
Conectar la toma de presin del tubo de difusin al manmetro de mercurio
Permitir la entrada de aire al tubo de difusin
Llenar la cuba con agua hasta un nivel ligeramente superior al tubo de difusin.
Encender control de temperatura y fijar T 25 C
Iniciar circulacin de aire: Abrir lentamente las vlvulas correspondientes
Ajustar flujo para que la cada de presin de mximo 3 mm de agua en el medidor de orificio
Se debe garantizar flujo laminar.
Estabilizados flujo y presin, dejar transcurrir de 10 a 15 min para alcanzar estado
estacionario.
13. Colocar el teodolito aproximadamente a un metro de distancia de la cuba y nivelarlo por
medio del tornillo de ajuste fino.
14. Si la burbuja del nivel de alcohol est centrada, el teodolito est nivelado, si no, colocar de
nuevo el teodolito.
15. Alinear el retculo del ocular con la parte inferior del menisco, tomando como referencia la
mitad del lente.

LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

16. Tomar la lectura superior de la trayectoria de difusin, que coincide con la parte final de la
seccin del capilar y no en la unin de este con el tubo.
17. Tomar la lectura del nivel del menisco en el capilar, poner en funcionamiento el cronmetro.
18. Registrar cada 15 min. el nivel del menisco en el capilar hasta completar mnimo 7 lecturas
19. Retirar agua de la cuba hasta alcanzar un nivel inferior al tubo de difusin.
20. Suspender el flujo de aire, cerrando vlvulas correspondientes.
21. Limpiar el capilar, a menos que no se necesite de nuevo. La limpieza se puede realizar con un
algodn y un alambre.

Colocar tapn y luego en la cuba centrado y posicin vertical

capilar
calentando
Adicionaro sustancia
a presin (acetona),al capilar gota a gota
Lavar capilar con agua y Secar
sustancia
escogida
Conectar manmetro al tub

Ajustar flujo para una cada de presin max 3 mm agua

Llenar cuba con agua hasta nivel superior al tubo de difusin
Iniciar circulacin de aire
Fijar temperatura 25C
Permitir entrada aire a tu

Colocar teodolito a 1 m de la cuba y nivelarlo

con tornillo
ajuste
fino est centrada
La burbuja
del nivel
de alcohol
Garantizar flujoEsperar
laminar10 a 15 min para estado estacionario

No

Si

Alinear retculo ocular con la parte inferior del menisco,

Tomardel
lectura
decapilar.
difusin, que coincide con la parte final de la seccin
Tomar lectura
nivelde
deltrayectoria
menisco en
Poner en funcionamiento cronometro

Tomar datos de nivel menisco en capilar cada15 min hasta 7 lecturas

Figura 2. Diagrama de flujo del proceso experimental

LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

5. EQUIPO DE COEFICIENTE DE DIFUSION

El equipo consta de una cuba donde se encuentra alojada cierta cantidad de agua, con el fin de mantener la
temperatura lo ms constante posible, es decir, no permitir muchas variaciones en sta, que afecten al gas
que est circulando, en este caso el aire. El tubo de vidrio que se encuentra perpendicular a la cuba, consta
de las siguientes conexiones: una para el termmetro, otra para la medicin de la presin que est
conectada con un manmetro de agua y la ltima localizada para el capilar. Para el flujo del gas (aire), se
encuentra un dispositivo que consta de tres pequeas columnas empacadas, de donde, el aire se dirige a
dos recuperadores de partculas, donde el flujo de aire despus de estos, se une en una sola corriente y
culmina su camino en la cuba para tener contacto con la sustancia problema (acetona) que se encuentra en
el capilar (ver figura 3).

Figura 3. Celda Arnold

6. TABLAS DE DATOS

t (min)
0
10
20
30
40
50
60
70
80
90
100

t C
Z1 (cm) h1(cm)

LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

Tabla 2: datos experimentales para la determinacin del coeficiente de difusividad

de la mezcla aire-acetona en una celda de Arnold

Temperatura del sistema (C)

Presin del sistema (mm Hg)
Cada de presin a travs del orificio(mm
H2O)
Diametro del capilar
Nivel de referencia z1 (cm)
Concentracin liq. (P/P)
Tabla 3: datos experimentales complementarios
7. ECUACIONES A UTILIZAR

Con los datos obtenidos de la tabla 2, se hace una regresin lineal de la diferencia de alturas contra el
tiempo ya que la ecuacin es:

z=

2 D AB P ( y A 1 y A 2 )
t
C AL RT
y BML

(19)

Donde la pendiente de la recta es:

2 D AB P y A 1 y A 2
=m
C AL RT
y BML

Y donde
tiempo t.

z=z2t z2

(20)

, siendo zto la altura tomada en el tiempo igual a 0, y z t la altura tomada en el

Como se puede ver en la ecuacin 19, existen muchos datos que no se tienen, lo cual hay que
determinarlos. Para hallar la concentracin de la sustancia en el lquido se usa la siguiente ecuacin:

LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

( 100xggsln ) ( gmLsln )1000 mLL

g
Peso molecular (
gmol )

%pureza
C AL =

(21)

Para determinar la composicin de los gases en la interfase se asume que el equilibrio tiene un
comportamiento ideal debido a la pureza del reactivo:

y A 1=

Pvap A
PT

(22)

La presin de vapor se halla con la ecuacin de Antoine, con las correspondientes constantes de la
acetona:

ln Pvap =A

B
T +C

(23)

La composicin de yA2 es igual a 0.

Ya teniendo estos datos se puede despejar el coeficiente de difusin de la ecuacin (19):

D AB=

y BML RT C AL z
2 ( y A 1 y A 2) t

(24)

Otras formas para hallar los coeficientes de difusin son con las ecuaciones semiempiricas

LABORATORIO DE PROPIEDADES TERMODINMICAS Y DE TRANSPORTE

La ecuacin de Ecuacin de Hirschfelder- Bird- Spotz que es la ecuacin (15); La ecuacin de Ecuacin
de Wilke Lee que es la ecuacin (16); Ecuacin de Fuller, Shettler y Giddings que es la ecuacin (17); y
el valor terico con la ecuacin (18).
El error del coeficiente de difusividad es:

%Error=

D AB teorico D AB experimental
100
D AB teorico

REFERENCIAS BIBLIOGRAFICAS
McCABE, W.; SMITH, J.; HARRIOT, P. Operaciones Unitarias en Ingeniera
Qumica. McGraw Hill. Sexta Edicin. Ciudad de Mxico. Mxico. 2002. Pags.
559 561. Apndices 18 y 19.
POLING, B.; PRAUSNITZ, J.; O'CONELL, J. The Properties of Gases and Liquids.
McGraw Hill. Fifth Edition. New York. USA. 2001.
TREYBAL, R. E., Operaciones de Transferencia de Masa. 2a edicin. McGraw
Hill. Ciudad de Mxico, Mxico. 1991.
BETANCOURT RAMIRO. Guas de Laboratorio de Operaciones Unitarias III.
Primera Edicin.
Universidad Nacional de Colombia Sede Manizales.
Manizales 2001.