UNIVERSIDAD DEL ISTMO

CAMPUS TEHUANTEPEC

TRANSFERENCIA DE MASA

Introducción a la Transferencia de Masa y

Transferencia de masa por difusión molecular

Lectura temática del libro: PRINCIPIOS DE TRANSFERENCIA DE MASA. Segunda Edición. Ricardo A.
Lobo Oehmichen. Universidad Autónoma Metropolitana Unidad Iztapalapa.

ELABORADO POR: LÓPEZ MARTÍNEZ FERNANDA

PROFESOR: M. EN C. SÁNCHEZ LÓPEZ JOSÉ LUIS

CARRERA: INGENIERIA QUÍMICA SEMESTRE: 5TO

A 28 de octubre de 2020. Juchitán de Zaragoza, Oaxaca.

El transporte de masa puede darse de dos maneras: Por difusión o por convección.
La difusión es causada por un gradiente de concentración mientras que la convección es originada
por el movimiento masivo de fluido. Se pueden dar dos tipos de convección, convección natural y
convección forzada, la convección forzada es causada por una fuerza externa mientras que la
convección natural se da debido a las diferencias de densidad causadas por las diferencias de
concentración o temperatura en el fluido.
La convección natural y la difusión siempre ocurren de manera simultánea, sin embargo, hay
factores que pueden hacer que la contribución por convección natural sea despreciable en
comparación con la difusión. La magnitud de la concentración del soluto y de la diferencia de
concentraciones son estos factores que juegan un papel importante, ya que en soluciones diluidas y
de densidad constante la convección natural es de magnitud insignificante comparada con la
difusión.
Descripción cuantitativa de la difusión
Para describir el proceso de difusión podemos analizarlo de dos formas, mediante una
representación del mecanismo en el cual el movimiento de los átomos o las moléculas es una
consideración básica o podemos ver a nuestro fluido como un continuo (estudiarlo
macroscópicamente) ignorando el mecanismo real por el cual se transportan las moléculas.
EL segundo caso es el que nos interesa como ingenieros químicos y nos lleva a considerar un tipo
de modelo fenomenológico, que es un modelo en el que solo importa la cuantificación del
fenómeno y no sus causas fundamentales. La descripción cuantitativa de los fenómenos nos lleva a
la elaboración de modelos matemáticos que relacionen las cantidades que deseamos conocer, con
variables y propiedades del sistema que podemos medir y calcular.

Ejemplo: Experimento para medir la rapidez de difusión.

Descripción del modelo: sistema a Temperatura y presión constantes. Tenemos dos recipientes
conectados mediante un capilar largo y una válvula que impide el paso del contenido de un
recipiente al otro, uno de los recipientes contiene oxigeno y el otro nitrógeno. Deseamos calcular la
rapidez de difusión que tendría el sistema una vez que abrimos la válvula.
Si tomamos mediciones de la concentración de O2 en el recipiente que originalmente no lo
contenía, veremos que a través del tiempo la concentración de O2 varia linealmente.
Mediante la pendiente de la línea podemos obtener la taza de cambio de la concentración respecto
al tiempo:
moles de O 2 transferidos CO 2 Vbulbo
( Flux difusivo de oxigeno )= =
tiempo t
El área del capilar por donde pasa el fluido puede causar variaciones en las mediciones, es
conveniente que nuestro experimento no dependa del área del capilar. Por tanto, se dividirá el flujo
de O2 por el área seccional del capilar. Obtenemos una nueva cantidad denominada densidad de
flujo o flux.
moles de O 2 cantidad transferidos CO 2 Vbulbo
( Flux difusivo de oxigeno )= =
0 densidad de flujo difusivo (tiempo)( Areadel capilar ) tAs

Podemos suponer que el flux difusivo depende de la diferencia de concentraciones entre los dos
recipientes.

Flux difusivo de O2=K (Diferencia de concentracion de O 2)

Además, el oxígeno debe transportarse a través del capilar desde un recipiente a otro. Entre mayor
sea la longitud L del capilar, más tardará el oxígeno en difundirse. Podemos suponer entonces que
el flux difusivo dependería de manera inversa de la longitud del capilar (dependencia lineal) de
modo que nuestra hipótesis sería que el flux difusivo es inversamente proporcional a la distancia a
lo largo de la cual ocurre la difusión.
Diferencia de concentracion de O 2 ∆ CO 2
Flux difusivo de O2=D =D
Longitud del capilar L
Donde a D se le denomina coeficiente de difusión o difusividad y procede de la llamada ley de Fick.
Ley de Fick
Una manera fenomenológica de cuantificar la difusión es representada por la ley de Fick
∂ CA
J A ,Z =−D AB
∂Z
Esta ecuación que representa la ley de Fick es aplicable a soluciones diluidas de densidad constante,
nos Permite conocer el flux difusivo en cada punto dentro del sistema, relaciona el flux difusivo con
el gradiente de concentración mediante el coeficiente de difusión D AB ( propiedad del sistema).

Podemos establecer el significado físico de la ley de Fick en los siguientes puntos:

1. Existe una tendencia natural (tendencia que depende del coeficiente de difusión) en la que el
soluto se difunde de una región con concentración menor a una con mayor concentración. 2. la
magnitud del flux difusivo es directamente proporcional a la magnitud del gradiente de
concentración 3. El flux difusivo y el gradiente de concentración tienen direcciones (signos)
opuestas.

El coeficiente de transferencia de masa

Existen problemas en donde no es posible dar información detallada sobre el flux o el perfil de
concentración, un ejemplo de este tipo de problemas es la disolución de pequeñas partículas de un
sólido A en agua (B).

Dentro del tanque que contiene el sólido y el

agua, se lleva también a cabo una agitación
provocando corrientes que pasan por encima de
cada partícula de A y arrastran al soluto
disuelto que se halla alrededor de ésta para
mezclarlo con el resto del líquido del tanque.
En este caso los mecanismos de transporte de
masa presentes son 3: la difusión, la
convección natural y la convección forzada.

El problema es complejo porque: (a) no conocemos la distancia exacta sobre la cual ocurre el
transporte de masa; (b) no conocemos los valores de la velocidad del líquido en las cercanías de
cada partícula del sólido, entonces para resolver este tipo de problemas es útil emplear el modelo
que usa coeficientes de transferencia de masa. Este modelo integra una constante de
proporcionalidad llamada coeficiente de transferencia de masa que engloba los diversos efectos que
no conocemos, este coeficiente solo depende de la agitación y de la viscosidad de la solución.
La cantidad total del sólido que se disuelve por unidad de tiempo (WA) se calcula suponiendo que
es
proporcional a la diferencia de concentración del soluto (A) que existe entre el líquido que toca la
superficie del sólido y la concentración promedio en el seno del líquido (CA1−CA2) y que es
también proporcional al área de las partículas sólidas (AS).

WA=KcAs ∆ CA =Kc ( CA 1−CA 2 ) As

Si dividimos entre el área de las partículas obtendremos el flux total de moles (NA) del soluto A
que se disuelve.
WA
=NA=Kc ∆ CA =Kc ( CA 1−CA 2 )
As
Fenómenos de transporte
Los fenómenos de transporte son fenómenos que pueden ser tratados de manera similar a como
tratamos el transporte de masa, es decir, podemos establecer un tratamiento unificado para los
fenómenos de transporte de momento, de energía y masa. Esto es así porque las leyes que los rigen
presentan en buena medida analogías entre sí, estas leyes son:
∂vy
τ zy =−μ … Ley de la viscosidad de newton(transporte de momentum)
∂z
∂T
q z =−k T … Ley de Fourier (Conduccion de calor)
∂z
∂C A
J A ,z =−D AB … Ley de Fick (difusion de masa)
∂z
descripción de las similitudes:

 τzy, qz y JA,z son los fluxes moleculares, de calor y masa.

 Las propiedades de los materiales están dadas por: μ (viscosidad), kT (conductividad
térmica) y DAB (coeficiente de difusión
 Los fluxes son provocados por un gradiente de potencial. El gradiente de velocidad origina
un flux de momentum, el gradiente de temperatura provoca un flux de calor y el gradiente
de concentración un flux de masa.
A pesar de sus similitudes estas leyes aun no pueden llamarse análogas, es necesario que las
constantes de proporcionalidad tengan las mismas dimensiones y que los gradientes que impulsan el
transporte expresen la misma idea de concentración.
Si multiplicamos y dividimos la expresión de la ley de viscosidad de newton por la densidad, y la
expresión para la ley de Fick por ρ C´ P, llegamos a las expresiones de los fluxes de momentum y
calor que sí son análogas a la ley de Fick:

−μ ∂ ( ρ v y ) ∂(ρ vy)
τ zy = =−v
ρ ∂z ∂z

−k T ∂ ( ρ C´ P T ) ∂ ( ρ C´ P T )
qz= =−∝
ρ C´ P ∂z ∂z
 Donde:

μ Difusividad de momentum longitud 2

v= =
ρ (
o viscosidad cinematica
[¿ ]
tiempo )
kT longitud 2
∝= =(difusividad termica)[¿ ]
ρ Ć P tiempo

Difusión en estado estacionario en una película estancada

Este caso de estudio aparentemente sencillo es de suma importancia porque describe el problema
típico de trasporte de masa y representa uno de los dos comportamientos limite de sistemas
difusivos. Nos interesa conocer el perfil de concentraciones, el flux difusivo y el flujo molar de un
soluto A dentro de un medio.
Pasos útiles para resolver este tipo de problemas:
1. descripción y comprensión de la situación física.
Sistema
tenemos un sistema en equilibrio, compuesto por dos soluciones perfectamente agitadas a cada
lado de una película delgada de volumen V hecha de un componente B de espesor Lz en
dirección z.
consideremos que la concentración del soluto en el lado izquierdo dentro de la película
permanece constante y tiene un valor CA0 y del lado derecho el calor es CAL.
Tenemos un área transversal constante LxLy=As por la que el soluto A se difunde en dirección
z, desde la región con mayor concentración (ubicada en z =0) hacia la de menor concentración
(en z=Lz). El sistema se encuentra a temperatura y presión constantes.
2. Hacer un balance de masa del compuesto de interés en un elemento diferencial de
volumen representativo del sistema, con el fin de obtener la ecuación diferencial del
flux como función de la distancia.

Balance de masa:

Rapidez de generación del Rapidez de acumulación

Flujo molar (moles/tiempo) Flujo molar (moles/tiempo) del soluto A en
- + soluto A por reacción en el =
del soluto A que entra por del soluto A que sale por
el volumen As∆ z
volumen As∆ z
difusión en el plano z difusión en el plano z+∆ z

Donde: Δz es el espesor del elemento diferencial de volumen y

AS es el área perpendicular a la dirección de la difusión.
Consideraciones: como no existe reacción química, podemos
despreciar el termino de generación. Adicionalmente el sistema
se encuentra en estado estacionario, por lo que el termino de
acumulación del Souto A es también cero.
Después de hacer estas simplificaciones y escribiendo en términos del flux tenemos:

moles moles
( A s J A , Z )|z− ( A s J A , Z )|z+∆ z=0 [ ¿ ] area ( area∗tiempo
=
tiempo )( )
Si dividimos entre el volumen de control (As∆ z ) y rearreglamos términos, obtenemos:

( A s J A ,Z )|z +∆ z− ( A s J A ,Z )|z J A , Z )|z+ ∆ z − ( J A , Z )|z

− [ As ∆ z ]
=−( (
∆z ) =0

Si tomamos el límite cuando Δz→0, obtenemos la derivada del flux difusivo de A con respecto a la
distancia:

J A , Z )|z+∆ z− ( J A , Z )|z −d J A , Z
lim ¿−
∆ z →0
(( ∆z )
=
dz
=0

La integración de esta expresión nos lleva a J A ,Z =cte , lo que nos dice que el flux difusivo es
constante e independiente de la distancia.

3. Emplear de la ley de Fick para relacionar el flux difusivo con la concentración, y aplicar las
condiciones de frontera para obtener el perfil de concentraciones.

Podemos obtener la ecuación gobernante del proceso al sustituir la derivada del flux de A
respecto a z en la ley de Fick.
−d J A , Z ∂ CA
=−D AB
dz ∂z
Al reordenar tenemos:

d 2 CA
=0 …(Ec . Diferencial ordinaria)
dz 2
Esta ecuación esta sujeta a las siguientes condiciones frontera:

CF1 En z=0 CA =CA 0

CF2 En z=Lz CA =CAL

Si integramos dos veces y evaluamos las constantes de integración en las CF. obtenemos el perfil
de concentraciones del soluto A dentro de la película.

z dCA −CA 0−CAL

CA =CA 0−( CAO−CAL ) ∴ =
Lz dz Lz
Podemos ver que la concentración varía linealmente con la distancia, el signo negativo en el
gradiente de concentración señala que la concentración decrece conforme aumenta la distancia.

Podemos escribir la ecuación de forma adimensional como:

CA−CAL Z d CA ¿
=CA =1− =1−z ∴
¿ ¿
=−1
CA 0−CAL Lz d z¿

Perfiles de concentración dimensional y adimensional de la difusión

en estado estacionario en una película estancada.

4. Obtener del flux y el flujo difusivos mediante la derivación del perfil de concentraciones.

Para

obtener el flux difusivo sustituimos la ecuación para el perfil de concentraciones en la

definición de la ley de Fick y realizamos la derivación indicada:
Ahora ya podemos cuantificar el flux difusivo en términos de variables conocidas y propiedades
del sistema, notamos también que es positivo, lo cual quiere decir que su dirección es la misma
que la de la coordenada z.

Otra cantidad importante es el flujo molar (o másico) del soluto que entra o sale del sistema, el
flujo molar de la película que estamos analizando es:

5. Obtención de la concentración, del flux y del flujo másico promedios.

El valor
medio o

promedio de una función f en fenómenos de transporte puede referirse a:

(a) el valor medio o promedio de la función en el volumen del sistema, (b) el valor medio o
promedio de la función sobre determinada área del sistema; y (c) el valor medio o
promedio de la función en un cierto intervalo de tiempo. Los valores promedio de los
incisos (a) y (b) se definen por:

Donde <f> es el valor promedio de la función f, V es el volumen sobre el cual se

promedia la función; dV es el elemento diferencial de volumen, A es el área sobre
la cual se promedia la función y dA es el elemento diferencial de área.

El promedio volumétrico sirve para encontrar la concentración promedio del soluto A en la

película y se define de la siguiente manera:

Comencemos por observar que el perfil de concentraciones no depende de x ni de y. El volumen

total de la película es V= LxLyLz = AsLz donde As = LxLy, y el elemento diferencial de volumen es
dV=ASdz. Entonces, usando la definición general del promedio volumétrico:

Sustituimos aqui la Ec. Para el perfil de concentraciones

y realizamos la integración.
Este resultado (el promedio aritmético entre el valor más alto y el más bajo de la concentración) se
debe a que el perfil de concentración es lineal y éste lo es porque el área AS a través de la cual
ocurre la difusión es constante para cualquier valor de z.

Dado que los valores en el flux y el flujo másico son constantes e independientes de z, los valores
promedio del flux difusivo y el flujo molar serán los mismos.

Difusión en estado estacionario a través de una película cilíndrica

descripción del caso de estudio:

una gota esférica de radio R1 de una sustancia A se evapora lentamente en nitrógeno (B).

Supongamos que la rapidez de evaporación de la gota está

controlada por la difusión de A en una película estancada de
B de espesor R2−R1.

La concentración de la sustancia A como vapor en la

superficie de la gota es CA1 y fuera de la película es CA2.
Consideremos que la concentración de A en la película es
baja, de modo que podemos suponer que el compuesto está
diluido, y que la evaporación ocurre en estado estacionario.

Para la resolución del problema se sigue una serie de pasos análogos a los que anteriormente se
presentaron, con la diferencia de que ahora se presenta una geometría esférica, es decir, El
elemento diferencial de volumen será un cascarón esférico de espesor Δr y de volumen 4πr^2Δr.

Balance de masa

Al dividir entre el elemento diferencial de volumen (4πr^2Δr) y tomando el límite cuando ∆ r → 0,

obtenemos la ecuación diferencial de (JA⋅r2) con respecto al radio.

Al integrar esta ecuación nos queda:

La ecuación señala que el flux difusivo JA,r es una función inversa cuadrada de la variable
independiente r, y no una constante.

Como paso siguiente sustituimos esta ecuación en la ley de Fick en coordenadas esféricas para
difusión radial ( ∂CA ) para relacionar la concentración
J A ,r =−DAB
∂r
con la variable independiente r.

Al integrar nos queda:

Utilizamos las siguientes condiciones de frontera

Para evaluar las constantes de integración y obtener el perfil de concentraciones, que también
podemos escribir en su forma adimensional.

La ecuación del flux difusivo se obtiene derivando la Ec. Para el perfil de concentraciones y
sustituyendo el resultado en la ley de Fick.

La rapidez de evaporación del sólido A es el producto del flux

difusivo multiplicado por el área evaluado en la superficie del sólido:

la concentración promedio en la película

esférica es:
Difusión en una película con disolución del soluto en las fronteras

A. Contacto entre fases que contienen al

menos dos componentes, uno de
ellos en común.
Descripción del sistema:
Tenemos una película de polímero B
que esta separando dos soluciones,
ambas contienen un soluto a en el
disolvente C. El polímero es insoluble
en las soluciones, el soluto A es poco
soluble en el polímero y el solvente C
es insoluble en el polímero.

Las soluciones están siendo agitadas de forma que las concentraciones permanecen uniformes en
todo el volumen.

Consideraciones: CA1 > CA2 , el sistema está en estado estacionario y a temperatura y presión
constante. Este problema se aborda con la metodología ya conocida, de hecho, nos plantea una
geometría y condiciones iniciales de las concentraciones anteriormente trabajadas, por tanto la
ecuación gobernante del sistema es:

d 2 CA
=0
dz 2
CA se refiere a la concentración del soluto A dentro del polímero, ya que el balance de materia lo
realizamos dentro de la película.

La diferencia con el caso anteriormente estudiado es que para que el soluto A pase de la solución
en el compartimento (1) al polímero B, primero debe disolverse en él y entonces, una vez que se
ha disuelto, comenzará a difundirse. De modo semejante, para que el soluto A salga del polímero
B a la solución en el compartimento (2), debe primero disolverse en ella.

Tenemos entonces el problema de conocer el valor de las concentraciones del soluto en las
fronteras de la película, CA0 y CAL. Para resolver este problema es necesario relacionar la
concentración del soluto A en las fronteras del polímero B, que desconocemos, con la
concentración del mismo soluto fuera de él (en las soluciones), que sí conocemos. La forma de
obtener esta relación es formular la hipótesis de que en las interfases solución-polímero, y sólo en
las interfases, existe equilibrio3. Si adoptamos esta hipótesis, podemos escribir las condiciones de
frontera de nuestro problema:
donde m es el coeficiente de partición, que es una propiedad de equilibrio del sistema polímero-
solución y que puede obtenerse a partir de curvas de equilibrio que se obtienen de datos
experimentales o que se calculan teóricamente a partir
de ecuaciones de estado o por los modernos métodos de
simulación molecular.

Entonces, si conocemos m y las concentraciones en las soluciones, la hipótesis de existencia de

equilibrio en la interfase nos permite conocer las
concentraciones en las fronteras del polímero.

Conociendo las concentraciones en la

frontera podemos obtener el perfil de concentraciones la
manera en que lo hemos venido haciendo, el perfil esta dado por:

Sustituyendo en la ley de Fick obtenemos la

expresión del flux difusivo:

Donde el grupo mDAB/Lz =PAB con frecuencia se denomina la permeabilidad del polímero al
soluto A. Como es evidente, entre mayor sea la permeabilidad, mayor será el flux difusivo.

Multiplicamos por el área transversal a la difusión para obtener la ecuación del flujo molar:

Contacto entre un compuesto puro y una fase con al menos 2 componentes.

Planteemos en mismo sistema, con la diferencia de que ahora el líquido en el compartimiento

es A puro. La condición de equilibrio en la interfase se seguirá cumpliendo de modo que el
único problema ahora es encontrar la concentración en z=0 dentro del polímero que estará en
equilibrio con el compuesto puro A.

Para que un compuesto puro y una solución de este se encuentren en equilibrio es necesario
que coexistan las dos fases a la misma temperatura y presión, y que el compuesto puro ya no
pueda disolverse más en la solución. Entonces se dice que la
solución está saturada y la concentración del soluto presente en
ella se denomina concentración de saturación.

Sabiendo esto podemos escribir entonces las condiciones

fronteras:
Integramos y evaluamos las constantes con estas condiciones frontera para obtener:

Con el estudio de los sistemas anteriores de transferencia de masa, podemos darnos cuenta
de que hay una base común en el procedimiento para su solución, variando algunas
condiciones o la geometría del caso.

Podemos tener incluso sistemas difusivos que incluyan reacciones químicas, aquí la difusión
de un compuesto ocurre simultáneamente con una reacción química en la que dicho
compuesto participa, es decir, la difusión y la reacción ocurren en paralelo. La diferencia de
estos sistemas con los sistemas difusivos en los que no ocurren reacción química es que, en
los primeros, el balance de masa ya incluye el término de rapidez de generación por reacción
química y por tanto requieren tratamientos enfocados a la presencia de la reacción.

Es importante mencionar que en ciertos procesos de transferencia de masa debemos

simplificar el modelo para su resolución, esto puede hacerse aplicando un análisis de orden de
magnitud a la ecuación diferencial gobernante. En este análisis tratamos de estimar el valor
promedio de una función en una cierta región de interés con la finalidad de descartar los
términos cuya aportación sea insignificante en comparación con el de otro termino en la
ecuación diferencial.

Por último es debido hacer énfasis en la importancia que

tiene establecer las propiedades en un sistema
difusivo, entender su comportamiento y las leyes
físicas que lo rigen.