Qué es el

transporte de masa ?

El transporte de masa es el fenómeno de transporte que se

dedica a estudiar el movimiento de las especies químicas en
función de:

• La difusividad de las especies (molecular)

• La velocidad del fluido (convectivo)
• Las reacciones químicas (otros, fuente/sumidero)
 Fases: Sólido, líquido, gas
• MEDIO Propiedad de transporte: Difusividad Ɗ [m2/s],[cm2/s]

• INTERFASE Área interfacial (Perpendicular al flujo)

Fuerza motriz: Diferencia de concentración

• DIFERENCIA DE Fracción másica de A (wA)
CONCENTRACION Fracción molar de A (xA)

Flujo Densidad de flujo

wA nA
Másico [kg/s] [kg/m2*s]
FLUJO [lb/s] [lb/ft2*s]
Flujo de un componente A
WA NA
Molar [mol/s] [mol/m2*s]
Ecuación de balance general
 Ecuación de balance general
Balances microscópicos
Sistema isotérmico, sin esfuerzos viscosos

Acumula = Entra – Sale  Genera

(molecular, convectivo) (otros)

Balance de masa global - Ecuación Continuidad (,v)

Masa
Balance por componente

Balance de energía mecánica, eléctrica, etc.

Energía
Balance de energía térmica (radiación)

Cantidad
Movimiento Balance de movimiento - fuerzas
 Ecuación de balance general
Balances microscópicos
Sistema isotérmico, sin esfuerzos cortantes

Acumula = Entra – Sale  Genera

(molecular, convectivo) (otros)

Masa Balance por componente

(Composición)

Versión másica
(Recomendado sin Rxn química)
wA = fracción másica de A [-]

Versión molar
(Con y sin Rxn química)

xA = fracción molar de A [-]

Ecuación de balance de masa por componente [1]
Acumula = Entra – Sale  Genera

Versión másica

Convectivo Molecular Otros

wA = fracción másica de A

[1] Bayron BIRD, Warren STEWART y Edwin LIGHTFOOT. Transport Phenomena. 2da Edición. Apendice B
Ecuación de balance de masa por componente [1]

Acumula = Entra – Sale  Genera

Versión molar

Convectivo Molecular Otros

c = densidad molar de A [mol/m3]
xA = fracción molar de A [-]
cxA = CA = concentración molar de A [mol/m3]
Cuando hay reacción química las moles no se conservan, por tanto este balance es
válido solo para mezclas binarias o mezclas de peso molecular similar

[1] Bayron BIRD, Warren STEWART y Edwin LIGHTFOOT. Transport Phenomena. 2da Edición. Apendice B
Ecuación de balance de masa por componente [1]

Acumula = Entra – Sale  Genera

Para mezclas binarias o mezclas de

peso molecular similar
v  v*
[1] Bayron BIRD, Warren STEWART y Edwin LIGHTFOOT. Transport Phenomena. 2da Edición. Apendice B
Condiciones de frontera típicas
Para transporte de Masa

CA, xA = constante
CA, xA = f (r,z,…)
Concentración
El compuesto A entra o sale a
una concentración conocida
CA = CA0

n.NA = constante NA
CA = CA0
Flujo n.DABCA = f (CA,xA)
n.NA = kte [mol/m2*s]
El compuesto A entra o sale

n.NA = 0
a un valor de flujo conocido

n.NA = 0
CA = 0

n.NA = 0
No flujo dCA/dz = 0
n.NA = 0 [mol/m2*s]
(Aislado) El compuesto A no fluye
a través de la frontera
Condiciones de frontera típicas
Para transporte de Masa

Coeficiente global Concentración

de transf. de masa de referencia

Flujo n.NA = kC*(CA – CA*)

convectivo CA = CA0
El compuesto A entra o sale
a un valor de flujo conocido
NA

n.NA = kC(CA-CA*)

NA.n = 0 CA = C A *
Simetría Existe un comportamiento
(Continuidad) espejo a lo largo de la
frontera de simetría

n.NA = 0
Condiciones de frontera típicas
Para transporte de Masa

Condición inicial
Es una condición aplicable solo cuando el problema opera en estado transitorio

La condición inicial define sobre toda la geometría

el valor de cada variable dependiente en el
instante que comienza el fenómeno (t=0)

CA (t=0) = CA0

CA, xA = 0
En t = 0 CA, xA = f (r,z,…)
etc…

CA (t=0) = 0
Fenómenos de transporte ofrece una
descripción matemática exacta de El problema es que no existe la
prácticamente todos los fenómenos matemática suficiente para
naturales e industriales resolver todos los problemas

Ecuación exacta
del fenómeno

Solución exacta Solución aproximada Ecuación aproximada a la

(métodos analíticos) (métodos numéricos) ecuación exacta
(Muy pocos casos) (COMSOL) (Números adimensionales)

Solución exacta o numérica

a la ecuación aproximada
EJEMPLO 1: BSL Ejercicio 18.5
Considere una columna de paredes húmedas utilizada industrialmente para la
absorción de un soluto A gaseoso dentro de una capa líquida B. Suponga que A es
ligeramente soluble en el líquido B y que el flujo entre el gas y el líquido es en
paralelo. Describa el perfil de concentración de A dentro del líquido B

Flujo en
paralelo
Gas Líquido

Flujo en
contracorriente
Líquido

Gas
EJEMPLO 1:
Describa el perfil de concentración de A dentro del líquido B

. . . .. . . DATOS:
. . .. . … Espesor de película  = 2 mm
.. . . .. . . Longitud película L = 100 cm
. . . .. . Velocidad máxima Vmax = 0.005 cm/s
. . . C.A0. . . Difusividad A en B DAB = 1x10-6 cm2/s
Líquido (B)

. Mezcla
. . . Gaseosa
.. . . Concentración inicial CA0 = 0.5 mol/cm3
Pared . . . (A)
.. . .
. . . . .. . .
. . . .. .
. . . .. . .
. . . . .. . .
. . . .. .
. . . .. . .


En la capa de líquido es
donde ocurre la
transferencia de masa
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 1

1) Definir la geometría del medio y las fronteras

x
z Medio: Líquido B

Fronteras: Interface Gas/Líquido (x=0)

Interface Líquido/Pared (x=)
Entrada de líquido (z=0)
Salida de Líquido (z=L)

Coordenadas: Cartesianas, Cilíndricas, Esféricas

(x, y, z)
Dimensiones: 1D, 2D, 3D

Ubicación del eje de coordenadas

Interface Gas/Líquido
y Entrada de líquido
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 1

2) Plantear las suposiciones del problema

Ecuación de balance general

Balances microscópicos
Sistema isotérmico, sin esfuerzos viscosos

Acumula = Entra – Sale  Genera

(molecular, convectivo) (otros)

Balance de masa global - Ecuación Continuidad (,v)

Masa
Balance por componente

Balance de energía mecánica, eléctrica, etc.

Energía
Balance de energía térmica (radiación)

Cantidad
Movimiento Balance de movimiento - fuerzas
SOLUCION EJEMPLO 1

2) Plantear las suposiciones del problema

Versión másica

Versión molar

x
Tipo de flujo: Laminar vs Turbulento
z Vz=Vmax(1-(x/)2)
x0 VzVmax

Acumulación: Estable vs Transitorio

kte Vmax
Convectivo: c, v = (vx, vy, vz)
kte Convección z >> Difusión z
Molecular: DAB, XA = (XAx, XAy, XAz)

Otros: rA
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 1

3) Depurar las ecuaciones de balance

Ecuación de balance de masa por componente

 2C A C A
cxA = CA D AB  vmax 0
vz  vmax x 2
z
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 1

4) Plantear las condiciones de frontera

Opciones de condiciones de frontera:

x
CA, xA = conocida Concentración conocida
z n.NA = f (CA,xA) Flujo
n.NA = kC*(CA – CA*) Flujo convectivo
n.NA = 0, dCA/dz = 0 No flujo
Simetría

Condiciones frontera:
En Z = 0: CA = 0
En X = 0: CA = CA0
En X = : dCA/dx = 0 A no se difunde en
A es ligeramente
soluble en B En X = : CA = 0 la pared
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 1
5) Obtener perfiles: Resolver las ecuaciones de balance + condiciones de frontera

 2C A C A
D AB  v 0
z
z , max
x 2

Difusión en x Convección en z

Cómo se resuelve?
Método de sustitución
 2C A C A
DAB  C A  f ( x, z )
x 2
z vmax 

x
Sustituyo las dos variables
(x,z) por una sola variable  :

(eta) 4 D AB ( z / vmax )

d 2C A dC A
 2  0 C A  f ( )
d 2
d
 d 2C A dC A x
 2  0 
d 2
d 4 D AB ( z / vmax )

Condiciones límite: Nuevas condiciones límite:

En Z = 0: = CA = 0 = CA = 0
En X = 0: =0 CA = CA0 =0 CA = CA0
En X = : = CA = 0

Resolviendo la ecuación diferencial…

Por tablas del BSL se sabe que:

Reemplazando:



C A  C1 e  d  C 2
2 2
d CA dC A
 2  0
d 2
d 0

Esta es la función error

 


2
e  d
2
Función error erf ( ) 

0

Evaluando las condiciones límite:

C A  C1erf ( )  C2
0
=0 CA = CA0 C A0  C1erf (0)  C2 C2=CA0

1
= CA = 0
0  C1erf ()  C2 C1=-C2
 C A  C1erf ( )  C2

C A  C A0  C A0 erf ( )
CA
 1  erf ( )
C A0

CA  x 
 1  erf  
C A0  4D ( z / v ) 
 AB max 

Perfil de composición de A en la dirección x y z

SOLUCION EJEMPLO 1
5) Obtener perfiles: Resolver las ecuaciones de balance + condiciones de frontera

Absorción columna pared mojada

Perfil 2D de concentración de un gas A difundiendose en un líquido B

CA  x 
 1  erf  
C A0  4 D AB ( z / v max ) 


Espesor de película  0.20 [cm]

Difusividad A en B D AB 1.00E-06 [cm2/s]
Longitud pelicula L 100 [cm]
Velocidad máxima V max 0.005 [cm/s]
Concentración inicial de A C A0 0.5 [mol/cm3]
Posición transversal en la película X 0.05 [cm]
Posición longitudinal en la película Z 100 [cm]
Concentración de A en (X,Z) CA 0.40 [mol/cm3]
Grafique el perfil 2D
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 1
6) Estimar flujos y concentraciones deseadas

❶ Deducciones lógicas:
Si el espesor de película aumenta, CA = ?
Si la longitud de película aumenta, CA = ?
Si la difusividad de A en B aumenta, CA = ?

❷ Cuán larga debe ser una torre de pared mojadas para que el compuesto A penetre
mínimo el 75% de la película de líquido con un espesor  =2 mm?

❸ Cuánto debería ser espesor de película ( ) para que la concentración promedio de A

dentro de la película sea igual a CA0/2 en Z = L?
EJEMPLO 2: BSL Ejercicio 18.4
En un recipiente de laboratorio de base cuadrada se tiene un líquido B sobre el cual
reposa un gas A. Se ha detectado que una pequeña cantidad de gas se difunde hacia
el interior dentro del líquido y reacciona para formar el compuesto AB. Estime el
perfil de concentración dentro del líquido.

A se difunde en B y
reacciona dentro de B DATOS :
A + B  AB DAB = 1x10-5 cm2/s
CA0 = 0.5 mol/cm3
L = 20 cm
Reacción de primer orden rA = -kCA
k = 5x10-6 s-1
CA0
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 2

1) Definir la geometría del medio y las fronteras

Medio: Líquido B

Fronteras: Interface Gas/Líquido (z=0)

Fondo del tanque (z=L)
CA0
x
Z=0
z Coordenadas: Cartesianas, Cilíndricas, Esféricas
(x, y, z)
L
Dimensiones: 1D, 2D, 3D

Z=L Ubicación del eje de coordenadas

Interface Gas/Líquido
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 2

2) Plantear las suposiciones del problema

Ecuación de balance general

Balances microscópicos
Sistema isotérmico, sin esfuerzos viscosos

Acumula = Entra – Sale  Genera

(molecular, convectivo) (otros)

Balance de masa global - Ecuación Continuidad (,v)

Masa
Balance por componente

Balance de energía mecánica, eléctrica, etc.

Energía
Balance de energía térmica (radiación)

Cantidad
Movimiento Balance de movimiento - fuerzas
SOLUCION EJEMPLO 2

2) Plantear las suposiciones del problema

Versión másica

Versión molar

Tipo de flujo: Laminar vs Turbulento

Acumulación: Estable vs Transitorio
kte
Convectivo: c, v = (vx, vy, vz)
kte
Molecular: DAB, XA = (XAx, XAy, XAz)
-k*CA Reacción de
Otros: rA primer orden
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 2

3) Depurar las ecuaciones de balance

Ecuación de balance de masa por componente

 2C A
cxA = CA
DAB  kC A  0
z
rA = -k*CA 2

Difusión en z Rxn Química

 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 2

4) Plantear las condiciones de frontera

Opciones de condiciones de frontera:

CA, xA = conocida
n.NA = f (CA,xA)
n.NA = kC*(CA – CA*)
n.NA = 0, dCA/dz = 0

Condiciones frontera:
En Z = 0: CA = CA0
En Z = L: dCA/dz = 0 A no se difunde en
la pared
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 2
5) Obtener perfiles: Resolver las ecuaciones de balance + condiciones de frontera
 2C A
DAB  kC A  0
z 2

Difusión en z Reacción química

Cómo se resuelve?
Por tablas del BSL:

Por lo tanto:

C A  C1 cosha  z   C2 sinha  z 

a 
2 k  k   k 
DAB C A  C1 cosh z   C2 sinh  z 
 DAB   DAB 
 C A  C1 cosha  z   C2 sinha  z 

Evaluando condiciones límites:

En Z = 0: CA = CA0
1 0
C A0  C1 cosha(0)  C2 sinha(0)  C1

En Z = L: dCA/dz = 0

 C A0 cosha  z   C2 sinha  z   ?
dCA d
dz dz
Por tablas del BSL

 C A0  a sinha  z   C2  a cosha  z 
dC A
dz
En Z = L: dCA/dz = 0

 C A0  a sinha  L   C2  a cosha  L  0
dC A
dz

C A0 sinha  L 
C2  
cosha  L 

Reemplazando las constantes halladas

C A  C1 cosha  z   C2 sinha  z 

sinh(a  L)
C A  C A0 cosha  z   C A0 sinha  z 
cosh(a  L)
Redefiniendo
(phi)   aL C A cosh 1   
(zeta)  
z 
C A0 cosh( )
L
 Perfil de concentración en un tanque con rxn química
C A cosh  1   

C A0 cosh( )
Profundidad del tanque L 20 [cm]
Constante de reacción k 1.00E-05 [s-1]
2
Difusividad de A en B DAB 5.00E-05 [cm /s]
3
Concentración de A en la interfase CA0 0.5 [mol/cm ]
Constante Fi (a*L)  8.94 [-]
constante a (k/DAB)^0.5 a 0.45 [cm-1]
Variable posición (Z) Variable posición (Z/L) Concentración de A
Z [cm]  [] CA [mol/cm3]
0 0 0.5000000
2 0.1 0.2044209
4 0.2 0.0835758
6 0.3 0.0341694
8 0.4 0.0139701
10 0.5 0.0057122
12 0.6 0.0023369
14 0.7 0.0009591
16 0.8 0.0004012
18 0.9 0.0001862
20 1 0.0001305

Perfil de concentración de A
0.6
Concentración de A [mol/cm3]

0.5

0.4

0.3

0.2

0.1

0.0
0 5 10 15 20 25
Profundidad del tanque [cm]
 Cómo se resuelve un problema típico de fenómenos?

1. Definir la geometría del medio + frontera

2. Plantear las suposiciones del problema (Fenómenos, términos de la ecuación balance)
3. Depurar las ecuaciones de balance a resolver (aplicar suposiciones)
4. Plantear las condiciones de frontera sobre la geometría
5. Obtener los perfiles resolviendo las ecuaciones de balance + condiciones frontera
(método analítico/numérico)
6. Estimar flujos y concentraciones deseadas (max/min/promedios, puntuales, planos,
etc.)
SOLUCION EJEMPLO 2
6) Estimar flujos y concentraciones deseadas

❶ Cuánto es la concentración promedio de A dentro del tanque?

❷ Cuánto es el flujo molar de A en:

La superficie del líquido (Z=0) ?
La mitad del tanque (Z=L/2) ?
El fondo del tanque (Z=L) ?

❸ En dónde se ubica el punto de concentración más alta? más baja?

❹ Cuánta es la conversión promedio de A dentro del tanque ?

Conversión = (CA0 - CA ) / CA0
❶ Cuánto es la concentración promedio de A dentro del tanque?
L 1
Por definición un promedio en
coordenadas cartesianas 1D es: CA 
0
C dz A

 C d
0
A
L 1
 dz 
1
d 
0 0

1C  1    d
A0 cosh 
  cosh 1   d
C A0 1
CA  
0 cosh( ) cosh( ) 0

Por tablas del BSL

tanh 
sinh 
C A0 1
CA   C A0   aL
cosh( )  
❷ Cuánto es el flujo molar de A en: La superficie del líquido (Z=0) ?

FLUJO = Molecular + Convectivo

NA  cDAB X A  C Av

En z = 0 : N A z 0  cDAB X A z 0   DAB C A
z z 0

Si Entonces:
C A cosh 1    11  dd  ( z / L) (  / L)

dC A
C
?CAA000
C
sinh
cosh
cosh 
1(z /L ) 
C A0 cosh( ) dz cosh()) dz
cosh( dzcosh( )

C A
En z = 0 : NA z 0
  D AB
z z 0

DAB C A0
N A z 0   tanh( )
L
 BSL Ejercicio 19.1.1
TAREA :
Un gas A se mueve a través de un reactor tubular en cuyo centro existe un sólido poroso B
(catalizador). El gas se difunde dentro del sólido y a medida que se difunde produce una
reacción química.
Encuentre el perfil de concentración de A dentro del sólido poroso asumiendo que el perfil
de la velocidad de flujo es plano (completamente chato). Desprecie los efectos del calor de
la reacción.

Catalizador Sólido
poroso

Conducción en z

Convección en z Reacción química

Fenómenos de transporte ofrece una
descripción matemática exacta de El problema es que no existe la
prácticamente todos los fenómenos matemática suficiente para
naturales e industriales resolver todos los problemas

Ecuación exacta
del fenómeno

Solución exacta Solución aproximada Ecuación aproximada a la

(métodos analíticos) (métodos numéricos) ecuación exacta
(Muy pocos casos) (COMSOL) (Números adimensionales)

Solución exacta o numérica

a la ecuación aproximada
 EJEMPLO 1 COMSOL:
Halle el perfil de concentración, la concentración promedio de A en el líquido, el flujo molar
de A en la superficie del líquido y la conversión de la reacción.
Asuma las siguientes condiciones:
SOLUCION:
k = 5x10-6 s-1
DAB = 1x10-5 cm2/s Suposiciones:
CA0 = 0.5 mol/cm3 Coordenadas rectangulares 1D
L = 20 cm Estado estable, isotérmico, sin convección
rA = -kCA (reacción de 1er orden)
Conversión = (CA0 - <CA> ) / CA0 Condiciones de frontera:

x
z

Z=0 CA = CA0

A + B  AB
L

Z=L dCA/dz = 0
(Aislado – No flujo)
 EJEMPLO 2 COMSOL:
Cuánto es la concentración promedio de A después de 1, 2 y 3 días de estar almacenado
el gas A dentro del tanque con el líquido B?
SOLUCION:
Resolver en estado transitorio
Planteamiento en Estado estable
Estado transitorio = +
Punto o condición inicial
CA (t=0) = 0 [mol/cm3]

x
z

Z=0 CA = CA0

A + B  AB
L

Z=L dCA/dz = 0
(Aislado – No flujo)
 EJEMPLO 3 COMSOL:
Resuelva el mismo problema del Ejemplo 1 pero en coordenadas 2D

SOLUCION:
Suposiciones:
Coordenadas rectangulares 2D
Estado estable, isotérmico, sin convección
Difusión solo en la dirección z

Condiciones de frontera:
x
z

Z=0 CA = CA0

A + B  AB
L

Z=L dCA/dz = 0
(Aislado – No flujo)
 EJEMPLO 4 COMSOL:
Considere una columna de paredes húmedas utilizada industrialmente para la absorción
de un soluto A gaseoso dentro de una capa líquida B. Describa el perfil de concentración
de A dentro del líquido B.

Condiciones de frontera: Asuma las siguientes condiciones:

CA = 0 . Espesor de película  = 2 mm
. Difusividad A en B DAB = 1x10-6 cm2/s
. Longitud pelicula L = 100 cm
. Velocidad máxima Vmax = 0.001 cm/s
. Concentración inicial CA0 = 0.5 mol/cm3
Líquido (B)
dCA/dx = 0

CA = CA0

A
L SOLUCION:
Suposiciones:
Estado estable e isotérmico
Coordenadas rectangulares 2D
No hay reacción química
Vz  Vmax

?
 EJEMPLO 5 COMSOL:
Si se arranca a operar el proceso de absorción (concentración de A en al película = 0)
entonces cuánto tiempo demoraría la columna en llegar al estado estable?

Condiciones de frontera:
CA = 0
Líquido (B)
dCA/dx = 0

CA = CA0

A
L

  x 

C A  C A0 1  erf  

  4 DAB ( L / vmax ) 
Fenómenos de transporte ofrece una
descripción matemática exacta de El problema es que no existe la
prácticamente todos los fenómenos matemática suficiente para
naturales e industriales resolver todos los problemas

Ecuación exacta
del fenómeno

Solución exacta Solución aproximada Ecuación aproximada a la

(métodos analíticos) (métodos numéricos) ecuación exacta
(Muy pocos casos) (COMSOL) (Números adimensionales)

Solución exacta o numérica

a la ecuación aproximada
Análisis dimensional
Técnica que busca reescribir las ecuaciones de balance general, condiciones
límite e iniciales con el objeto de proveer números adimensionales que
caractericen los fenómenos
Ecuación Navier-Stokes Número de Reynolds

D Dv  Convectivo
 v     2v Re 
Dt  Molecular

Para el balance de masa por componente se puede hacer lo mismo

D
Balance masa por componente C A  DAB 2C A
Dt
Si reagrupo DAB y redefino D 1
CA   CA
2
la variable CA se puede Dt PeM
obtener:
Número de Peclet (Masa)
t , v, C A Variables adimensionales lv0 Transp. Convectivo
PeM 
PeM Número adimensional D AB Transp. Difusivo
Numeros adimensionales en el transporte de masa:

lv0 Transp. Convectivo

Peclet PeM 
D AB Transp. Difusivo
(Masa)

PeM  Re* Sc
 Transp. Viscoso
Schimdt Sc 
D AB Transp. Difusivo

k l kl Transp. masa en la interfase

Sherwood Sh  x  c
cDAB D AB Transp. masa de A en el medio B

kx, kc son coeficientes de transferencia de masa

Grasoff, etc.
 Balance de cantidad Balance de masa
de movimiento
Velocidad
Perfil Concentración y flujo molar
Presión
(2D, 1D, 0D) Concentración y flujo másico
Esfuerzo cortante

Promedio Velocidad
Máximo Concentración y flujo molar
Presión
Mínimo Concentración y flujo másico
Esfuerzo cortante

Números Peclet de masa

Adimensionales Reynolds Schimdt,
Sherwood, etc.

Aplicación
Números Factor fricción = f(Re) Coef. Transf. Masa = f(Pe, Sc, Sh, Re, ...)
Adimensionales
Coeficientes de transferencia de masa

FLUJO = Molecular + Convectivo

NA  cDABX A  C Av
(Analogía a la Ley de Fick)
Flujo = Propiedad*Gradiente
NA  cDABX A  C Av Para sistemas
complejos N A  k x C A
Se requiere conocer
mol m mol
el perfil de concentración
m2*s s m3

kx = f (Pe, Sc, Re, etc..)

(DAB, v, l, D, , etc…)
Los coeficientes de transf. de masa pueden ser estimados teórica y empíricamente
(lo más común)
Para qué sirven estos coeficientes?
1. Para hallar el flujo o la densidad de flujo
2. Para escalar procesos donde intervengan fenómenos de transporte
3. Para diseñar equipos donde intervengan fenómenos de transporte
Transferencia de A a través de B (sin difusión de B en A)
Unidades de los coeficientes k?

N A  k g PA moles transferidas

tiempo*área*presión
N A  k y YA moles transferidas
tiempo*área*fracción_mol
Gases
N A  kc C A moles transferidas
tiempo*área*(moles/vol)
n A  k y YA moles transferidas
tiempo*áarea*(masa A/ masa B)

N A  k L C A moles transferidas
tiempo*área*(moles/vol)
Líquidos moles transferidas
N A  k x X A tiempo*área*fracción_mol
Usos de los coeficientes de transferencia de masa
Ejemplo 1:

x
Flujo sobre una placa plana Fluido (B)
z Placa
L

ReL < 105 Sh  0.66 Re L1/ 2 Sc1/ 3

ReL > 105 Sh  0.36 Re L 0.8 Sc1/ 3

Lv   kc L
donde: Re L  Sc  Sh 
 D AB D AB
Usos de los coeficientes de transferencia de masa
Ejemplo 2:

Flujo sobre una esfera

Sh  2.0  0.4(Re D1/ 2  0.06 Re D 2 / 3 )Sc0.4

donde:
Dv 
Re D 


La menor
El fluido puede ser gotas de D0 entonces Sh  2 resistencia está en
líquido, burbujas de gas, la interfase
partículas pequeñas de sólido
Usos de los coeficientes de transferencia de masa
Ejemplo 3:

Flujo en una tubería

ReD < 2100

Flujo laminar Sh  8 / 3

ReD > 2100

Flujo turbulento
Sh  0.023 Re D 0.8 Sc1/ 3

donde:
Dv 
Re D 

Usos de los coeficientes de transferencia de masa
Ejemplo 4:

Empaques de anillos Rasching y sillas Berl, etc.

Lechos empacados Aplicado a destiladores, absorbedores, etc.

0.36
Si se transfiere algún kGas PB Sc 2/3
 dsG' 
 1.195  
compuesto dentro del gas
G 
 Gas (1   ' ) 
0.45
Si se transfiere algún k Liq d S  dsL' 
 25 .1 
0.5
ScLiq
  Liq 
compuesto dentro del líquido
DALiq
Shullman y colaboradores (1955)
kgas Coef. transferencia masa en el gas
Kliq Coef. transferencia masa en el líquido
ds diámetro equivalente del empaque
 fracción de vacío del lecho
L’ vel. superficial másica del líquido
G’ vel. superficial másica del gas
G vel. superficial molar del gas
Ejemplo 4:
Se tiene un lecho empacado a escala laboratorio compuesto de Sillas Berl de 38 mm (1.5 in)
que es utilizado para la absorción de vapor de benceno (A) que viene en un gas de refinería
a una concentración de 0.25 mol/cm3. El lecho opera haciendo pasar un disolvente acuoso
que posee la propiedad de solubilizar selectivamente el vapor de benceno.
Hallar el flujo global molar que se transfiere dentro del líquido de la torre.

1er paso Hallar coeficiente de transferencia de masa de líquido (kLiq.)

0.45
k Liq d S  dsL' 
 25 .1 
0.5
ScLiq
DALiq   Liq 

2do paso Calcular flujo molar con N A  k Liq C A

Datos:

Líquido entra

Flujo = 0.31 [kg/s]

 = 834 [kg/m3]
 = 2 cpoise
DALiq = 4.71x10-6 [cm2/s]
CA0 = 0
 = 0.717
d=38 cm

Gas entra

Flujo = 81.6 [g/s]

Peso mol = 11 [g/mol]
 = 0.472 [kg/m3] N A  k Liq C A  ?
 = 0.01 cpoise
DAGas = 0.13 [cm2/s] C A  (C A0  0)
P = 2 [atm]
CA0 = 0.25 mol/cm3
1er paso Hallar coeficiente de transferencia de masa de líquido (kLiq.)

[mol/m2*s]*[m3/mol]
+ Robert TREYBAL. Operaciones de transferencia de masa. Segunda edición. New York: McGraw Hill (1991)
2do paso Calcular flujo molar con N A  k Liq C A

N A  k LiqC A0

NA  1.16x10-4 [mol/m2*s][m3/mol] * 0.25 [mol/cm3]

NA  29 [mol/m2*s] de Benceno

2,262 [g/m2*s] de Benceno

Si el área interfacial de líquido (área de transferencia de masa) es 100 m2

Entonces:

Flujo molar de A = 2,900 [mol/s] de Benceno

Flujo másico de A = 226,2 [g/s] de Benceno
Con cuál tipo de empaque obtendría mayor flujo de benceno transferido?

Diámetro de la Coeficiente de
Tamaño nominal
esfera equivalente transferencia de masa
ds k Liq
[mm] [in] [m] [m/s] x 104
13 0.50 0.0177 1.98
Anillo Rasching 25 1.00 0.0356 1.35
cerámico 38 1.50 0.0530 1.09
50 2.00 0.0725 0.91
25 1.00 0.0130 2.35
Anillo Rasching
38 1.50 0.0543 1.07
carbón
50 2.00 0.0716 0.92
13 0.50 0.3162 0.41
Silla Berl
25 1.00 0.0320 1.43
cerámica
38 1.50 0.0472 1.16

kLiq mas alto representa el empaque que favorece mayor transferencia de masa
Por tanto, el Anillo Rasching de carbón
de 1 in de tamaño es el mejor
 Coeficiente de transferencia de masa en un lecho empacado
0.36
k Gas PB Sc 2/3
 dsG ' 
 1.195  
G 
 Gas (1   ' ) 
Presión del gas PB 2.00E+05 [Pa]
Diámetro de la esfera equivalente ds 0.0472 [m]
Fracción de vacío  0.72 [-]

Flujo másico de gas n Gas 0.08 [kg/s]

Peso molecular del gas M Gas 1.10E-02 [kg/mol]
Diámetro del lecho d 0.38 [m]
Area transversal del lecho a 0.11 [m2]
Vel. superficial molar del gas (nGas/a*MGas) G 65.07 [mol/m2*s]
Vel. superficial másica del gas (nGas/a) G' 0.72 [kg/m2*s]
Densidad gas  0.472 [kg/m3]
Viscosidad del gas  Liq 1.00E-05 [Pa*s]
Difusividad de A en el gas D ALiq 1.30E-05 [m2/s]
Número de Schimdt Sc 1.63E+00 [-]

Coeficiente de transferencia de masa k Gas 9.56E-06 [mol/m2*s]*[m3/mol]

 Aplicaciones de la
Transferencia de masa, calor y
cantidad de movimiento
 Separar un compuesto de una mezcla utilizando
Evaporación la diferencia de volatilidad (Temp. ebullición)

Hay transferencia de masa , calor y cantidad de movimiento

Ejemplo:
A un evaporador se alimentan una solución que contiene un % conocido de azúcar
(38%w). Del evaporador se espera retirar parte del agua de la solución y obtener
una solución concentrada en azúcar (74%w).

Dónde ocurre la transferencia?

Cuanta energía se requiere para este equipo ?
Aplicaciones de la transferencia de masa

Separar un o varios compuestos de una mezcla

Destilación utilizando la diferencia de volatilidad en
operación en contracorriente
Ejemplo:

Vapor
Concentrado
en el liviano
Líquido

Resistencia Agua Fria

(Calentamiento) (Enfriamiento)

Solución agua-etanol Destilado

(20% etanol) Solución agua etanol
(80% etanol)
Aplicaciones de la transferencia de masa

Ejemplo:
Producción de licores Producción de gasolinas
Aplicaciones de la transferencia de masa
Hay transferencia de masa , calor y cantidad de movimiento
Destilación Dónde ocurre la transferencia?

Torre de platos
Torre de empacada
Aplicaciones de la transferencia de masa

Destilación

Aquí ocurre la
transferencia de
Torre de platos masa
Aplicaciones de la transferencia de masa

Separar un compuesto de una mezcla

Absorción gaseosa utilizando el contacto con un
líquido (solubilidad)
Hay transferencia de masa , calor y cantidad de movimiento
Ejemplo:
En un absorbedor se alimentan una mezcla gaseosa que contiene 20%w de acetona y
el resto de un gas portador. La mezcla gaseosa se trata con una corriente de agua pura.
Se produce un gas de descarga libre de acetona y una solución de acetona en agua.

Dónde ocurre la transferencia?

Cuanta energía se requiere para este equipo ?
Aplicaciones de la transferencia de masa

Adsorción/ Separar un compuesto de una mezcla fluida

utilizando el contacto con un sólido
Intercambio iónico
Aplicaciones de la transferencia de masa

Separar un líquido de una mezcla sólida

Secado utilizando la diferencia de volatilidad
Aplicaciones de la transferencia de masa

Separar un compuesto de una mezcla sólido-líquido

Lixiviación en un disolvente líquido selectivo

Preparación de café Lixiviados industriales (Aceites y Biodiesel)

Rellenos sanitarios
Aplicaciones de la transferencia de masa

Otras aplicaciones industriales:

• Filtración
• Tamizado
• Sedimentación
• Cristalización
• etc…
EJEMPLO 1:
Considere una columna de paredes húmedas utilizada industrialmente para la absorción de
un soluto A gaseoso dentro de una capa líquida B. Describa el perfil de concentración de A
dentro del líquido B.

SOLUCION:
x
.CA0 . . . . . . Suposiciones:
z
. . .. . … Estado estable e isotérmico
.. . . . . . . Coordenadas rectangulares 2D (despreciar efectos en Y)
. . . .. . Líquido fluyen en régimen laminar (Vz=Vmax(1-(x/)2)
. . . .. . . A es ligeramente soluble en B (Vz  Vmax)
Líquido (B)

. Mezcla
. . . Gaseosa
.. . . No hay reacción química (Gas y líquido inertes)
Pared . . . (A) .. . .
Condiciones límite:
. . . . .. . .
En Z = 0: CA = 0
. . . .. .
En X = 0: CA = CA0
. . . .. . .
En X = : dCA/dx = 0 A no se difunde en
. . . . .. . .
En X = : CA = 0 la pared
. . . .. .
. . . .. . . Se puede deducir entonces que:

Vx=Vy=0
En la capa de líquido es Trans.Masa de A en x  0
donde ocurre la
transferencia de masa Trans.Masa de A en z = 0
rA = 0
Definiciones importantes…
Densidad de Flujo total
FLUJO = Molecular + Convectivo

Densidad flujo

Masa
[kg/m2*s]
[lb/ft2*s] nA  jA   Av
Velocidad másica
promedio

Mol [mol/m2*s]
NA  J *
A  C Av *
Velocidad molar
promedio
Densidad de Flujo total
FLUJO = Molecular + Convectivo

Densidad flujo

Masa
[kg/m2*s]
[lb/ft2*s] nA   Dj AAB wA   Av
Ley de Fick (Masa)

Mol [mol/m2*s]
NA  cDAB
*
J
A
x A  C Av *

Ley de Fick (Mol)

Densidad de Flujo total
FLUJO = Molecular + Convectivo

Densidad flujo

Masa
[kg/m2*s]
[lb/ft2*s] nA   Dj AAB wA   Av
m
 Av  wA  n
Ley de Fick (Masa)

 1

Ley de Continuidad

Mol [mol/m2*s]
NA  cDAB
*
J
A
x A  C Av *

m
Ley de Fick (Mol) C Av*  x A  N
 1

Ley de Continuidad
Densidad de Flujo total
FLUJO = Molecular + Convectivo

Densidad flujo

Masa
[kg/m2*s]
[lb/ft2*s] nA   Dj AAB wA  x A N A  ANvB 
Ley de Fick (Masa) Ley de Continuidad
(Mezcla Binaria)

Mol [mol/m2*s]
NA  cDAB
*
AJ x A   C v
wA n A  AnB*

Ley de Fick (Mol) Ley de Continuidad

(Mezcla Binaria)
EJEMPLO 2:
En un recipiente de laboratorio se tiene un líquido B sobre el cual reposa un gas A. Se ha
detectado que una pequeña cantidad de gas se difunde hacia el interior dentro del líquido y
reacciona para formar el compuesto AB. Estime:
El perfil de concentración, la concentración promedio de A en el líquido y el flujo molar de A
en la superficie del líquido
SOLUCION:
A se difunde en B y Suposiciones:
reacciona dentro de B Estado estable e isotérmico
A + B  AB Coordenadas rectangulares 1D (despreciar efectos X e Y)
Reacción de primer orden rA = -kCA No hay convección (Fluido estancado)
Las cantidades AB son pequeñas (Solución seudobinaria)
x
z Condiciones limite:
En Z = 0: CA = CA0
A no se difunde en
En Z = L: dCA/dz = 0 la pared del fondo
CA0
Z=0
Se puede deducir entonces que:
Vz=Vx=Vy=0
L Trans. masa de A en x = 0
Trans. masa de A en z  0
Z=L rA  0
Preguntas típicas de ingeniería:

El gas dentro del tanque está compuesto de dos gases (A y C).

Suponga que el gas A es soluble en el líquido y puede reaccionar con éste mientras que
C ni es soluble ni reactivo. Con cuánta cantidad de gas inerte C debería llenarse el
tanque para disminuir la concentración de A en el fondo del tanque a un 50% menos
que cuando el gas solo está compuesto de A ?

Suponga que A y C son solubles pero solo el gas A reacciona, cuánto sería la cantidad de
gas inerte C para obtener el mismo efecto que en la pregunta anterior?
Coeficientes de transferencia de masa
Interfaciales

NA0

Zona Laminar Zona Turbulenta

v Gradiente de
concentración
varía el 99%
NA>>0

El flujo cerca de las fronteras es aproximadamente

laminar y es donde ocurre la mayoría de la transferencia

Flujo = Propiedad*Gradiente Por análisis dimensional Flujo = Propiedad*Gradiente

C A
n A   DAB Para sistemas n A  k x A C A
z complejos
 sinh( a  L)
C A  C A0 cosh a  z   C A0 sinh a  z 
cosh(a  L)

 cosh a  z   sinh a  z 
CA sinh(aL)
C A0 cosh(aL)

C A cosh a  z  cosh(aL)  sinh(aL) sinh a  z 


C A0 cosh(aL)
C A cosh a  z  cosh(a  L)  sinh(a  L) sinh a  z 

C A0 cosh(a  L)

C A cosh a  z  a  L 

C A0 cosh(a  L)

C A cosh ( / L)(  L)   


C A0 cosh(aL)
C A cosh      

C A0 cosh(aL)

C A cosh 1   

C A0 cosh( )
Ejemplos

Estime el coeficiente de transferencia de masa kc para el caso de la película de

líquido de la torre de absorción

WA 
4 D AB vmax
WL C A0  0  kc AC A
L

Por tanto,
WA 4 D AB vmax
kc 
AC A 
4 D AB vmax
WL C A0  0 WL C A0  0  L
L

Por definición: kc L 4 DAB vmax L

Sh  
D AB L D AB
 4 D AB vmax L2 4vmax L 4  Lv max    
     
L D AB 2  D AB     D AB 

Sh 
4
Re Sc1/ 2

 BSL
Ley Fick
Ecuación de balance (conducción + convección)
Nooo Tres casos: binaria D kte, binaria cD kte, Binaria vel. Cero (pag. 584)

Ejemplos:
Nooo Difusión a través de un líquido estancado – absorción? (pag. 545) No se puede por Ec
Difusión con reacción química homogénea (pag. 554) Si se puede
Difusión a través de un líquido – absorción (pag. 558) Si se puede
Deshumidificador (pag. 602) No se puede por Ec. Balance

Nooo Multicomponentes (pag. 586)

Análisis dimensional – Schmidt y Grashof (pag. 600)
Nooo Multivarialbes (pag. 612)
Nooo Capa límite (pag. 633)
Nooo Transporte en flujo turbulento (pag. 657)
Coeficientes de transferencia de masa (pag. 676)
Treybal

Tipos OPUS (sol- liq, liq-liq, etc)

Coeficientes de transferencia de masa (pag. 50)
Laminar (Pag. 55)
Turbulento (pag. 60)
Analogía masa y calor (75)
Casos sencillos (82)

Difusividad en sólidos (difusividad efectiva)

 Ecuación de balance de masa por componente [1]

Acumula = Entra – Sale  Genera

másica

molar

Para mezclas binarias o mezclas

Velocidad molar
promedio de peso molecular similar
v  v*
Velocidad másica
promedio

[1] Bayron BIRD, Warren STEWART y Edwin LIGHTFOOT. Transport Phenomena. 2da Edición. Apendice B
 sinh(a  L)
 cosha  z   sinha  z 
CA
C A0 cosh(a  L)

C A cosha  z cosh(a  L)  sinh(a  L) sinha  z 


C A0 cosh(a  L)

C A coshaL  az) cosha( L  z )  L

 
C A0 cosh(a  L) cosh(a  L) L

 kL2   
cosh 1   
z
  z   D AB  L  
cosh aL1    CA  
CA   L  
 C A0  kL2 
C A0 cosh(a  L) cosh 
 D AB 
 
 . . . .. . .
. . .. . …
.. . . .. . .
. . . .. .

Solución de amina
. . . .. . .

Pared Tubo interior

. . . . .. . .
Gas Natural
. . . .. . .
L (H2S)
. . . . .. . .
. . . .. .
. . . .. . .
. . . . .. . .
. . . .. .
. . . .. . .

 A se difunde en B y
reacciona dentro de B
A + B  AB
Reacción de primer orden rA = -kCA

CA0

L
Método de Sustitución

 2C A C A
DAB  vmax
x 2
z
 2C A C A
DAB  C A  f ( x, z )
 z v max 
1
x 2

x
Redefino:
 C A  f ( )
4 DAB ( z / vmax )

   
2
1 1
    A
x 4 DAB ( z / v)  x  4 DAB ( z / vmax )

 x
 B
 ( z / v) 2 4 D ( z / v)3 / 2
Por regla de la cadena:
C A dC A 
 
 ( z / vmax ) d  ( z / vmax ) 2

C A dC A 
  3
x d x

 2C A   C A    dC A   d 2C A   
2
        4
x 2
x  x  x  d x  d 2  x 
4 en 1

 2C A C A d 2C A   
2
C A
D  DAB    5
x 2
z v  d  x 
2
 ( z / vmax )
2 en 5

d C A    C A 
2 2
dC A
DAB      6
d  x 
2
 ( z / vmax ) d  ( z / vmax )
A B

Derivadas ordinarias

A en 6
d C A    
2 2
dC A
DAB    
d  x 
2
d  ( z / vmax )
d 2C A 1 dC A 
DAB   
d 4 D( z / v)
2
d  ( z / vmax )

B en 6
d 2C A 1 dC A x
DAB    
d 4 DAB ( z / vmax )
2
d 2 4 DAB ( z / vmax ) 3 / 2
 d 2C A 1 dC A x
DAB    
d 2 4 DAB ( z / vmax ) d 2 4 DAB ( z / vmax ) 3 / 2

Agrupando y simplificando
d 2C A dC A 4 DAB ( z / vmax ) x
  
d 2
d 2 4 D AB ( z / vmax ) 3 / 2 DAB

d 2C A dC A x
  
d 2
d D AB z / vmax
x

4 DAB ( z / vmax )
d 2C A dC A
   2
d 2
d
d 2C A dC A
 2  0
d 2
d
Ecuación diferencial ordinaria de 2do orden