UNIVERSIDAD NACIONAL AUTONOMA DE

MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
DEPARTAMENTO DE INGENIERIA QUIMICA

L-9 DETERMINACIÓN DE COEFICIENTES DE DIFUSIÓN

JAIRO PÉREZ MONTES

Carrera de Ingeniería Química, Facultad de Estudios Superiores Zaragoza, Universidad

Nacional Autónoma de México
Batalla 5 de Mayo S/N, Esq. Fuerte de Loreto, Col. Ejército de Oriente, Iztapalapa.
C.P.:09230, Ciudad de México, Teléfono:(55) 56230731 FAX:(55)57736331

RESUMEN

En la realización de este protocolo experimental se inicio con la

familiarización y verificación del funcionamiento del equipo. Posteriormente se
llenó el capilar de la celda a la mitad de la celda del fluido empleando el
embudo y el vaso de precipitados, manteniendo cerrada la llave para drenar y
abierta la llave de paso, para eliminar prácticamente todas las burbujas que se
hayan formado. De esta manera se continúo con el llenado.
Se mantuvo constante la temperatura del baño (temperatura ambiente para la
primera determinación). Una vez realizado lo anterior como lo marca el
protocolo experimental se espero el tiempo necesario para alcanzar el
equilibrio térmico (20 min aproximadamente). Después se registro la altura del
nivel del líquido en el tubo capilar en las condiciones iniciales.
Prosiguieron las observaciones de la altura del líquido en el capilar cada 30
minutos por un período de 2.5hrs.
Para concluir la experimentación se realizo el drenado del agua de recirculación
y el fluido de trabajo.

Palabras Clave: Celda de Arnold, equilibrio térmico, coeficiente de difusión,

difusión, transporte molecular.

ABSTRACT

In the realization of this experimental protocol began with the

familiarization and verification of the operation of the team. The cell capillary
was then filled to the middle of the fluid cell using the funnel and the beaker,
keeping the key closed to open and open the stopcock, to eliminate all the
bubbles that have formed. In this way the filling is continued.
The bath temperature is kept constant. Once the above was done as would be
the experimental protocol, the necessary time to reach the thermal equilibrium
was expected (approximately 20 minutes). After recording the liquid level in
the capillary tube in the initial conditions.
Observations of the height of the liquid in the capillary continued every 30
minutes for a period of 2.5 hours.
 UNIVERSIDAD NACIONAL AUTONOMA DE
MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
DEPARTAMENTO DE INGENIERIA QUIMICA

To conclude the
experiment, the progress of the recirculation water and the working fluid were
made.

Keywords: Arnold cell, thermal equilibrium, diffusion coefficient, diffusion,

molecular transport.

1. INTRODUCCIÓN

El coeficiente de difusión es un parámetro fundamental para describir el

transporte de una especie en otra mediante el fenómeno de difusión en una
misma fase, impulsada por la diferencia de concentración de una especie
particular y el movimiento aleatorio de sus moléculas en dirección opuesta a
ese gradiente de concentración, para lo cual se han desarrollado modelos de
transferencia de masa y correlaciones para una gran variedad de sistemas que
pueden estar presentes a nivel de laboratorio o industrial, los cuales pueden
ser binarios o multicomponentes, así como pueden también incluir las
interacciones de la fase en donde se desarrolla la difusión con el medio que la
rodea.

2. ANTECEDENTES

Los estudios en el área de transferencia de masa datan de los mismos tiempos

de origen de la ingeniería como ciencia; y en la medida en que crecieron las
necesidades para el desarrollo de tecnologías en el terreno de los procesos
químicos, se profundiza su estudio a fin de caracterizarla. Cuantificar el
fenómeno de transferencia de masa depende de muchos elementos y
parámetros interrelacionados entre sí, por ejemplo: fases donde ocurre,
número de especies que intervienen, viscosidad de la mezcla, concentración,
etc.
Aunque la transferencia de masa por difusión y convección no son eventos
aislados, los trabajos realizados al respecto se basan en modelos donde la
intervención de uno u otro proceso se obvia por conveniencia; en lo que
respecta a difusividad las ecuaciones desarrolladas por Fick (primera y
segunda ley de Fick) dan respuesta a una innumerable cantidad de procesos
difusionales, pero por la complejidad de cómo otros procesos se llevan a cabo
se hace necesario abordarlos de diferente manera; por ejemplo, sustentándose
en la teoría cinética de los gases, importantes trabajos como los de
Hirschfelder-Bird-Spotz (según Bird, 1998) modelan procesos difusionales con
poco margen de error; sin embargo, la relación propuesta por estos
investigadores y muchas otras ecuaciones trabajadas con el mismo marco
teórico, para la determinación de difusividades en gases, no se pueden
emplear si no se conocen los diámetros de las moléculas.
Existen diversos trabajos basados en las leyes de Fick, adaptados a situaciones
particulares como por ejemplo: evaporación de un líquido puro en un gas
estancado o mezclas de líquidos miscibles. Con relación a la determinación de
 UNIVERSIDAD NACIONAL AUTONOMA DE
MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
DEPARTAMENTO DE INGENIERIA QUIMICA

difusividades en gases dos

trabajos muy didácticos son presentados por Kwon, et al, (2002) y Nirdosh et

al, (2000). En sus documentos proponen un desarrollo matemático apoyados

en la primera ley de Fick junto con un balance de masa, cuando se evaporan
líquidos orgánicos puros (como pentano, hexano, n-heptano, acetona y otros)
en tubos abiertos; sus resultados experimentales comparados con los teóricos
presentan un error máximo de
10%.
Amplios estudios para la determinación de coeficientes de transferencia de
masa se han hecho en columnas tipo pared mojada. Gilliland, Sherwood y
Woerstz, Chilton y Colburn entre otros, presentan trabajos rigurosos donde se
han obtenido resultados interesantes que han podido ser correlacionados. El
sistema agua-aire ha sido el más estudiado (Perry, 1972 ; Badger et al, 1964).

3. OBJETIVO

Determinar experimentalmente el coeficiente de difusión de líquidos en aire.

Determinar la variación del coeficiente de difusión con respecto a la
temperatura.
Comparar los resultados experimentales con los obtenidos a partir de
correlaciones teóricas.

4. HIPÓTESIS

Si se aborda correctamente el protocolo experimental, determinando

experimentalmente tanto el coeficiente de difusión de líquidos en aire como la
variación del coeficiente de difusión con respecto a la temperatura, entonces al
realizar una comparación de los resultados experimentales con los datos
obtenidos a partir de correlaciones teóricas se observara que la variación entre
unos datos y otros será mínima.

5. MARCO CONCEPTUAL DE REFERENCIA

Difusión.
La difusión es un fenómeno de transferencia de masa que causa la distribución
de especies químicas generando uniformidad en un medio. Este fenómeno se
describe bajo el estudio de la transferencia de masa, que para una especie
química describe la evolución de su concentración en el tiempo y espacio. Al
haber una diferencia de concentración de la especie en el medio, se genera un
gradiente de concentración que será la fuerza impulsora del proceso de
difusión hasta conseguir una concentración uniforme, en donde las moléculas
disueltas estarán igualmente distribuidas en la solución.
 UNIVERSIDAD NACIONAL AUTONOMA DE
MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
El DEPARTAMENTO DE INGENIERIA QUIMICA

modelo de
transferencia de masa por difusión describe el proceso mencionado mediante
la Ley de Fick:
dC A
J a =−D AB
dZ

En donde J a corresponde al flux difusional definido como las moles de una

especie A (que se difunden en dirección contraria al gradiente de
concentración) por unidad de área transversal por tiempo.
dCA
El término corresponde al gradiente de concentración en el sentido Z, la
dZ
cual tiene un valor negativo debido al cambio de mayor a menor concentración
para la especie A.

Coeficiente de difusión.
El término D AB de la ley de Fick corresponde al coeficiente de difusión o
difusividad de la especie A sobre la especie B (unidades de área sobre tiempo),
la cual es una medida de que tan rápido puede A difundirse en B. Es pertinente
aclarar que para un par de especies gaseosas A y B: D AB=D BA , y es una
cantidad fuertemente dependiente de la temperatura, mientras que estos
coeficientes de difusión no serán iguales si ambas especies son liquidas, así
como tampoco la cantidad D AB en fase liquida no será igual que D AB en fase
gaseosa.
Existen varias maneras de calcular los coeficientes de difusión dependiendo de
las características del sistema y de las especies involucradas. Para los sistemas
gaseosos binarios pueden estimarse mediante la ecuación de Hirschfelder
usando los parámetros de Lennard-Jones (aplicable a sustancias no polares a
presiones menores de 25 atm):
3
1 1 12
0.00185∗T 2∗( + )
PM A PM B
D AB= 2
P∗σ AB ∗❑D

Siendo T la temperatura absoluta (K), P la presión absoluta (atm), PM el peso

molecular de la respectiva sustancia (g/mol) y los parámetros con expresiones
reportadas en la literatura: σ AB el diámetro de colisión entre A y B (Amstrongs),
y ❑D la integral de colisión.

Celda de difusión de Arnold.

La celda de Arnold es un montaje experimental para determinar un coeficiente
de difusión binario para un sistema gaseoso. Consta de un tubo angosto que se
llena parcialmente con un líquido puro A a temperatura constante usando un
baño de agua y a presión constante. En la parte superior del montaje se hace
fluir por un tubo un gas B, el cual debe tener una solubilidad sobre A ( D AB )
despreciable, asi como debe ser inerte químicamente con A, por lo que se toma
 UNIVERSIDAD NACIONAL AUTONOMA DE
MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
al DEPARTAMENTO DE INGENIERIA QUIMICA
gas
de la especie B como
estancado asumiendo que su solubilidad en A es prácticamente cero.
El componente A se vaporiza y se difunde en la fase gaseosa de B. Se ha
definido un volumen de control en el que ocurre la difusión de A limitado entre
las alturas Z1 (la superficie del líquido de la especie A) y Z2 (la frontera entre
tubo angosto vertical y el tubo horizontal por donde fluye la especie B (Figura
1).
Es de destacar que en la altura Z1 la fracción molar de A ( Y A ) es máxima e
igual a la presión de vapor de A ( PV ) dividido la presión total del sistema (P). En
la altura Z2:Y A =0 si se asume que en tal punto B fluye a tal velocidad que ha
arrastrado cualquier cantidad de A a esta altura.

Fig 1. Esquema general de la celda de Arnold.

Puede determinarse que para este tipo de sistema el flux de A es constante a

partir del modelo de transporte molecular (consideración de Ley de Fick junto
con la ecuación de continuidad).
La ecuación de continuidad está definida como:

dC A
∇ N A+ −R A =0
dT

En donde los tres términos del lado izquierdo son respectivamente: el cambio
del flux en los bordes del volumen de control, la rapidez de acumulación de A
en el volumen de control y la rapidez de producción química de A dentro del
volumen de control. Al no haber acumulación en el volumen de control ni
reacción química, la ecuación de continuidad se reduce a:

∇ N A=0
 UNIVERSIDAD NACIONAL AUTONOMA DE
MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
Que DEPARTAMENTO DE INGENIERIA QUIMICA
en
coordenadas
cartesianas:

dC A
∇ N A= =0
dT

Por lo cual N A en la dirección Z es constante para este sistema.

6. MATERIAL Y EQUIPO

MATERIAL
(Para un sistema)
 Mangueras de hule látex
 Vasos de precipitado de 100 ml
 1 Cronómetro

 1 Termómetro escala 0 -100 °C con graduación 1 °C

 1 Embudo de tallo corto
 1 Pinza de Mohr

REACTIVOS
 50 ml Acetona
 50 ml Etanol

EQUIPO
 1 Baño de temperatura constante
 1 Aparato para determinación de coeficientes de difusión (equivalente a
celda Arnold)
 1 Recirculador

SERVICIOS
 Electricidad
 Agua
 Aire

7. METODOLOGÍA EXPERIMENTAL

Familiarizarse con el aparato.

 UNIVERSIDAD NACIONAL AUTONOMA DE
MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
DEPARTAMENTO DE INGENIERIA QUIMICA

Fig. 2 Equipo para determinación de coeficientes de difusión.

Mantener cerrada la llave Cerrar la llave de paso cuando el

Llenar el capilar de la capilar se encuentre
v celda. para drenar y abierta la
aproximadamente a 1 cm de la
llave de paso. salida del agua.

Temperatura ambiente Mantener constante la

Eliminar todas las burbujas
para la primera. temperatura del baño
que se hayan formado.

Iniciar el flujo de agua en

Esperar para alcanzar el equilibrio las entradas
térmico (20min.) correspondientes

Anotar la altura del nivel

del líquido en el tubo
capilar en las condiciones
iniciales.
 UNIVERSIDAD NACIONAL AUTONOMA DE
MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
DEPARTAMENTO DE INGENIERIA QUIMICA

8. RESULTADOS

T (°C) 24
ambiente
ρ etanol 0.789
(g/cm3)
PM etanol 46.07
(g/mol)
PM aire 29
(g/mol)
Cm 45.0456
C 0.94
YA 0.94
YB 0.06
Tabla 1 Condiciones iniciales

Temperatur tiempo Altura Z

a (K) (s) (cm)
296.15 0 0
296.15 30 0.01
297.15 60 0.02
296.15 90 0.03
295.65 120 0.035
 UNIVERSIDAD NACIONAL AUTONOMA DE
MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
DEPARTAMENTO DE INGENIERIA QUIMICA

295.15 150 0.045

Tabla 2 Datos experimentales

1.00E-03
DAB 9.00E-04
Temperat DAB teórico DAB 8.00E-04
ura (K) (cm2/s) experimental
7.00E-04
(cm2/s)
296.15 1.20E-05 0 6.00E-04
DAB TEORICO
296.15 1.20E-05 3.33E-06 5.00E-04

DAB
DAB EXPERIMENTAL
297.15 1.20E-05 6.66E-06 4.00E-04
296.15 1.20E-05 0.00001 3.00E-04
295.65 1.20E-05 1.02E-05
2.00E-04
295.15 1.20E-05 9.35E-04
Tabla 3 DAB experimental y teorico 1.00E-04
0.00E+00 Temperetaura (K)
295 295.5 296 296.5 297 297.5
Grafico 1. DAB exp & teorico VS T°

CORRELACION GILLILAND Correlación GILLILAND

Temperatura (K) DAB
1.40E-01
experimental
296.15 1.38E-01 1.39E-01
296.15 1.38E-01 1.38E-01 DAB EXPERIMENTAL
DAB

297.15 1.39E-01 1.37E-01

296.15 1.38E-01 1.36E-01
295.65 1.37E-01 295 296 297 298
295.15 1.37E-01 Temperatura (K)
Tabla 4 Correlación Gilliliand

CORRELACIÓN FULLER, SCHETTLER Y GIDDINGS correlación FULLER,

Temperatura (K) SCHETTLER Y GIDDINGS
296.15 1.58E-01 1.60E-01
296.15 1.58E-01
1.59E-01
297.15 1.59E-01 DAB experimental
1.58E-01
DAB

296.15 1.58E-01
295.65 1.57E-01 1.57E-01
295.15 1.57E-01 1.56E-01
Tabla 5 Correlación Fuller, Schettler y Giddings 295 296 297 298
Temperatura (K)
 UNIVERSIDAD NACIONAL AUTONOMA DE
MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
DEPARTAMENTO DE INGENIERIA QUIMICA

CORRELACIÓN
DAB
CINÉTICA DE EDO CORRESPONDIENTES 5.80E-06
Temperatura (K) DAB experimental 5.78E-06
5.76E-06
296.15 5.75E-06 5.74E-06
DAB

DAB
296.15 5.75E-06 5.72E-06
297.15 5.78E-06 5.70E-06
5.68E-06
296.15 5.75E-06 5.66E-06
295.65 5.71E-06 295 296 297 298
295.15 5.71E-06 Temperatura (K)
Tabla 6 Correlación de cinética de edo.

9. ANÁLISIS DE RESULTADOS

En la tabla 3 grafico 1 se puede apreciar el comportamiento de los coeficientes

de difusión a diferentes temperaturas.
También se observa que los datos teóricos son diferentes de los
experimentales.
Mientras que con las correlaciones posteriores se aprecia gran similitud en el
comportamiento gráfico, aunque la tendencia no es lineal se aprecia que hay
cierta proporcionalidad entre la temperatura y el coeficiente de difusión ( D AB ¿

10. CONCLUSIÓN
En los gases la difusividad aumenta con la temperatura y disminuye al
aumentar la presión, también disminuye al aumentar el peso molecular, el
comportamiento anterior es similar en líquidos, en contraste la presión no
afecta considerablemente.
La velocidad promedio de la molécula depende de la temperatura. Como la
molecula viaja en una trayectoria en la cual se mueve durante cierto tiempo
“rapidez de la difusión”, solo es una pequeña fracción de la longitud de su
trayectoria real.Por esta razón la rapidez de difusión es muy pequeña aunque,
podría aumentar con un descenso de presión que reduciría el numero de
choques y un incremento en la temperatura, que aumentaría la velocidad
molecular.
 UNIVERSIDAD NACIONAL AUTONOMA DE
MEXICO
FACULTAD DE ESTUDIOS SUPERIORES
ZARAGOZA
11. DEPARTAMENTO DE INGENIERIA QUIMICA

BIBLIOGRAFÍA

1. Forero, M. (2003). Guía práctica difusividad. Universidad Nacional de

Colombia. Pp. 1-10
2. Betancourt, R. (2001). GUÍAS DE LABORATORIO DE OPERACIONES
UNITARIAS III. 1ra edición. Universidad nacional de Colombia sede
Manizales. Pp. 7-26
3. Welty, J. (2005). Fundamentos de Transferencia de Momento, Calor y
Masa.
4. Multiphysics Cyclopedia, «Diffusion Coefficient,» 2015. [En línea].
Available: https://www.comsol.com/multiphysics/ diffusion-coefficient.
5. McCabe, W (2005). Operaciones Unitarias en Ingeniería Química.