DIFUSIÓN MOLECULAR EN

FLUIDOS
 Ley de difusión de Fick
dz

CA,z CA,z+dz
A

C A x A  moles 
J A   DAB  cDAB JA   
z z  tiempo * Area 
 Agua Etanol
Kg kmol kg kmol
Inicialmente 100 5.55 100 2.17
 Agua Etanol
Kg kmol kg kmol
Inicialmente 100 5.55 100 2.17
Finalmente 44.08 2.45 55.92 3.10
Agua
Etanol 44.08 0.96 55.92 1.21
Total 88.16 3.41 111.84 4.31

Flujo de A Flujo total * x A ( Flujo de A  Flujo de B) * x A

 
Area Area Area
Flujo de A  Flujo de A Flujo de B 
   * xA   N A  N B  * xA
Area  Area Area 
 N  N A  NB
Flux de A resultante del movimiento total= N xA
El movimiento de A esta formado por dos partes:
• Movimiento de A resultante del movimiento global de materia.
• Movimiento de resultante de la difusión.

cA c A
N A  Nx A  J A  N A  ( N A  N B )  DAB
c z
cB cB
N B  NxB  J B  N B  ( N A  N B )  DBA
c z
cA c A cB cB
N A  N B  ( N A  N B )  DAB  ( N A  N B )  DBA
c z c z
(c A  c B ) c A cB
N A  NB  (N A  NB )  DAB  DBA
c z z
c A  cB  c
1

c c A cB
N A  N B  ( N A  N B )  DAB  DBA
c z z
c A cB
 DAB  DBA J A  JB
z z
c A  (c  c A )
 DAB  DBA
z z
Si c es constante (flujo equimolar):
c A   ( c )  (c A )   (c A )
 DAB  DBA      DBA DAB  DBA
z  z z  z
Difusión molecular en estado estacionario en
fluidos sin movimiento y en flujo laminar
Corriente de B

Gas B

Liquido
A
NA es constante
Difusión molecular en estado estacionario
en fluidos sin movimiento y en flujo
laminar
Si la ecuación 2.4 se aplica al caso de la difusión en el sentido de z únicamente, con
NA y NB constantes (estado estacionario) y considerando DAB constante:

cA c A
N A  Nx A  J A  N A  ( N A  N B )  DAB (2.4)
c z
c A cA 1
DAB  ( N A  N B )  N A   ( N A  N B )c A  cN A 
z c c
1 1
c A  z
( N A  N B )c A  cN A cDAB
CA z
1 1
C ( N A  N B )c A  cN A c A  z cDAB z (2.20)
A1 1
 CA z
1 1
ln ( N A  N B )c A  cN A   z
(N A  NB ) C
cDAB z
A1 1

1  ( N A  N B )c A  cN A  1
ln    ( z  z1 )
( N A  N B )  ( N A  N B )c A1  cN A  cDAB

Si multiplicamos la parte interna de logaritmo por 1 en la forma:

1  c( N A  N B )
1  c( N A  N B )

z  z1 1  (c A c )  N A ( N A  N B ) 
 ln 
cDAB ( N A  N B )  (c A1 c)  N A ( N A  N B ) 
Si z toma el valor límite de la parte superior z=z2
Y si Z es la distancia total del recorrido: Z= z2 – z1

Z 1  (c A 2 c )  N A ( N A  N B ) 
 ln  
cDAB ( N A  N B )  (c A1 c)  N A ( N A  N B ) 

Si multiplicamos ambas parte de la ecuación por:

N AcDAB
Z

NA cDAB  (c A2 c)  N A ( N A  N B ) 
NA  ln   (2.22)
(N A  NB ) Z (
 A1c c )  N A ( N A  N )
B 
 Para gases
Si se trata de gases y aplicamos la ley de gases ideales:

PT c A PA / RT PA
c n
V     yA
RT c PT / RT PT

NA PT DAB  y A 2  N A ( N A  N B ) 
NA  ln   (2.26)
( N A  N B ) RTZ  y A1  N A ( N A  N B ) 
Para utilizar estas ecuaciones, debe conocerse la relación entre NA y NB.
Esta relación generalmente se establece de acuerdo a condiciones que no
dependen del mecanismo de difusión.
 Hay tres clases de relación entre NA y NB:
• Difusión en estado estacionario de A a través del no difundente B.
• Contradifusión equimolar en estado estacionario.
• Contradifusión no equimolar en estado estacionario.

Difusión en estado estacionario de A a través del no difundente B

En este caso tenemos que NB=0 y NA es constante:

NA
1
N A  NB
Retomando la ecuación 2.22
NA cDAB  (c A 2 c)  N A ( N A  N B )  cDAB  (c A 2 c)  1
NA  ln    ln  
(N A  NB ) Z (
 A1c c )  N A ( N A  N )
B  Z (
 A1c c )  1 
cDAB 1  (c A 2 c)  cDAB  (c  c A 2 ) c cDAB  cB 2 
NA  ln    ln   ln  
Z  1  (c A1 c)  Z  (c  c A1 ) c  Z  cB1 
Entonces:
cDAB  cB 2 
NA  ln  
Z  cB1 
PT PB1 PB 2
Por ejemplo si B es un gas: c cB1  cB 2 
RT RT RT
PT DAB  PB 2 
NA  ln  
RT Z P
 B1 

PB 2  PT  PA2 PB1  PT  PA1

PB 2  PB1  ( PT  PA2 )  ( PT  PA1 )  PA1  PA2
PA1  PA2
1
PB 2  PB1  
 
PT DAB  PA1  PA 2   PB 2  PT DAB  ( PA1  PA 2 ) 
NA    ln   
RT Z  PB 2  PB1   PB1  RT Z  PB 2  PB1 
 ln P P  
 B2 B1 
 PB 2  PB1 PT DAB ( PA1  PA 2 )
PB , M  NA 
ln PB 2 PB1  RT Z PB , M
 Solución numérica
cA c A
N A  ( N A  N B )  DAB
c z
En estado estable NA= constante y NB=0
P c A PA / RT PA
c  Vn     yA
RT c P / RT P
PA DAB PA
NA  NA 
PT RT z
DAB PA PA
 NA  NA
RT z PT
PA RTN A  PA 
   1
z DAB  PT 
Método de Runge kutta:

k1  h f ( xi , yi ) k2  h f h k
( xi  , yi  1 )
2 2

k3  h f h k2 k 4  h f ( xi  h, yi  k3 )
( xi  , yi  )
2 2
1
yi 1  yi  (k1  2k 2  2k3  k 4 )
6
xi 1  xi  h
 Suponer un valor para NA

Integrar desde Z=0 hasta Z=0.002

Comparar PA calculado con PA,z2

No
PA,calculado  PA, z 2  Error

Si

Final de la iteración
 Difusión de B – no equimolar
En este caso la difusión de B se da en el sentido contrario y

N B  N A  1
Por ejemplo:

CH 4

CH 4  C  2H 2
2H 2
N B  2 N A
 Retomando la ecuación 2.26

NA PT DAB  y A 2  N A ( N A  N B ) 
NA  ln 
( N A  N B ) RTZ  y A1  N A ( N A  N B ) 

NA NA NA 1
  
( N A  N B ) ( N A  N A ) N A (1   ) (1   )
 1 
y 
1 PT DAB  A 2 (1   ) 
NA  ln  
(1   ) RTZ  y  1 
 A1
(1   ) 

 y A 2 (1   )  1 
1 PT DAB  (1   )  1 PT DAB  y A 2 (1   )  1
NA  ln    ln 
(1   ) RTZ  y A1 (1   )  1  (1   ) RTZ  y A1 (1   )  1 
 (1   ) 
 1 PT DAB  PT (1   )  1
PA 2

NA  ln  PA1 
(1   ) RTZ  PT (1   )  1 

1 PT DAB  PA 2 (1   )  PT 
NA  ln 
(1   ) RTZ  PA1 (1   )  PT 

CH 4  C  2H 2

N B  2 N A   2

1 PT DAB   PA 2  PT  PT DAB  PT  PA 2 
NA  ln    ln 
 1 RTZ   PA1  PT  RTZ  PT  PA1 
 Contradifusión equimolar
En este caso: NB  N A
N A  NB  0
No podemos aplicar la ecuación 2.22

NA cDAB  (c A 2 c)  N A ( N A  N B ) 
NA  ln   (2.22)
(N A  NB ) Z  (c A1 c)  N A ( N A  N B ) 
cA c A
N A  ( N A  N B )  DAB
c z
c A NA
N A   DAB c A   z
z DAB
CA z
NA NA
C c A   DAB Z z C A  C A1  
DAB
( z  z1 )
A1 1

(C A  C A1 )
N A   DAB
( z  z1 )
 DAB
Si z=z2 y CA=CA2, con Z=z2-z1 NA  (C A 2  C A1 )
Z

Para el caso de gases, considerando comportamiento de gases ideales:

 DAB
NA  ( PA2  PA1 )
RTZ