CAPÍTULO I.

DIFUSIÓN EQUIMOLAR

Autor: María Gabriela Gómez C.

En un sistema binario, la difusión de A ocurre por gradientes de concentración de A.

Este fenómeno se llama difusión ordinaria, para distinguirlo de difusión por presión
(movimiento de A por un gradiente de presión), difusión térmica (movimiento de A
originada por un gradiente térmico) y difusión forzada (movimiento de A originado por
fuerzas externas diferentes sobre A y B).

I.1. Definiciones de concentraciones, velocidades y flujos de masa

Las concentraciones de varias especies en un sistema multicomponentes, pueden ser
expresadas de cuatro formas diferentes:
• Concentración másica: ρi (masa de la especie i por unidad de volumen de
solución).
• Concentración molar: ci (moles de la especie i por unidad de volumen de
solución)
• Fracción másica wi (concentración másica de la especie i dividida por la densidad
de la solución).
• Fracción molar xi (concentración molar de la especie i dividida por la densidad
molar de la solución).
En una mezcla de varias especies químicas que están difundiendo, ellas (las especies)
se mueven a diferentes velocidades. Si se llama ui a la velocidad de la especie i con
respecto a un sistema de coordenadas fijo (estático), se puede definir la velocidad
promedio local de masa como:

∑ρ u i i n
ρi ui n
u =
m i =1
=∑ = ∑ wi ui (1.1)
1=1 ρ
n

∑ρ
i =1
i
i =1

donde ρu corresponde a la velocidad local a la cual la masa pasa a través de una sección
transversal, perpendicular a la velocidad u. La velocidad local puede ser medida por
ejemplo con un tubo Pitot.
De forma similar, se puede definir la velocidad promedio local molar como:

∑c u i i n
ci ui n
u =
M i =1
n
=∑ = ∑ xi ui (1.2.)
c
∑c i =1
i
1=1 i =1

y también se puede definir la velocidad promedio local volumétrica como:

 uv = ∑ Vi ci ui (1.3)

donde Vi es el volumen molar parcial de la especie i.

En sistemas en flujo (donde los fluidos “fluyen”), a veces conviene definir la
velocidad de la especie i con respecto a un sistema de coordenadas fijo. De esta forma se
definen las denominadas velocidades de difusión de las siguientes formas:
ui − um : velocidad de difusión de la especie i con respecto a la velocidad promedio um
ui − u M : velocidad de difusión de la especie i con respecto a la velocidad promedio u M
Por otra parte, el flujo másico o molar de la especie i es un vector cuantitativo que
denota la masa o los moles de la especie i que pasa a través de un área unitaria por unidad
de tiempo.
Los flujos másico y molar se pueden definir en función del sistema de referencia
utilizado como se muestra a continuación:
• Relativo a un sistema de coordenadas fijo (estático):

ni = ρi ui  M tL2  (flujo másico) (1.4)

N i = ci ui  mol tL2  (flujo molar) (1.5)

• Relativo a la velocidad promedio másica local:

m i j = ρi ( ui − u m )  M tL2  (flujo másico) (1.6)

m J i = ci ( ui − u m )  mol tL2  (flujo molar) (1.7)

• Relativo a la velocidad promedio molar local:

M ji = ρi ( ui − u M )  M tL2  (flujo másico) (1.8)

M J i = ci ( ui − u M )  mol tL2  (flujo molar) (1.9)

La expresión (1.9) puede ser reescrita:

M J i = ci ( ui − u M ) = ci ui − ci u M = N i − ci u M ⇒ N i = M J i + ci u M (1.10)

La ecuación anterior se le conoce como la ecuación de balance de masa por flujo

global de mezcla. En esta ecuación se pueden distinguir dos contribuciones o razones por
las cuales se produce la transferencia de masa, a saber:
 Ni = M J i + ci u M (1.11)
 Flujo molar de la   Flujo difusivo de la   Flujo convectivo de la 
     
 especie i = especie i + especie i 
 por unidad de area   por unidad de area   por unidad de area 
     

El término convectivo se puede escribir también de la siguiente forma:

n n
c c n n
ci u M = ci ∑ xi ui = ci ∑ i ui = i ∑ ci ui = xi ∑ N i = xi N T (1.12)
i =1 i =1 c c i =1 i =1

donde N T es el flujo molar total por unidad de área de flujo.

Sustituyendo la Ec. (1.11) en (1.10) resulta:

n
Ni = M J i + xi ∑ N j (1.13)
j =1

I.2. Concepto de difusión

La difusión molecular o transporte molecular es la transferencia de masa o
movimiento individual de las moléculas a través de un fluido por medio de movimientos
individuales y aleatorios de las mismas. Los flujos difusivos aparecen cuando hay un
sistema de referencia que se desplaza a una velocidad de referencia promedio u m , u M , uv ;
es decir, se pueden definir expresiones como las representadas en las ecuaciones
comprendidas entre la (1.6) y la (1.9).
Con el fin de definir la contribución difusiva en el flujo neto de un componente, se
requiere de una ecuación constitutiva. En el transporte de cantidad de movimiento, la
ecuación constitutiva que se define regularmente es la ecuación de Newton, mientras que
en el transporte de calor, se define la ley de Fourier a saber:

 2
τ = −2µΓ −  κ −  ( ∇ ⋅ u ) δ ← ley de Newton
 3 (1.14)
q = − k ∇T ← ley de Fourier

donde κ es la viscosidad dilatacional, Γ es el tensor de tasa de deformación (es un tensor

simétrico), δ el delta de Kronecker el cual tiene como características:

1 , i = j
δ ij = 
0 , i ≠ j

La forma del tensor de tasa de deformación se presenta en el anexo A.

Los efectos de la viscosidad dilatacional, en la fluidodinámica son difíciles de
detectar y usualmente se desprecian. Se ha demostrado que esta propiedad es
exactamente igual a cero en el caso de gases monoatómicos ideales, mientras que se cree
que para otros fluidos κ << µ . Adicionalmente es claro que para el caso de fluidos
incompresibles, donde se cumple que ( ∇ ⋅ u ) = 0 , ese parámetro desaparece al eliminarse
el término que lo contiene en la ecuación de Newton. Por ello, la ecuación de Newton se
puede simplificar a:

 1 
τ = −2 µΓ − ( ∇ ⋅ u ) δ  (1.15)
 3 

Tomando en consideración la analogía existente entre los fenómenos de transporte, Fick

(1855) propuso una expresión, restringida a sistemas isotérmicos o isobáricos:

m i j = − D∇ρi (1.16)

donde D es el coeficiente de difusividad. El signo negativo de la ecuación anterior se

define por convención al igual que en las expresiones mostradas en (1.14), donde se
quiere que el tensor de esfuerzos, el flujo de calor o el flujo difusivo de masa sean
positivos, ya que el término de gradiente que se define respectivamente en cada uno, es
decreciente (posee un valor negativo).
La Ec. (1.16) fue reescrita por Groot (1951) quedando de la siguiente forma:

m i j = −ρD∇wi (1.17)

En el caso donde la densidad es constante, las Ecs. (1.16) y (1.17) son equivalentes.
En el caso cuando la densidad no es constante, la ecuación de Fick no es correcta, por lo
que se propuso la ecuación (1.17). La ecuación (1.17) también puede ser escrita en
términos de flujo molar:

M J i = −cD∇xi (1.18)

Sustituyendo entonces la ecuación constitutiva de flujo difusivo en la Ec. (1.13) se

obtiene:

n
N i = −cD∇xi + xi ∑ N i (1.19)
i =1

Con el objetivo de poder pasar de un sistema de referencia a otro, Hooyman et al.

(1953) propusieron la siguiente expresión:

1 (1 − χi )
Ji = − D ∇xi (1.20)
V (1 − xi )
donde xi es la fracción molar, χi es la fracción de peso reportada en la tabla 1.1 y V es el
volumen molar de la mezcla.

Tabla 1.1. Fracción de peso usada en la ecuación Hooyman et al. (1953)

Velocidad de referencia χi
um wi
uM xi
uv ciVi

I.3. Cálculo del coeficiente de difusividad

El coeficiente de difusividad es función de la temperatura y presión tal como lo es la
viscosidad y la conductividad térmica que son las propiedades usadas en cantidad de
movimiento y transferencia de calor respectivamente.
Para una mezcla binaria, gaseosa a baja presión se encuentra que:

Si p ↓ => DAB ↑

Si T ↑ => DAB ↑
Si xi ↑ o xi ↓ => DAB ≈ constante

Para el cálculo del coeficiente de difusividad se han propuesto diversas correlaciones

dependiendo de la fase, del intervalo de presión y/o temperatura en la que se encuentre el
sistema.

• Cálculo del coeficiente de difusividad para gases:

En una mezcla a baja presión, se desarrolló una expresión producto de combinar la
teoría cinética de los gases y los correspondientes argumentos de estado:

b
pDAB  T 
= a  (1.21)
 1
12
1   T T 
 cA cB 
( pcA pcB ) (TcATcB )
13 5 12
 + 
 MA MB 
donde:
DAB : coeficiente de difusividad de A en B [cm2/s]

pci : presión crítica del componente i [atm]

Tci : temperatura crítica del componente i [K]

Mi : peso molecular del componente i [g/mol]

p: presión del sistema [atm]

T: temperatura del sistema [K]
Para un par de gases no polares: a = 2.745×10-4 b = 1.823
-4
Para agua con un gas no polar: a = 3.640×10 b = 2.334
Si los gases son no polares y se dispone de los parámetros de Lennard-Jones (1962),
el método de la teoría cinética da mejores resultados (Ec.1.21).
Para un sistema gaseoso a alta presión (el coeficiente de difusión prácticamente no
cambia con la presión). Así, para un sistema gaseoso a baja densidad:

 1 1 
T + 
 MA MB 
cDAB = 2, 2646 ×10−5 (1.22)
σ2AB ΩD

donde si se aproxima c por gas ideal ( c = p RT ) :

 1 1 
T3  + 
 MA MB 
cDAB = 1,858 ×10−27 (1.23)
p σ2AB ΩD
donde:
DAB : coeficiente de difusividad [m2/s]
T: temperatura [K]
Mi: peso molecular [kg/kgmol]
p: presión absoluta [atm]
σAB: diámetro de colisión [m]
ΩD: integral de colisión [-]
En un sistema binario, las constantes de fuerza de Lennard-Jones se evalúan usando
las expresiones:
σ A + σB
σ AB = (1.24)
2
 ε AB ε ε
= A B (1.25)
k k k

donde k es la constante de Boltzman (1,3806×10-23 J K-1).

El procedimiento para evaluar la difusividad es:
- Se lee de la tabla 2-3 (pág. 21 del Hines & Maddox) los valores de σi y εi k para
cada componente puro.
- Se evalúan las expresiones (1.24) y (1.25).
- Se halla el inverso del valor obtenido de la ecuación (1.25) y luego este número se
multiplica por la temperatura del sistema en grados Kelvin; luego, con este último
valor se lee de la tabla 2-4 (pág. 23 del Hines & Maddox) el valor de Ω D .
- Se sustituyen todos los valores, ahora ya conocidos, en las ecuaciones (1.22) o
(1.23).
En el caso, en el que los componentes a trabajar no aparezcan en la tabla 2-3, se
utilizan las siguientes correlaciones:

σi = 0,833 Vci1 3 , σi = 1,18 Vbi1 3 (1.26)

εi εi
= 0, 77Tci , = 1,15Tbi (1.27)
k k
donde:
Vci : volumen crítico de la especie i [m3/gmol]
Tci : temperatura crítica de la especie i [K]
Vbi: volumen molar en el punto normal de ebullición de la especie i [m3/gmol]
Tbi : temperatura del punto de ebullición normal de la especie i [K]
σi: diámetro de colisión de la especie i, [10-10 m]
Hirschfelder et al. (1949) recomendaron para extrapolar los datos experimentales para
rangos de presión moderados:

32
p1  T2  ΩT1
DAB ,T2 , p2 = DAB ,T1 , p1   (1.28)
p2  T1  ΩT2
• Coeficiente de Difusividad para Líquidos:
Para el cálculo del coeficiente de difusión en líquido se han propuesto una serie de
correlaciones, a saber:
- Ecuación de Stokes-Einstein (Frenkel, 1946): supone que las partículas son esferas
duras que se mueven a una velocidad uniforme en un continuo, bajo una fuerza
unitaria F.
kT
D= (1.29)
6πrµ

donde:
D: coeficiente de autodifusión
k: constante de Boltzman (1.3806×10-23 J K-1).
T: temperatura absoluta
r: radio de la partícula que difunde
µ: viscosidad dinámica
Esta ecuación predice, en varios casos, los coeficientes de difusión con el orden de
magnitud correcto.
- Correlación de Sutherland-Einstein (1905): introduce una mejora en la ecuación
anterior, donde modifica empíricamente, la ley de Stokes.

kT  βr + 3µ 
D=   (1.30)
6πrµ  βr + 2µ 

donde β es el coeficiente de fricción por deslizamiento entre el átomo y el medio.

Para soluciones no electrolíticas diluidas, Wilke y Chang (1955) propusieron la
siguiente ecuación:

1,17 ×10−13 ( ξ B M B ) T
12

D 0
AB = (1.31)
VA0,6µ
donde:
0
DAB : coeficiente de interdifusión en soluciones diluidas [m2 / s]

µ: viscosidad de la solución [mPa·s]

VA: volumen molar del soluto en el punto de ebullición normal [m3 / kgmol]
MB: peso molecular del solvente [kg / kgmol]
T: temperatura absoluta [K]
ξB: factor de asociación del solvente B (tiene valores de: 1,0 para solventes no
polares como el benceno, éter o hidrocarburos alifáticos, 1,5 para el etanol, 1,9 para el
metanol y 2,6 para el agua.
- Sitaraman et al. (1963) propusieron una ecuación menos restrictiva que la de Wilke
y Chang, la cual se recomienda cuando el agua es el soluto:

0,93
−14  M B1 2 ∆H B1 3T 
D 0
AB = 16, 79 × 10  0,3 
(1.32)
 µ BVA ∆H A 
12

donde:
∆Hi: calores latentes de vaporización del soluto y solvente en sus puntos normales
de ebullición [J/kg]
µB: viscosidad del solvente en centipoises [mPa·s].
Para soluciones binarias no electrolíticas concentradas, donde la difusividad depende
de la concentración:

0  ∂ log γ A 
DAB = DAB 1 +  (1.33)
 ∂ log cA 

donde:
0
DAB : coeficiente de difusión independiente de la concentración

γA: coeficiente de actividad

- Vignes (1966) propuso una ecuación empírica (a partir de la ecuación anterior) que
ajusta para soluciones ideales y no ideales:
  ∂ log γ A 
DAB = ( DAB ) (D )
xB xA
1 +
0 0
BA  (1.34)
 ∂ log xA 

- Leffler y Cullinan (1970) modificaron la ecuación anterior para ajustar mejor los
datos experimentales, obteniendo:

xA  ∂ log γ A 
DAB µ m = ( DAB µB ) (D µ A ) 1 +
0 xB 0
BA  (1.35)
 ∂ log xA 

- Para soluciones electrolíticas diluidas, Nernst (1888) propuso (a dilución infinita):

RT λ 0+ λ 0− Z − + Z + −14 λ 0+ λ 0− Z − + Z +
0
DAB = = 8,931 × 10 T (1.36)
F 2 λ 0+ + λ 0− Z − Z + λ 0+ + λ 0− Z − Z +

donde:
0
DAB : coeficiente de difusión a dilución infinita [m2/s]
F: constante de Faraday [A s / g-equiv]
λ 0+ : conductancia catiónica a dilución infinita [(A/cm2)(cm/V)(cm3/g-equiv]

λ 0− : conductancia aniónica a dilución infinita [(A/cm2)(cm/V)(cm3/g-equiv]

λ 0+ + λ 0− : conductancia electrolítica a dilución infinita [(A/cm2)(cm/V)(cm3/g-equiv]

Z+: valencia catiónica
Z-: valencia aniónica
T: temperatura absoluta [K]
Los valores de la conductancia se presentan en la tabla 2-7 (pág. 35 Hines y Maddox).
- Para metales líquidos:
- Ecuación de Walls y Upthegrove (1964). Los autores modificaron la ecuación
de Stokes-Einstein. La ecuación es la siguiente:

kT
DAA = (1.37)
2πr ( 2b + 1) µ
donde:
DAA: coeficiente de autodifusión
b = 0.419
r: radio iónico univalente de Pauling
k: constante de Boltzman
- Ecuación de Hines et al. (1975):

k 2T 2 12 − Z
DAA = (1.38)
hK Z

donde:
K: constante de fuerza de Waser y Pauling (1950)
Z: sus valores se presentan en la tabla 2-8 (pág. 37 del Hines y Maddox) para
metales líquidos.
h: constante de Planck (6,6262×10-34 J·s)
- Ecuación de Hines y Walls (1979):

12
1  σ  T µ m 12 − Z
DAA = j 2 υs  L  (1.39)
6  σ S  Tm µ 2

donde:
T: temperatura absoluta [K]
Tm: temperatura del punto de fusión [K]
j: distancia de fluctuación media [m]
µ: viscosidad absoluta [Pa s]
µm: viscosidad en el punto de fusión [Pa s]
σL, σS: conductividades del líquido y sólido en el punto de fusión [Ω-1·m-1]
νs: frecuencia vibracional del sólido calculada a partir de la temperatura Dabye θD
(presentado en la tabla 2-8 de la pág. 37 del Hines y Maddox) [s-1] (νs = θD k / h)
 Cuando cambia la temperatura, el coeficiente de difusión cambia de la forma:

 DAB
0
µ  DAB
0
µ
  =   (1.40)
 T T1  T T2

En rangos de temperatura moderados, los coeficientes de difusión se relacionan

con la temperatura a través de la expresión de Arrhenius:

 −Q 
DAB = D0 exp   (1.41)
 RT 

donde:
D0: factor de frecuencia [m2/s]
Q: energía de activación [energía/mol]
Los datos para algunos elementos metálicos se presentan en la tabla 2-9 (pág. 39
del Hines y Maddox).
• Coeficiente de Difusión para Sistemas Multicomponentes
- Para un sistema gaseoso:
n
N i − xi ∑ N j
j =1
Dim = − (1.42)
n xi N j − x j N i
∑
j =1 Dij

donde para una mezcla estática Nj = 0. Entonces, la expresión para una mezcla
estática de gases A, B, C y D, donde A es el que difunde (Wilke, 1950):

1 − xA
DAm = (1.43)
xB DAB + xC DAC + xD DAD
 Para un sistema líquido: Leffler y Cullinan (1970) propusieron la siguiente ecuación
para un sistema líquido ternario, donde A es el soluto diluido que difunde a través de
la mezcla concentrada de B y C:

µ ABC = ( DAB µB ) (D µC )
0 0 xB 0 xC
DAm AC (1.44)

• Coeficiente de Difusión en Sólidos

- Sherby y Simnad (1961) desarrollaron una ecuación empírica para predecir los
coeficientes de autodifusión en sólidos:

 T 
DAA = D0 exp  − ( K 0 + V ) m  (1.45)
 T 
donde:
V: valencia
Tm: temperatura del punto de fusión [K]
K0 = 14 para metales bcc (γ-U, Na, α-Fe, Li, Nb, β-Ti, Cr, Ta, β-Zr), 17 metales fcc
(Pt, Au, Ag, Cu, β-Co, γ-Fe, Ni, Pb) y hcp (Mg, Zn, Cd, α-Ti, α-Zr, complejo hex. β-
U y grafito) y 21 para la estructura de diamante (Ge, β-Sn).
En la ecuación (1.43), el factor exponencial corresponde directamente a la energía de
activación Q entre la temperatura T. Para diversos sólidos, Askill (1970) proporciona
una tabla con datos de la energía de activación, el factor pre-exponencial y la
temperatura Tm.
- Para sólidos porosos:
- Difusión ordinaria: el diámetro de poro del material es grande en comparación
con la trayectoria libre media de las moléculas del gas.
- Difusión tipo Knudsen: el diámetro del poro del material es menor en
comparación con la trayectoria libre media de las moléculas del gas.
- Difusión superficial: las moléculas se absorben sobre la superficie del material
y después se transportan hacia los lugares de menor concentración.
Para la difusión de Knudsen se tiene:
 Ud
DA, K = (1.46)
3

donde:
U: velocidad de la molécula de gas
d: diámetro del poro
Para poros rectos, redondos, la difusividad se evalúa como:

12
 T 
DA , K = 97r   (1.47)
 MA 
donde:
r: radio del poro [m]
T: temperatura [K]
MA: peso molecular del componente A
El radio del poro se puede evaluar como:
2ε 2V
r= = p (1.48)
SρB S

donde:
S: área superficial del sólido poroso [m2/kg]
ρB: densidad de bulto de la partícula sólida [kg/m3]
Vp: volumen de poro específico de una partícula sólida [m3kg]
ε: porosidad del sólido
Con el objeto de considerar el camino tortuoso de la molécula y la porosidad del
material, se define una difusividad de Knudsen efectiva:

ε
DA , K , e = DA , K (1.49)
τ

donde:
 τ : factor de tortuosidad
Como también en este tipo de sistemas la difusión de bulto también contribuye al
transporte molecular total a través de los poros, la difusividad efectiva de bulto se
puede escribir como:

ε
DAB ,e = DAB (1.50)
τ

Entonces, para el caso de autodifusión o contra difusión equimolar, la difusividad

efectiva en un material poroso se puede escribir como:

1 1 1
= + (1.51)
DA,e DA, K ,e DAB ,e

Algunos valores de porosidad y tortuosidad son dados en la tabla 2-11 (pág. 47

del Hines y Maddox).

I.4. Relaciones de Maxwell-Stefan

Si se considera un sistema binario gaseoso, del balance de mometum que describen
las colisiones entre las moléculas de la especie 1 y moléculas de la especie 2, se puede
obtener (Taylor y Krishna, 1993):

∇p1 = − f12 x1 x2 ( u1 − u2 ) (1.52)

donde f12 es un parámetro empírico análogo al factor de fricción o coeficiente de arrastre.

Por conveniencia, se puede definir un coeficiente inverso de arrastre:

D12 = p f12 (1.53)

Sustituyendo la Ec. (1.53) en (1.52):

 ∇p1 x x (u − u )
d1 = =− 1 2 1 2 (1.54)
p D12

siendo:
d1: fuerza impulsora de la difusión de la especie 1 en una mezcla binaria de gas
ideal
D12 : difusividad de Stefan-Maxwell
De forma similar, puede plantearse la ecuación para la especie 2:

∇p2 x x (u − u )
d2 = =− 1 2 2 1 (1.55)
p D21

Si se considera la presión constante a lo largo de la trayectoria de difusión, las

expresiones se pueden escribir como:

x1 x2 ( u1 − u2 )
∇x1 = − (1.56)
D12

x1 x2 ( u2 − u1 )
∇x2 = − (1.57)
D21

Para sistemas multicomponentes, la ecuación (1.54) puede ser reescrita como:

n xi x j ( ui − u j )
d i = −∑ , i = 1, 2… , n − 1 (1.58)
j =1 Dij

La expresión (1.58) puede escribirse en término de los flujos molares Ni:

n xi N j − x j N i
di = ∑ , i = 1, 2 …, n − 1 (1.59)
j =1 cDij
o también en término de los flujos de difusión Ji:

n xi J j − x j J i
di = ∑ , i = 1, 2 …, n − 1 (1.60)
j =1 cDij

Es importante mencionar que sólo n-1 ecuaciones son independientes porque la suma
de los di debe suma cero. Por ello se puede demostrar que para una mezcla de gas ideal
(Taylor y Krishna, 1993):
Dij = D ji (1.61)

Por otra parte, también se puede demostrar que para un sistema binario de gases
ideales se cumple que:

J i = −cD12 d1 = cD12∇x1 (1.62)

En el caso de fluidos no ideales, la fuerza impulsora de difusión debe definirse en

función de los gradientes de potencial químico:

xi
di = ∇ T , p µi (1.63)
RT

Los subíndices T y p se usan para enfatizar que el gradiente debe calcularse a

condiciones de temperatura y presión constante. La Ec. (1.61) puede escribirse en
término de los coeficientes de actividad, γ i :

n
d i = ∑ Γij ∇x j = Γ∇xi (1.64)
j =1

donde la matriz del factor termodinámico, Γ viene dada por la expresión:

 ∂ ln γ i
Γij = δij + xi (1.65)
∂x j
T , p, ∑

donde δij es el delta de Kronecker que se define como:

δij = 1 , para i = j
δij = 0 , para i ≠ j

El signo de sumatoria en la ecuación (1.65) se usa para indicar que la diferenciación

con respecto a la fracción molar debe ejecutarse de forma que la suma de las fracciones
molares siempre sea igual a la unidad.
Si se manipulan las ecuaciones (1.59), (1.60) y (1.64) se obtiene:

n xi J j − x j J i n xi N j − x j N i
Γ∇xi = ∑ =∑ , i = 1, 2, … n − 1 (1.66)
j =1 cDij j =1 cDij

Para una mezcla de gases ideales, se tiene:

n xi J j − x j J i n xi N j − x j N i
∇xi = ∑ =∑ , i = 1, 2, … n − 1 (1.67)
j =1 cDij j =1 cDij

Se puede encontrar una relación entre la difusividad de Maxwell-Stefan y el

coeficiente de difusividad de Fick. La relación en un sistema binario se puede expresar
como:

D12 = ΓD12 (1.68)