Laboratorio de Ingeniería Química II 0

Difusión de acetona en el aire

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

“DETERMINACIÓN DEL COEFICIENTE DE DIFUSIÓN DE LA ACETONA EN AIRE A 40°C”

PROFESOR : Mg. Cesario Condorhuamán Ccorimanya

ALUMNOS : Alvarado Huanca Manuel Augusto

Martínez Meza Diego Antonio 12070190
Núñez Aguilar Giuseppe Luigi 12070154
Yagua Farfán Emily Lucero 13070198

FECHA DE
REALIZACIÓN : 5 de mayo del 2017

FECHA DE
ENTREGA : 5 de junio del 2017
 Laboratorio de Ingeniería Química II 1
Difusión de acetona en el aire

Ciudad Universitaria, junio del 2017

RESUMEN
En el presente informe se quiere determinar el coeficiente de difusividad de la acetona en aire,
en la que se utilizó el “equipo de difusión de gases Armfield” a una temperatura de operación
de 40°C. La práctica se desarrolla a las condiciones ambientales de 20 °C y a 756mmHg.

Se emplea el Método de Winklemann y resulta un coeficiente de difusividad para la acetona

en aire a una temperatura de 40°C de 2.8719 x 10-6 m2/s que, comparado con la difusividad
de la acetona en aire a 40°C proporcionado por datos experimentales del manual, con un
valor de 12.4 x 10-6 m2/s, se obtiene un % de desviación de 76%.
 Laboratorio de Ingeniería Química II 2
Difusión de acetona en el aire

INTRODUCCION
El movimiento de una especie química desde una región de alta concentración hacia una
región de baja concentración puede observarse al dejar caer un sustancia química A a través de una
mezcla binaria A y B debido a una gradiente de concentración de A.
Este fenómeno se fundamenta en la ley de Fick de la difusión. Cuando la gradiente se
mantiene mediante el suministro continuo de los componentes 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. Por ejemplo cuando se remueve amoniaco de un gas por
medio de absorción en agua en una columna empacada, en cada punto de la columna una gradiente
de concentración en la fase gaseosa provoca la difusión del amoniaco a la interfase gas-líquido , donde
se disuelve , y un gradiente en la fase líquida provoca la difusión en la masa líquida. Cuando se extrae
un soluto de un líquido, los gradientes se invierten; aquí la difusión conduce al soluto desde la masa
líquida a la interfase y de ahí hacia el interior de la fase gaseosa.
En algunas otras operaciones de transferencia de masa, tales como lixiviación y la adsorción,
tiene lugar la difusión en estado no estacionario, y los gradientes y flujos disminuyen con el tiempo en
cuanto se alcanza el equilibrio. En la destilación, el componente menos volátil se difunde a través de
la fase líquida hacia la interfase y desde esta hacia el vapor. En la lixiviación, la difusión del soluto a
través de la fase sólida va seguida de la difusión dentro del líquido. En la extracción de líquido, el
soluto se difunde a través de la fase de refinado hacia la interfase y después dentro de la fase
extracto. Así podemos observar que la difusión cumple un rol importante en las operaciones de
transferencia de masa.
Aunque la causa habitual de la difusión es un gradiente de concentración, en ocasiones la
difusión también se origina por un gradiente de actividad, como la ósmosis inversa, por un gradiente
de presión, un gradiente de temperatura o por la aplicación de un campo de fuerza externa como en
el caso de la centrífuga. La difusión molecular inducida por la temperatura es la difusión térmica, y
debido a un campo externo es la difusión forzada. Sin embargo en este informe solo se estudiará la
difusión por causa de una gradiente de concentración.
 Laboratorio de Ingeniería Química II 3
Difusión de acetona en el aire

PRINCIPIOS TEORICOS

1. FUNDAMENTO DE LA DIFUSIÓN
La difusión es el movimiento, bajo la influencia de un estimulo 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.

2. PRIMERA LEY DE DIFUSIÓN DE FICK

En la siguiente ecuación general para la difusión unidimensional:
dC A
J A   D AB
dz
Donde:
JA=flux molar del componente A, Kg mol/m2h o lb mol/f2 h
DAB= coeficiente de difusividad del componente A en B
dCA/dz =gradiente de concentración del componente A en la dirección de la difusión
Esta ley se basa en tres decisiones:
1. El flux está en moles por unidad de área por tiempo.
2. La velocidad de difusión es relativa a la velocidad volumétrica promedio.
3. El potencial impulsor está en términos de concentración molar ( moles del componente A
por unidad de volumen).

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
 Laboratorio de Ingeniería Química II 4
Difusión de acetona en el aire

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  cD AB (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 )  cD AB (2)
z
Flux molar Flux molar que Flux de difusión
relativo a una resulta del flujo relativo a la veloc.
superficie fija global molar promedio

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.
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
 Laboratorio de Ingeniería Química II 5
Difusión de acetona en el aire

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. DEPENDENCIA DE LAS DIFUSIVIDADES CON RESPECTO A LA TEMPERATURA Y LA PRESIÓN

Para mezclas gaseosas binarias a baja presión, D AB es inversamente proporcional a la presión,
aumenta con la temperatura y es casi independiente de la composición para una pareja de gases
dada. La siguiente ecuación para estimar D AB a bajas presiones se dedujo a partir de una combinación
de argumentos de la teoría cinética y de los estados correspondientes.:
pD AB T
 a( )b
( p cA p cB ) 1/ 3
(TcATcB ) 5 / 12
(1 / M A  1 / M B ) 1/ 2
TCATCB

Donde:
DAB = cm2/s
P=atm
T=K
A= 3.640x10^-4 y b=2.334 , adimensionales. Para parejas de que constan de H 2O y un gas no polar.
A altas presiones y en estado líquido, el comportamiento de DAB es más complicado. La situación más
simple y que se comprende mejor es la de autodifusión ( interdifusión de moléculas de las mismas
especies químicas). Pero no se abordará este tema por no ser objeto de estudio.
1
Lo mejor es disponer de valores experimentales para estimar las difusividades; y cuando esta
disponible tal información para el sistema de interés deberá utilizarse directamente. Sin embargo, con
 Laboratorio de Ingeniería Química II 6
Difusión de acetona en el aire

frecuencia no se dispone de los valores deseados y es preciso estimarlos a partir de correlaciones

publicadas. A veces e dispone de un valor para un conjunto de condiciones de temperatura y presión;
entonces las correlaciones resultan útiles para predecir, a partir de un valor conocido, los valores
deseados en otras condiciones.
Una sencilla teoría para los gases muestra que la difusividad es proporcional al producto de la
velocidad molecular media u y la trayectoria libre media λ.
1
Dv  u.
3
Algunas ecuaciones empíricas predicen la difusividad con la variación de la temperatura,
considerándose que la presión no afectaría hasta un valor de 10 atm. Pero va a depender mucho de
las consideraciones de trabajo para tomar en cuenta ecuaciones empíricas.

4. METODO DE WINKELMANN
Consideremos el sistema de difusión que se presenta en la fig. N° 1. El liquido A se está
evaporando en el seno del gas B, e imaginemos que mediante un artificio es posible mantener el nivel
del liquido 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 una mezcla gaseosa ideal. Finalmente se supone que la
solubilidad de B en el liquido 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 x A2, 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:
cD AB 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)
 Laboratorio de Ingeniería Química II 7
Difusión de acetona en el aire

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 N Az.
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 
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 A 2 ) (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  CT 
N A  D A   (11)
 L  C BM 
Donde:
D: Difusividad (m2/s).
CA: Concentración de saturación en la interface (Kmol/m 3)
 Laboratorio de Ingeniería Química II 8
Difusión de acetona en el aire

L: Distancia efectiva de transferencia de masa (mm)

CBM: Media logarítmica de concentración de vapor (Kmol/m 3)
CT: Concentración molar total = C A  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  xx A 2 ) S (12)
M A dt ( z 2  z1 (t )) x BM

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 ecuacion (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 A CT t
Lo LdL    L  C BM   tdt
0
(14)

 2 MD  C ACT 
L2  L20    t (15)
  L  C BM 
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  2 Lo)( L  Lo)

 Laboratorio de Ingeniería Química II 9
Difusión de acetona en el aire

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

t    C    C 
  L  BM  L  L0    L BM  L0 (17)
 L  L0   2MD  C ACT   MDC A CT 
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 C T 
Donde:

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

*Kmol x Volumen = 22.414 m3/Kmol



 C B1  C T

 P  PV 
 C B 2   a 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.
 Laboratorio de Ingeniería Química II 10
Difusión de acetona en el aire

(a)

Aire
Referencia
z=z2

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

∆L=∆z

L = L; t = t

NA|L

Acetona
 Laboratorio de Ingeniería Química II 11
Difusión de acetona en el aire

(b)

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

disminuye muy lentamente cuando se evapora.
 Laboratorio de Ingeniería Química II 12
Difusión de acetona en el aire

DETALLES DEL EQUIPO Y PROCEDIMIENTO EXPERIMENTAL

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.
 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.
 Microscopio.
 Vernier (con una lectura de hasta dos decimales).
 Controlador de temperatura digital.

Fig. N°3: Esquema del equipo

 Laboratorio de Ingeniería Química II 13
Difusión de acetona en el aire

PROCEDIMIENTO EXPERIMENTAL

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

bomba de aire 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.
 Laboratorio de Ingeniería Química II 14
Difusión de acetona en el aire

TABULACION DE DATOS Y RESULTADOS

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

Tabla N° 2: Datos Experimentales de la práctica .Mediciones de la Altura del Capilar

Tiempo (t) t (seg) t (ks) L (mm) Lo (mm) L-Lo (mm) t/(L-Lo)

(ks/mm)
00:00:00 0 0 34.6 34.6 0 -
2:57:32 177.32 0.17732 34.7 34.6 0.1 1.7732
16:43:28 1003.28 1.00328 35 34.6 0.4 2.5082
24:53:34 1493.34 1.49334 35.3 34.6 0.7 2.133342857
34:47:84 2087.84 2.08784 35.5 34.6 0.9 2.319822222
44:50:26 2690.26 2.69026 35.9 34.6 1.3 2.069430769
54:56:47 3296.47 3.29647 36.2 34.6 1.6 2.06029375
65:06:17 3906.17 3.90617 36.5 34.6 1.9 2.055878947
85:00:92 5100.92 5.10092 37.3 34.6 2.7 1.88922963
125:28:80 7528.8 7.5288 38.4 34.6 3.8 1.981263158
145:25:82 8725.82 8.72582 38.9 34.6 4.3 2.029260465
165:04:14 9904.14 9.90414 39.4 34.6 4.8 2.0633625
185:10:24 11110.24 11.11024 39.9 34.6 5.3 2.096271698
204:52:92 12292.92 12.29292 40.4 34.6 5.8 2.119468966

Tabla N°3: Datos experimentales proporcionados por el 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
 Laboratorio de Ingeniería Química II 15
Difusión de acetona en el aire

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°4: Concentraciones obtenidas de la práctica para calcular el coeficiente de difusión

temperatura 40
CB1 (kmol/m3) 0.0389
CB2 (kmol/m3) 0.0174
CT (kmol/m3) 0.0389
CBM (kmol/m3) 0.0267
CA (kmol/m3) 0.0215

Tabla N° 5: Coeficientes de Difusividad de acetona-aire calculados por el Método de Winkleman

Coeficiente de difusividad DAB(m2/s) (40ºC) % Desviación

Experimento 2.8719 x 10-6
76
Manual 12.4 x 10-6
 Laboratorio de Ingeniería Química II 16
Difusión de acetona en el aire

I. ANALISIS Y DISCUSIÓN DE RESULTADOS

1. Para hallar la pendiente de la ecuación WINKLEMANN, solo se tomaron los puntos que
representan el estado estacionario a una temperatura de 40°C, por lo que fueron omitidos los
puntos en estado no estacionario. Del punto de máximo valor de t/(L-Lo) en la gráfica, podemos
deducir como el punto de mínima difusión a una menor temperatura.

2. El valor del coeficiente de difusión hallado experimentalmente difiere del valor hallado por el
manual y esto pudo deberse al poco tiempo en el que se trabajó.

3. Se observó que a medida que pasaba el tiempo las alturas del menisco del capilar que contenía
acetona aumentaba, esto quiere decir que el visor del microscopio estaba invertido.
 Laboratorio de Ingeniería Química II 17
Difusión de acetona en el aire

II. CONCLUSIONES

1. El valor del coeficiente de difusividad para la acetona en aire a una temperatura de 40°C es
2.8719 x 10-6 m2/s.

2. El porcentaje de desviación entre el valor experimental y el valor hallado con los datos del
manual resulta 76%.

III. RECOMENDACIONES

1. Realizar la mayor cantidad de lecturas, pero a largos intervalos de tiempo, esto nos permite
observar mejor la linealidad de la curva “L-L0 vs t/(L-L0)”.
2. Realizar la experiencia a varias temperaturas para observar el comportamiento frente a la
variación de esta, pero a temperaturas menores al punto de ebullición.
3. 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.
4. 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.

IV. BIBLIOGRAFIA
1. BASMADJIAN, DIRAN, “Mass Transfer: Principles and Applications”, pp.: 91-95.
 Laboratorio de Ingeniería Química II 18
Difusión de acetona en el aire

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.

 Laboratorio de Ingeniería Química II 19
Difusión de acetona en el aire

APÉNDICE
EJEMPLO DE CÁLCULO
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):

t  L,A  C BM    L , A C BM 
   L  L 0    L 0
L  L 0 2M A D AB  CACT   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.0756( L  Lo)  1.6933
( L  Lo)

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

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

 1  273.15 K Kmol
CT   3   0.0389
 22.4m / Kmol  313.15 K m3

2. 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

Siendo estos:
 A= 16.6513
 Laboratorio de Ingeniería Química II 20
Difusión de acetona en el aire

 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.7 mmHg  56.2kPa  0.555atm

3. 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  C T  0.0389
m3
Pa  PV
CB2   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.0174
Kmol
3
1atm m m3
4. 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
0.0389 m3
ln
0.0173
5. 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.0215
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:
t /( L  L0 )  0.0756( L  L0 )  1.6933
 Laboratorio de Ingeniería Química II 21
Difusión de acetona en el aire

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

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

m2
D A exp  2.8719 x10  6
s

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.0175( L  L0 )  1.637

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

ks 1000s 10 6 mm 2  L, A  C BM 
0.0175     
mm 2 1ks 1m 2 2 M A D AB  C A CT 
1. Cálculo de la Concentración Molar Total, CT

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

2. 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:
 Laboratorio de Ingeniería Química II 22
Difusión de acetona en el aire

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
3. 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
4. Cálculo de la Concentración Media Logarítmica del Vapor de Acetona (CBM)
CB 1  C B 2
CBM 
C
ln B1
CB 2
0.0389  0.174 kmol
CBM   0.0267
0.0389 m3
ln
0.0174

5. 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 mm 2  L, A  C BM 
0.0175     
mm 2 1ks 1m 2 2M A D AB  C A CT 
 Laboratorio de Ingeniería Química II 23
Difusión de acetona en el aire

Donde:
L,A: Densidad de la Acetona líquida = 767.88kg/m 3
M A: Peso molecular de la Acetona = 58.08kg/kmol

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

790kg / m 3  0.0267  m
3
1
D AB    
2  58.08kg / kmol  0.0215  0.0389  Kmol 0.0175  10 9 s / m 2
m2
D Amanual  12.4 x10  6
s

Calculo del porcentaje de desviación:

DAmanual  DA exp
% DESVIACION   100%
DAmanual

12.4 x10 6  2.8719 x10 6

% DESVIACION   100%
12.4 x10 6

%DESVIACION = 76%
 Laboratorio de Ingeniería Química II 24
Difusión de acetona en el aire

GRÁFICAS
Gráfica Nº 1: representación de los datos experimentales, con la utilización de acetona a 40ºC

Gráfica Nº 2: representación de los datos experimentales que siguen una tendencia lineal, con
la utilización de acetona a 40ºC

Gráfica Nº 3: representación de los datos experimentales proporcionados por el manual, con

la utilización de acetona a 40ºC
 Laboratorio de Ingeniería Química II 25
Difusión de acetona en el aire

Gráfica Nº 4: comparación entre grafica obtenida mediante data experimental y grafica

obtenida con datos del manual.