UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

(Universidad del Perú, Decana de América)

FACULTAD DE QUÍMICA, INGENIERÍA QUÍMICA

Escuela Académico-Profesional de
Ingeniería Química

DEPARTAMENTO ACADÉMICO DE OPERACIONES UNITARIAS

Laboratorio de Ingeniería Química II

“DIFUSION ACETONA-AIRE A 40°C”

PROFESOR : Dr. Condorhuaman Ccorimanya, Cesario

ALUMNOS : GONZALES CARAMANTIN, Carlos

MAYTA TICONA, Humberto
MEJIA MAITA, Frank
MONTENEGRO ATUNCA, Diego

Ciudad Universitaria, junio del 2017

 DIFUSION

TABLA DE CONTENIDO

Página

I.RESUMEN 3
II.INTRODUCCIÓN 4
III.PRINCIPIOS TEÓRICOS 5
IV.DETALLES EXPERIMENTALES 13
V.TABLAS DE DATOS Y RESULTADOS 15
VI.DISCUSIÓN DE RESULTADOS 17
VII.CONCLUSIONES 18
VIII.RECOMENDACIONES 18
IX.BIBLIOGRAFÍA 19
X.APÉNDICE I: EJEMPLO DE CÁLCULOS 20
XI.APÉNDICE II: GRÁFICAS 25

GRUPO E 2
 DIFUSION

RESUMEN

El objetivo del presente informe fue determinar el coeficiente de difusividad de la

acetona a través del aire a una temperatura constante de 40°C. Así mismo, comparar la
difusividad obtenida, con la difusividad de la acetona a 40°C proporcionada por el manual
y la teorica encontrada.

Para esta práctica se utilizó el equipo didáctico Armfield, el método empleado

para obtener el coeficiente de difusión fue el de Winklemann obteniendo el valor del
coeficiente de difusividad para la acetona a una temperatura de 40°C es 0.03172 cm2/s. y
comparado con la difusividad de la acetona a 40°C, que tiene un valor de 0.11068 cm2/s
según el manual, además el teórico hallado es de 0.11591 𝑐𝑚2 /s, según el
Yaws’Handbook.

GRUPO E 3
 DIFUSION

INTRODUCCIÓN
La difusión es el movimiento, debido a un estímulo físico, de un componente a
través de una mezcla. La principal causa de la difusión es la existencia de un gradiente de
concentración del componente que difunde. Este movimiento es aprovechado en las
operaciones de transferencia de materia. Por ejemplo, en absorción de gases el soluto se
difunde a través de la fase gaseosa hacia la interfase y a través de la fase líquida desde la
interfase. En destilación el componente menos volátil difunde a través de la fase líquida
hacia la interfase y desde ésta hacia el vapor. El componente menos volátil difunde en
sentido contrario y pasa a través del vapor hasta el seno del líquido. En lixiviación la
difusión del soluto a través de la fase sólida va seguida de la difusión en el seno del
líquido. En extracción líquido-líquido el soluto difunde a través de la fase de refinado
hacia la interfase y después hacia el interior de la fase extracto. También se aprovecha en
otras operaciones de transferencia de materia tales como cristalización, humidificación,
secado, etc.
La difusividad o coeficiente de difusividad, DAB, es una medida de la movilidad
de difusión. En la literatura es posible encontrar valores del coeficiente de difusividad
para diferentes pares de sustancias así como ecuaciones basadas en la teoría cinética y
empíricas que nos permiten calcular dichos valores. Sin embargo, Winklemann
proporciona un método experimental para determinar la difusividad de líquidos volátiles
a través de gases.
La importancia de la difusión radica en que es una operación de bajo costo de
operación, de inversión y de instalación lo cual permite una mayor eficiencia mecánica y
por tanto el consumo de potencia es menor. Así por ejemplo entre un difusor y un molino
de la misma capacidad en relación a la caña a procesar, muestra que un difusor cuesta
más o menos el mismo precio. Pero si esta comparación es hecha en relación a las
toneladas de sacarosa extraídas por hora, el difusor tiene ventajas incalculables.
La aplicación de la difusión es también importante en numerosas industrias como:
la industria azucarera, aplicación de tintes penetrantes, en la industria alimentaria
(envases termoplásticos), en el sector de la construcción (penetración de oxígeno en
membranas termoplásticas para edificaciones).

GRUPO E 4
 DIFUSION

PRINCIPIOS TEÓRICOS
1. FUNDAMENTOS DE LA DIFUSIÓN
La difusión es el movimiento, bajo la influencia de un estímulo físico, de un
componente individual a través de una mezcla. La causa más frecuente de la difusión es
un gradiente de concentración del componente que difunde. Un gradiente de
concentración tiende a mover el componente en una dirección tal que iguale las
concentraciones y anule el gradiente. Cuando el gradiente se mantiene mediante el
suministro continuo del componente de baja y alta concentración, existe un flujo en
estado estacionario del componente que se difunde. Esto es característico de muchas
operaciones de transferencia de masa. La transferencia de masa puede tener lugar tanto
en fase gaseosa como liquida o en ambas simultáneamente; por ejemplo, cuando se
remueve amoniaco de un gas por medio de absorción en agua en una columna empacada,
en cada punto en la columna un gradiente de concentración en la fase gaseosa provoca
la difusión del amoniaco a la interface gas-liquido, donde se disuelve, y un gradiente en
la fase liquida provoca la difusión en la masa liquida. Cuando se extrae un soluto de un
líquido, los gradientes se invierten; aquí la difusión conduce al soluto desde una masa
liquida a la interface y de ahí al interior de la fase gaseosa.
Aunque la causa habitual de la difusión es un gradiente de concentración, en
ocasiones la difusión también puede ser originada por un gradiente de actividad, como
en la osmosis inversa, por un gradiente de presión, por un gradiente de temperatura o
por la aplicación de una fuerza externa como en el caso de una centrífuga.
La difusión no está restringida a la transferencia molecular a través de capas
estacionarias de sólido o fluido. También tiene lugar cuando se mezclan fluidos de
diferentes composiciones.

2. DIFUSION MOLECULAR
Si una solución es completamente uniforme con respecto a la concentración de
sus componentes, no ocurre ninguna alteración; en cambio, si no es uniforme, la solución
alcanzara la uniformidad espontáneamente por difusión, ya que las sustancias se moverán
de un punto de concentración elevada a otro de baja concentración. La rapidez con la cual
el soluto se mueve en cualquier punto y en cualquier dirección dependerá, por tanto, del
gradiente de concentración en ese punto y en esa dirección. Para describir

GRUPO E 5
 DIFUSION

cuantitativamente este proceso, se necesita una medida apropiada de la rapidez de

transferencia.
La rapidez de transferencia puede describirse adecuadamente en función del flujo
molar, o moles/(tiempo)(área), ya que el área se mide en dirección normal a la difusión.
Sin embargo, aunque una solución no uniforme contenga solo dos componentes, estos
deberán difundirse si se quiere alcanzar la uniformidad. Surge entonces la necesidad de
utilizar dos fluxes para describir el movimiento de un componente: N, el flux relacionado
con un lugar fijo en el espacio, y J, el flux de un compuesto con relación a la velocidad
molar promedio. El primero es importante al aplicarse el diseño del equipo; el segundo
es característico a la naturaleza del componente. Por ejemplo, un pescador estaría más
interesado en la rapidez con la cual nada un pez en contra de la corriente para alcanzar el
anzuelo (análogo a N); la velocidad del pez con relación a del arroyo (análogo a J) es
característica natatoria del pez.
Así, la difusividad, o coeficiente de difusión, DAB de un componente A en solución
en B, que es una medida de la movilidad de difusión, se define como la relación de su
flux JA y su gradiente de concentración.
x A
J A  cDAB (1)
z
Esta es la primera ley de Fick, en este caso para la dirección z. Esta ecuación
establece que la especie A difunde (se mueve con relación a la mezcla) en la dirección
decreciente a la fracción molar de A (de la misma forma que el calor fluye por conducción
en dirección de temperatura decreciente). El flux molar NA relativa a coordenadas
estacionarias:
x A
N A  x A ( N A  N B )  cDAB (2)
z
Flux molar Flux molar que Flux de difusión
relatio a resulta del relativo a la
una flujo global veloc. molar
superficie promedio
fija

Esta ecuación indica que el flux de difusión NA con relación a un eje de

coordenadas estacionarias, es la resultante de dos magnitudes vectoriales: el vector xA(NA
+ NB), que es el flux molar de A que resulta del movimiento global del fluido, y el vector
JA que es flux de difusión relativo al a velocidad molar promedio.

GRUPO E 6
 DIFUSION

Por lo tanto, los términos de flujo global y de difusión de la ecuación (2) son del
mismo sentido para la especie A (puesto que se difunde a favor de la corriente) y de
sentido contrario para la especie B (debido a que B se difunde en contracorriente).
Las unidades de la difusividad DAB son cm2seg-1 o m2hr-1. Obsérvese que la viscosidad
cinemática ν y la difusividad térmica α también tiene las mismas unidades. La analogía
de estas tres magnitudes se deduce de las siguientes ecuaciones para densidades de flujo
de masa, cantidad de movimiento y energía, en sistemas unidimensionales:
d
 J Ay   D AB (  A ) , (Ley de fick para  constante) (3)
dy

d
  yx   ( v x ) , (Ley de newton para  constante) (4)
dy

d
 q y   ( C p T ) , (Ley de Fourier para C p constante) (5)
dy

Estas ecuaciones establecen respectivamente, que:

 El transporte de materia tiene lugar a causa de una gradiente de concentración.
 El transporte de cantidad de movimiento tiene lugar a causa de una gradiente de
concentración de cantidad de movimiento.
 El transporte de energía se lleva a cabo por una gradiente de concentración de
energía.

Estas analogías no pueden aplicarse a problemas bi y tridimensionales, puesto que  es

una magnitud tensorial con nueve componentes, mientras que JA y q son vectores con tres
componentes.

3. METODO DE WINKELMANN
Consideremos el sistema de difusión que se presenta en la fig. N° 1. El líquido A
se está evaporando en el seno del gas B, e imaginemos que mediante un artificio es posible
mantener el nivel del líquido en z=z1. La concentración de la fase gaseosa, expresada en
fracción molar, exactamente en la interface liquido-gas es xA1. Eso quiere decir que xA1 es
la relación entre la presión de vapor de A y la presión total, suponiendo que A y B forman

GRUPO E 7
 DIFUSION

una mezcla gaseosa ideal. Finalmente se supone que la solubilidad de B en el líquido A

es despreciable.
Por la parte superior del tubo (para z=z2) circula lentamente una corriente de
mezcla gaseosa A-B cuya concentración es xA2, de forma que la fracción molar de A en la
parte superior de la columna permanece constante e igual a xA2. Se supone que todo el
sistema se mantiene a presión y temperatura constantes y que los gases A y B se
comportan como ideales.
Cuando el sistema alcanza el estado estacionario, existe un movimiento neto de A
alejándose de la superficie de evaporación, mientras que para el vapor de B se tiene que
NBz=0. Por lo tanto, puede utilizarse NAz correspondiente a la ecuación (2).
Despejando NAz, se obtiene:
cDAB dx A
N Az   (6)
1  x A dz
Aplicando un balance de materia en estado estacionario a un incremento de altura ∆z, se
llega a:
SN Az z  SN Az z  z
0 (7)

En la que S es el área de la sección transversal de la columna. Las distribuciones

de concentración se representan en la fig. N°1. El examen de esta curva nos indica que la
pendiente dxA/dz no es constante, a pesar de serlo la densidad de flujo molar NAz.
Si bien los perfiles de concentración son útiles para describir los procesos de
difusión, lo que generalmente interesa en los cálculos ingenieriles es la concentración
media o la densidad de flujo de materia en una superficie. Por ejemplo, la concentración
media de B comprendida entre z=z1 y z=z2 es:
x B 2  x B1
x B ,med  (8)
ln( x B 2 / x B1 )
Es decir que el valor medio de xB es la media logarítmica de los valores extremos.
La velocidad de transferencia de materia en la interface liquido-gas, es decir la velocidad
de evaporación, se obtiene utilizando la ecuación (6):
cD AB dx A cD AB dx B cD AB x 
N Az    ln  B 2  (9)
z  z1
1  x A1 dz z  z1 x B1 dz z  z1 ( z 2  z1 )  x B1 

GRUPO E 8
 DIFUSION

Las ecuaciones (8) y (9) pueden combinarse para obtener otra expresión de velocidad de
transferencia de materia.
cD AB
N Az  ( x A1  x A2 ) (10)
z  z1
( z 2  z1 ) x BM
Esta ecuación indica la forma en que está relacionada la velocidad de transferencia de
materia con una fuerza motriz característica de concentración xA1-xA2. Por lo tanto
podemos reordenar la ecuación (10), haciendo las siguientes restricciones:
 Usando un sistema unidimensional (en la dirección z).
 La fuerza motriz característica de concentración xA1-xA2 sea CA la cual es la
concentración en la interface.
 Las fracciones molares para la especie B sean reemplazados en términos de
concentración CB1, CB2 y CBM.

 C  C 
N A  D A  T  (11)
 L  C BM 
Donde:
D: Difusividad (m2/s).
CA: Concentración de saturación en la interface (Kmol/m3)
L: Distancia efectiva de transferencia de masa (mm)
CBM: Media logarítmica de concentración de vapor (Kmol/m3)
CT: Concentración molar total = CA  CBM (Kmol/m3)
Sin embargo el sistema mostrado también puede trabajarse en un estado de cuasi-estado
de equilibrio de difusión; debido a que la especie A, el nivel de líquido disminuye muy
lentamente debido a su evaporación.
Primero, en lugar de mantener líquido-gas en una interface de altura constante, que
permita el nivel de líquido como para disminuir la evaporación producto, tal como se
muestra en la Fig N°1 (ya que el líquido se evapora muy lentamente), podemos utilizar el
método de cuasi-estado de equilibrio.
Lo primero es equiparar el flujo molar de evaporación de "A" de la fase líquida con el
ritmo molar de "A" entrar en la fase gaseosa.
 A dz1 cD AB
S  ( x A1  xxA2 ) S (12)
M A dt ( z 2  z1 (t )) x BM

GRUPO E 9
 DIFUSION

Aquí  A es la densidad del líquido puro A, MA es el peso molecular y z2-z1(t) es la distancia

de la interface que ha descendido en un tiempo t. En el lado derecho de ecuación (12) se
utiliza el estado de equilibrio de flujo de evaporación para evaluar la actual altura de
columna de líquido (esto es casi constante estado de aproximación).
Uno puede utilizar esta prueba para obtener la difusividad a partir de las mediciones del
nivel de líquido en función del tiempo. Como en el caso de la ecuación (12), se puede
arreglar esta ecuación en términos apropiados para su interpretación. Entonces:

  L  dL  C  C 
   D A  T  (13)
 M  dT  L  C BM 
De acuerdo a la figura N°3(a):
L
 MD  C ACT t
Lo LdL    L  C BM  tdt
0
(14)

 2MD  C ACT 
L2  L20    t (15)

 L  CBM 
Nota: L y Lo no pueden medirse exactamente, pero L - L0 se puede medir acertadamente
usando un catetómetro.
 L2  Lo 2  ( L  Lo)( L  Lo)  ( L  Lo  2Lo)( L  Lo)

  
L  L0 L  L0  2L0    2MD  CACT t (16)
  L  CBM 
Linealizando la expresión tenemos:

   C    C 
L  L0    L BM
t
  L  BM  L0 (17)
L  L0   2MD  C ACT   MDC ACT 
Donde: M: Peso molecular (kg/Kmol)
t: tiempo, s
Sea “s” la pendiente de la ecuación (22):

s
 LCBM 
(18)
2MDC ACT 
Despejando “D” tenemos:
 L C BM 
D (19)
2sMC A CT 

1
GRUPO E
0
 DIFUSION

Donde:

  1  Tabs
CT   
 Kmol.Vol  Ta

 *Kmol x Volumen = 22.414 m3/Kmol

 C B1  CT

P P 
 CB 2   a V CT
 Pa 

Fig. 1: (a) Difusión de A en estado estacionario a través de B inmóvil. (b) Forma en

que se distorsiona el perfil de concentración debido a la difusión de A.
(a)

Aire
Referencia
z=z2

NA|L+∆L
L = Lo=z1; t = 0

∆L=∆z

L = L; t = t

NA|L

Acetona

1
GRUPO E
1
 DIFUSION

(b)

Fig. N° 2: Evaporación con cuasi-estado estacionario de difusión. El nivel de

líquido disminuye muy lentamente cuando se evapora.

1
GRUPO E
2
 DIFUSION

DETALLES EXPERIMENTALES
1. DETALLES DEL EQUIPO
El equipo de trabajo es el DIFFUSION OF A GAS APPARATUS ARMFIELD, en
el cual determinaremos el coeficiente de difusión de un gas por evaporación desde su
superficie liquida. El equipo cuenta con las siguientes partes:
 Termómetro de mercurio (rango de 10 a 200°C).
 Tubo capilar en forma de “T”, el cual estaba lleno con acetona.
 Recipiente de vidrio con agua desmineralizada, que servirá como medio (baño
maría).
 Bomba de aire (compartimiento azul).
 Microscopio.
 Vernier (con una lectura de hasta dos decimales).
 Controlador de temperatura digital.

Fig. N°3: Esquema del equipo

1
GRUPO E
3
 DIFUSION

2. PROCEDIMIENTO EXPERIMENTAL

1. Se vierte la acetona en el capilar, previamente lavado (con la misma acetona). Se

enciende la bomba de aire (caja azul) el cual tenía conectado una manguerilla que
se conecta al tubo capilar.

2. Fijar la temperatura de trabajo (Para este caso 40°C) con el control de

temperatura.

3. Ajustar la altura del microscopio de modo que se visualice el capilar. Luego se

empieza a bajar poco a poco hasta encontrar el menisco de la acetona.

4. Ajustar la posición de la lente para obtener una visión clara y definida del menisco
(la imagen esta invertida). Cabe mencionar que la mesa donde está el equipo de
difusión no debe sufrir ningún tipo de movimiento o golpe pues esto puede
alterarlas condiciones de trabajo establecidas.

5. Fijar una altura en el menisco del líquido de trabajo (Lo), la escala del vernier
debe ubicarse en una posición adecuado la cual servirá de referencia, en ese
instante es el tiempo cero, y a partir del cual se realizan las mediciones en el
vernier cada 10 a 15 minutos, los cuales aumentaran en de tiempo conforme
avance el experimento.

1
GRUPO E
4
 DIFUSION

TABLAS DE DATOS Y RESULTADOS

Tabla N° 1: Propiedades la Acetona
Temperatura
40
(°C)
Pv (kPa) 56.21
3
ρ (kg/m ) 767.88
M (kg/kmol) 58.08

Tabla N° 2: Datos Experimentales de Mediciones de la Altura del Capilar

T°(C) t acumulado (s) t (ks) L(mm) Lo (mm) L - Lo (mm) t/(L-Lo) (ks/mm)

40 0 0 34 0 34 0

40 1215 0 30.1 30.1 0 0

40 2310 1.095 30.5 30.1 0.4 2.7375

40 3465 2.25 31.3 30.1 1.2 1.875

40 4677 3.462 31.9 30.1 1.8 1.923333333

40 5900 4.7 32.5 30.1 2.4 1.958333333

40 7410 6.195 33.2 30.1 3.1 1.998387097

40 8314 7.103 33.6 30.1 3.5 2.029428571

Tabla N° 3: Datos Experimentales Empleados para el cálculo de la difusión

T°(C) t (ks) L(mm) Lo (mm) L - Lo (mm) t/(L-Lo) (ks/mm)

40 2.25 31.3 30.1 1.2 1.875

40 3.462 31.9 30.1 1.8 1.923333333

40 4.7 32.5 30.1 2.4 1.958333333

40 6.195 33.2 30.1 3.1 1.998387097

40 7.103 33.6 30.1 3.5 2.029428571

1
GRUPO E
5
 DIFUSION

Tabla N°4: Dato del manual del equipo Armfield Limited (Acetona 40°C)
t (ks) L-Lo (mm) t/(L-Lo) (ks/mm)
0.00 0 0.000
3.60 2.20 1.636
7.20 4.20 1.714
11.16 6.30 1.771
15.90 8.80 1.807
19.98 10.80 1.850
23.40 12.40 1.887
78.78 34.50 2.233
83.52 36.10 2.313
87.24 37.30 2.339
91.80 38.90 2.360
97.32 40.80 2.385
101.10 42.00 2.407

Tabla N°5: Concentraciones halladas para calcular el coeficiente de difusión

C T (kmol/m3) 0.0389

C A (kmol/m3) 0.0215

C B1 (kmol/m3) 0.0389

C B2 (kmol/m3) 0.0174

C BM (kmol/m3) 0.0267

Tabla N°6: Resultados pendiente y coeficiente de difusión

m 0.0649

D AB (m2/s) 3.172E-06

D AB (cm2/s) 0.03172

D AB manual 0.11068

D AB teórico 0.11591

1
GRUPO E
6
 DIFUSION

DISCUSIÓN DE RESULTADOS

 En la tendencia lineal para determinar la pendiente “m” de la ecuación WINKLEMANN, se tuvo

que omitir alturas tomadas en el proceso experimental, debido a que no todos los puntos tienen
dicha tendencia, esto se da principalmente en los primeros puntos de la recta, esto se puede
atribuir a que al inicio no se llega al estado de equilibrio, esto se observa en las gráficas Nº 1, 2,3
y 4.

 Se trabaja en un sistema estacionario (a presión y temperatura constantes.); sin embargo el calor

de absorción no ha sido tomado en cuenta pero debe mencionarse que este calor tiende a
incrementar la temperatura interfacial. Dicho incremento no puede ser medido directamente,
pero si puede estimarse considerándose la transferencia simultanea de calor y masa dentro de
la película liquida interfacial.

 Los cálculos matemáticos son desarrollados por el método de WINKLEMANN, para lo cual se ha
asumido en primer término, como un proceso estacionario. Al mismo tiempo la difusión
ocasiona el descenso de altura de la acetona es por eso que también se asume condiciones de
un cuasi estado de equilibrio.

 Se utiliza tubos de pequeños diámetros (capilares), ya que un diámetro de mayor tamaño

implicaría una difusión en 2 direcciones, la radial y longitudinal.

1
GRUPO E
7
 DIFUSION

CONCLUSIONES
 El valor del coeficiente de difusividad para la acetona en aire a una temperatura de 40°C
es 0.03172 cm2/s, para la experiencia realizada.

 El valor del coeficiente de difusividad para la acetona en aire hallado en la experiencia,

presenta una desviación de 71,34 % con respecto al teórico encontrado, en el
Yaws’Handbook y con respecto al manual la desviación es de 72.63%.

RECOMENDACIONES

 Es recomendable que la experiencia realizada para hallar el coeficiente de difusión de la

acetona en el aire se realice a distintas temperaturas para observar su comportamiento
frente a la variación de esta asegurándose que las temperaturas trabajadas sean menores
al punto de ebullición.

 Realizar la mayor cantidad de lecturas, pero a largos intervalos de tiempo, ya que la

acetona a 40ªC demora en difundirse, esto nos permite observar mejor la linealidad de la
curva “L-L0 vs t/(L-L0)”.

 Realizar la experiencia a varias temperaturas para observar el comportamiento frente a

la variación de esta, asegurándose que las temperaturas trabajadas sean menores al punto
de ebullición.

 La lectura del descenso de la acetona debe realizarse de tal modo que se fije primero la
altura del liquido en el vernier para que justo en el momento que llegue a dicha altura se
tome el tiempo, ya sea para el tiempo cero o para las próximas lecturas.

 Evitar apoyarse o mover sobre la mesa que sostiene el equipo de difusión pues alteraría
las condiciones de trabajo como la concentración dentro del capilar.

1
GRUPO E
8
 DIFUSION

BIBLIOGRAFIA

1. BASMADJIAN, DIRAN, “Mass Transfer: Principles and Applications”, pp.: 91-

95.

2. BENÍTEZ, JAIME; “Principles and Modern Applications of Mass Transfer

Operations”; Wiley; 2da edición; USA; 1948; pp.: 18-19.

3. MCCABE, WARREN, “Operaciones unitarias en la ingeniería química”,

McGraw-Hill, 4ta edición, pp.: 678-679.

4. PERRY, ROBERT; “Chemical Engineers Handbook”, McGraw-Hill; 8va

edición; pp: 2-43, 2-138, 2-141, 5-51, 5-52.

5. Manual del Equipo Armfields CERA-Diffusion of a gas apparatus.

6. YAWS, CARL L., “Yaws' handbook of thermodynamic and physical properties

of chemical compounds”, knovel, 2003

1
GRUPO E
9
 DIFUSION

APÉNDICE

EJEMPLO DE CÁLCULO
1. Calculo de la difusividad a partir de los datos experimentales a 40°C
(experimentales)
A partir de los datos experimentales obtenidos del sistema en el cual se trabajo, se
traza la gráfica t/(L-Lo) vs. (L-Lo) obteniéndose una tendencia lineal (ver grafica N°1.
Los términos de la ecuación de recta obtenida se reemplazan con los de la ecuación (17):
 L,A  C BM    L, A C BM 
L  L 0   
t
  L 0
L  L 0 2M A D AB  C A CT   M A D AB C A C T 
Las ecuaciones obtenidas a partir de los valores experimentales son (para la lectura Nº1):
t
 0.0205( L  Lo)  1.61
( L  Lo)

2. Cálculo de la Concentración Molar Total, CT

 1  Tabs
CT   
 Kmol.Vol  Ta
Kmol.Vol: Volumen Molar a condiciones normales = 22.4 m3/Kmol
Tabs = 273.15K
Ta = 313.15K
 1  273.15 K Kmol
CT     0.0389
 22.4m / Kmol  313.15 K
3
m3

3. Cálculo de la Presión de vapor de la Acetona

Haciendo uso de la ecuación de Antoine:
B
ln( PV )  A 
C  T(K )
Donde:
PV: Presión de vapor de la acetona
Los valores para A, B y C son aceptados únicamente para la acetona los cuales deben
estar dentro de:
241  T(K )  350

2
GRUPO E
0
 DIFUSION

Siendo estos:
 A= 16.6513
 B= 2940.46
 C= -35.93
Reemplazando en la ecuación de Antoine, a 40°C <> 313.15K:
2940.46
ln( PV )  16.6513 
 35.93  (313.15K )

PV  421.7mmHg  56.2kPa  0.555atm

4. Cálculo de la Concentración Molar del Aire en la Superficie del Líquido (CB2)
y en la Corriente de Aire (CB1)
Kmol
C B1  CT  0.0389
m3
Pa  PV
CB2   CT
Pa

Pa = Presión de atmosférica (atm).

PV = Presión de vapor de la acetona (atm).

C B2 
1  0.555atm  0.0389 Kmol  0.0173 Kmol
1atm m3 m3
5. Cálculo de la Concentración Media Logarítmica del Vapor de Acetona (CBM)
C B1  C B 2
C BM 
C
ln B1
CB2
0.0389  0.0173 kmol
C BM   0.0267 3
0.0389 m
ln
0.0173
6. Cálculo de la Concentración de Saturación del Acetona en la Interface (CA)
PV
CA   CT
Pa
0.555atm Kmol Kmol
CA   0.0389 3
 0.0216
1atm m m3

Luego de realizar los cálculos de todas las concentraciones necesarias para aplicar el
método de WINKELMANN, se procede a graficar los datos realizados
experimentalmente. Tomando en cuenta la tendencia lineal se obtiene la ecuación:

2
GRUPO E
1
 DIFUSION

t /( L  L0 )  0.0205( L  L0 )  1.61
Igualando está pendiente con la de ecuación (18), tenemos:

ks 1000s 10 6 mm 2  L, A  C BM 
0.0205     
mm 2 1ks 1m 2 2M A D AB  C A CT 
Donde:
L,A: Densidad de la Acetona líquida(40ºC) = 767.88 kg/m3
MA: Peso molecular de la Acetona = 58.08kg/kmol
Despejando la Difusividad, DAB y reemplazando valores:
767.88kg / m 3   m
3
0.0267 1
D AB    
2  58.08kg / kmol  0.0216  0.0389  Kmol 0.0205  10 9 s / m 2

m2 5
D AB  1.02 x10
s
𝐷𝑎𝑐𝑒𝑡𝑜𝑛𝑎 𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 ( 40º𝐶) = 0.102 𝑐𝑚2 /𝑠
7. Cálculo de la Difusividad de Acetona con los datos del manual
Tomando en cuenta la tendencia lineal del grafico N°4; se obtiene la ecuación:
t /( L  L0 )  0.0186( L  L0 )  1.638
Igualando está pendiente con la de ecuación (18), tenemos:

ks 1000s 10 6 mm2  L, A  C BM 
0.0186 2
  2
  
mm 1ks 1m 2M A D AB  C A CT 
8. Cálculo de la Concentración Molar Total, CT
 1  Tabs
CT   
 Kmol .Vol  Ta
Kmol.Vol: Volumen Molar a condiciones normales = 22.4 m3/Kmol
Tabs = 273.15K
Ta = 313.15K
 1  273.15K Kmol
CT     0.0389
 22.4m / Kmol  313.15K
3
m3

9. Cálculo de la Presión de vapor de la Acetona

Haciendo uso de la ecuación de Antoine:
B
ln( PV )  A 
C  T(K )
Donde:
PV: Presión de vapor de la acetona

2
GRUPO E
2
 DIFUSION

Los valores para A, B y C son aceptados únicamente para la acetona los cuales deben
estar dentro de: 241  T(K )  350
Siendo estos:
 A= 16.6513
 B= 2940.46
 C= -35.93
Reemplazando en la ecuación de Antoine, a 40°C <> 313.15K:
2940.46
ln( PV )  16.6513 
 35.93  (313.15K )
PV  421.6mmHg  56.21kPa  0.555atm
10. Cálculo de la Concentración Molar del Aire en la Superficie del Líquido (CB2)
y en la Corriente de Aire (CB1)
Kmol
CB1  CT  0.0389
m3
Pa  PV
CB 2   CT
Pa
Pa = Presión de atmosférica (atm).
PV = Presión de vapor de la acetona (atm).
CB2 
1  0.555atm  0.0389 Kmol  0.0173 Kmol
1atm m3 m3
11. Cálculo de la Concentración Media Logarítmica del Vapor de Acetona (CBM)
CB1  CB 2
CBM 
C
ln B1
CB 2
0.0389  0.174 kmol
CBM   0.0267 3
0.0389 m
ln
0.0174
12. Cálculo de la Concentración de Saturación del Acetona en la Interface (CA)
PV
CA   CT
Pa
0.553atm Kmol Kmol
CA   0.0389 3
 0.0215
1atm m m3

Igualando la pendiente con la de ecuación (18), tenemos:

ks 1000s 10 6 mm2  L, A  C BM 
0.0186 2
  2
  
mm 1ks 1m 2M A D AB  C A CT 
Donde:
L,A: Densidad de la Acetona líquida = 767.88kg/m3
MA: Peso molecular de la Acetona = 58.08kg/kmol
2
GRUPO E
3
 DIFUSION

Despejando la Difusividad, DAB y reemplazando valore en la ecuación (19):

767.88kg / m 3   m
3
0.0267 1
D AB    
2  58.08kg / kmol  0.0215  0.0389  Kmol 0.0186  10 9 s / m 2
m2
D AB  1.135 x10 5
s
𝐷𝑎𝑐𝑒𝑡𝑜𝑛𝑎− 𝑑𝑎𝑡𝑜𝑠 𝑚𝑎𝑛𝑢𝑎𝑙 ( 40º𝐶) = 0.11068𝑐𝑚2 /𝑠
13. Calculo del coeficiente de difusividad teórico, según el Yaws’Handbook:
𝐷𝑎𝑐𝑒𝑡𝑜𝑛𝑎 𝑡𝑒𝑜𝑟𝑖𝑐𝑜 = 6.7556𝑥10−7 𝑇 2 + 3.2888𝑥10−4 𝑥𝑇 − 0.0532
Esta relación es con respecto a la temperatura, en el handbook se encontró los
coeficientes para dicha ecuación. Reemplazando para T=40°C=313K.

𝐷𝑎𝑐𝑒𝑡𝑜𝑛𝑎 𝑡𝑒𝑜𝑟𝑖𝑐𝑜 = 6.7556𝑥10−7 (313)2 + 3.2888𝑥10−4 𝑥(313) − 0.0532

𝐷𝑎𝑐𝑒𝑡𝑜𝑛𝑎 𝑡𝑒𝑜𝑟𝑖𝑐𝑜 = 0.11591 𝑐𝑚2 /𝑠

14. Valores de difusividad para cada caso:

𝐷𝑎𝑐𝑒𝑡𝑜𝑛𝑎 𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 ( 40º𝐶) = 0.03172 𝑐𝑚2 /𝑠

𝐷𝑎𝑐𝑒𝑡𝑜𝑛𝑎 𝑡𝑒𝑜𝑟𝑖𝑐𝑜 ( 40º𝐶) = 0.11591 𝑐𝑚2 /𝑠
𝐷𝑎𝑐𝑒𝑡𝑜𝑛𝑎 𝑚𝑎𝑛𝑢𝑎𝑙 ( 40º𝐶) = 0.11068 𝑐𝑚2 /𝑠

2
GRUPO E
4
 DIFUSION

GRAFICAS

Grafica N°1: Datos del manual del equipo Armfield Limited (Acetona 40°C)

Datos del manual

3

2.5

2
t/(L-Lo)

1.5

0.5

0
0 5 10 15 20 25 30 35 40 45
L-Lo

Grafica N°2: Datos del manual del equipo Armfield Limited, ajustados de manera
lineal, por lo cual se elimina el primer punto.

Datos del Manual (lineal)

y = 0.0186x + 1.6378
3.000
R² = 0.9971
2.500

2.000
t/(L-Lo)

1.500

1.000

0.500

0.000
0 5 10 15 20 25 30 35 40 45
L-Lo

2
GRUPO E
5
 DIFUSION

Grafica N°3: Datos experimentales con el equipo Armfield Limited (Acetona 40°C)

Datos Experimentales
3

2.5

2
t/(L-Lo)

1.5

0.5

0
0 0.5 1 1.5 2 2.5 3 3.5 4
L-Lo

Grafica N°4: Datos experimentales con el equipo Armfield Limited, ajustados de

manera lineal, por lo cual se elimina los primeros puntos. (Acetona 40°C)

Datos Experimentales (lineal) y = 0.0649x + 1.8011

R² = 0.9957
2.04
2.02
2
t/(L-Lo) (ks/mm)

1.98
1.96
1.94
1.92
1.9
1.88
1.86
0 0.5 1 1.5 2 2.5 3 3.5 4
L-Lo (mm)

Series1 Linear (Series1)

2
GRUPO E
6