Unidad Ii - Lobo

PRINCIPIOS DE TRANSFERENCIA DE MASA/ 2a.
EDICIÓN
6
TRANSFERENCIA DE MASA
EN LA INTERFASE
En los capítulos anteriores hemos analizado con algún detenimiento los fenómenos de
transporte de masa que ocurren dentro de una fase que se encontraba estacionario o en flujo
laminar. Para múltiples aplicaciones en ingeniería, sin embargo, es necesario relacionar esta
información con la rapidez de la transferencia de especies químicas de una fase a otra, la
cual ocurre a través de una interfase* . El contacto entre dos fases es frecuente en
operaciones industriales; tiene por objeto transferir uno o más compuestos de una fase a
otra y a menudo se da en condiciones en las que al menos una de las fases fluye. A la
transferencia de masa desde una interfase hacia una fase en movimiento o entre dos fases
inmiscibles se le denomina transferencia de masa convectiva.
* Una interfase es la superficie frontera entre dos fases. En general, las propiedades de cada fase en las
regiones muy cercanas a la interfase difieren de las que prevalecen en su seno.
R. LOBO/2007
266 PRINCIPIOS DE TRANSFERENCIA DE MASA/ 2a. EDICIÓN
Así, por ejemplo, la operación de secado (interfase sólido-gas) transfiere un compuesto

presente en el sólido hacia un gas que fluye sobre él; la adsorción ( interfase sólido-gas o
sólido-líquido) tiene por objeto retener en un sólido una o varias especies presentes en el
gas o en el líquido, los cuales generalmente están movimiento; la absorción (interfase gas-
líquido) un flujo de gas transfiere uno o varios compuestos hacia un flujo de un líquido; la
extracción (interfase líquido-líquido) transfiere compuestos de una fase líquida y los
concentra en otra fase líquida inmiscible; etc. La transferencia de masa entre fases también
puede tener por objeto promover una reacción química, como cuando se tranfiere un
compuesto presente en un gas hacia la superficie de un catalizador sólido, donde el
compuesto se transformará químicamente en otro.
Generalmente la transferencia de masa que ocurre en las interfases y las condiciones a las
que se lleva a cabo no pueden ser observadas experimentalmente con facilidad. Ello nos
obliga a echar mano de modelos empíricos aproximados como el del coeficiente de
transferencia de masa, introducido en el Capítulo 1, para tratar de cuantificar la rapidez con
la que los compuestos se transfieren, entender cómo esta transferencia se lleva a cabo y
cuáles variables le afectan, De hecho, el estudiante recordará que en los capítulos 2, y 3
hemos usado el modelo del coeficiente de transferencia de masa como condición de
frontera, pero no hemos profundizado en la naturaleza de este coeficiente.
Quizá la mayor dificultad para una descripción rigurosa de la transferencia de masa en las
interfases proviene de las condiciones de flujo que prevalecen en su vecindad. Cuando el
flujo sobre la interfase es laminar, estas condiciones son bien conocidas y el mecanismo de
transporte de masa que actúa es la difusión. En esta situación, es posible obtener una
descripción fenomenológica más o menos rigurosa de la difusión y la convección forzada;
así lo hicimos en los casos de la absorción de un gas en una película líquida descendente, y
en el de la disolución de un sólido hacia una película que analizamos en las Sec. 4.6 y 4.7,
respectivamente. Sin embargo, como bien sabemos, el mecanismo de transporte de masa
difusivo es lento. Con frecuencia se requiere aumentar la rapidez de la transferencia de
masa, lo cual demanda que los fluidos se muevan en flujo turbulento. Cuando la
transferencia de masa ocurre en interfases en cuya vecindad prevalecen condiciones de
flujo turbulento, los fenómenos a analizar se complican sustancialmente y la descripción
matemática del flujo es muy compleja o imposible de obtener. Es en estos casos cuando el
modelo del coeficiente de transferencia de masa resulta de mayor utilidad.
Aunque apenas nos iniciamos en el estudio de la transferencia de masa convectiva, es

necesario decir que ésta guarda muchas semejanzas con la transferencia de calor
convectiva. Seguramente el estudiante recuerda que para cuantificar la transferencia de
calor entre una superficie sólida y un fluido en movimiento se utiliza la llamada ley de
enfriamiento de Newton, que usa un coeficiente de transferencia de calor (h) y una
diferencia de temperatura. Sabemos que este coeficiente y, como veremos, también el de
transferencia de masa, son función de las características del flujo del fluido, de las
propiedades de éste y de la geometría del sistema. De modo que, aunque el tema que
iniciamos es "nuevo", tenemos muchos antecedentes que nos facilitan su estudio.
R. LOBO/2007
PRINCIPIOS DE TRANSFERENCIA DE MASA/ 2a. EDICIÓN 267
En este capítulo y el siguiente estaremos analizando problemas en los que la transferencia

de masa ocurre por difusión y convección forzada. Aquí nos concentraremos en los
aspectos más prácticos, que nos ayudan en los cálculos de la transferencia de masa
convectiva; en el Capítulo 7 estudiaremos los aspectos teóricos, que nos auxilian para
comprender más a fondo los fenómenos que ocurren en ella. Así, en la Sec. 6.1
introducimos una definición empírica del coeficiente de transferencia de masa. En la Sec.
6.2 veremos que este coeficiente empírico tiene un equivalente preciso en problemas en los
que usamos la ley de Fick para cuantificar la transferencia de masa que ocurre por difusión
y convección en un medio estacionario y/o en flujo laminar. Para contrastar los
mecanismos de transferencia de masa en flujo laminar y en flujo turbulento, en la Sec. 6.3
haremos una descripción cualitativa de la transferencia de masa en una interfase fluido-
sólido en cuya vecindad se dan condiciones de flujo turbulento. Enseguida estudiaremos
dos métodos para obtener los valores del coeficientes de transferencia de masa para
situaciones determinadas: en la Sec. 6.4 veremos cómo mediante experimentos y el análisis
dimensional es factible obtener correlaciones que permiten el cálculo del coeficiente de
transferencia de masa; y en la Sec. 6.5 lo haremos a partir de las denominadas analogías
entre las transferencia de momentum, calor y masa. En estas secciones analizaremos
únicamente la transferencia de masa desde una interfase hacia una fase fluida en
movimiento. Finalmente, en la Sec. 6.6 extenderemos nuestro análisis a situaciones de
mucha importancia práctica para el cálculo de la transferencia de masa que ocurre entre dos
fases.
6.1 Definición del coeficiente de transferencia de masa en una fase

Para definir el coeficiente de transferencia de masa nos interesa relacionar de manera
sencilla el flujo de masa de un compuesto que sale de una interfase hacia un fluido en
movimiento con características del sistema. Esperaríamos que la rapidez con la que la masa
se transfiere desde la interfase al fluido dependa del área de la interfase y de una diferencia
de concentración en la vecindad de la interfase. Es común definir el factor de
proporcionalidad k (= coeficiente de transferencia de masa) entre el flujo de masa, el área y
la diferencia de concentración como:
WAi = k A∆CA (6.1-1)
En esta ecuación, WAi es el flujo molar del compuesto que sale de la interfase, A es un área
característica indefinida y ∆CA es una diferencia de concentración característica, también
indefinida. Es importante destacar que en tanto A y ∆CA no estén claramente definidos, k
tampoco lo estará. Esto es fundamental porque existen múltiples posibles definiciones de k,
según se seleccionen A y ∆CA.
Pasemos ahora a precisar la definición de k que adoptaremos en este texto. Refirámonos a

la Figura 6.1-1, donde hemos tomado una especie de "fotografía instantánea", y hemos
R. LOBO/2007
"congelado" un proceso típico de transferencia de masa interfacial. En este sistema existe un

sólido A poco soluble que se disuelve en contacto con un líquido bien agitado. La
concentración del soluto A en el líquido en inmediato contacto con el sólido es CAi y la
concentración promedio del mismo en el seno del líquido es CAb. El área de la interfase es
Ai. Entonces, podemos concretar nuestra definición del coeficiente de transferencia de
masa que utilizaremos:
(
WA i = kc Ai CA i − CA b ) (6.1-2)
Al dividir esta ecuación por el área interfacial Ai, obtenemos:
(
N A i = k c CA i − CA b ) (6.1-3)
interfase
solución bien
sólido A agitada
CAi
concentración
promedio
solución bien CAb

agitada
N Ai
fluctuaciones de
concentración
Figura 6.1-1. Transferencia de masa desde una interfase hacia una solución bien mezclada.
Observemos que ahora el coeficiente de transferencia de masa es kc, donde el subíndice c es

un recordatorio de que hemos decidido expresar la diferencia de concentración en unidades
molares. El flux total de A relativo a la superficie interfacial NAi sale de la interfase en una
dirección perpendicular a ésta. La concentración en el líquido en inmediato contacto con el
sólido CAi es aquella que se encuentra en equilibrio con el mismo compuesto en la otra fase
(hipótesis de equilibrio en la interfase); la concentración promedio1 CAb es aquélla "lejos"
1 Para el significado preciso de CAb debemos consultar el Apéndice C. Además, conviene recordar que en flujo
turbulento (como sería el caso de una solución bien agitada) existen pequeñas fluctuaciones al azar en el valor de
las propiedades del fluido y se dan alrededor de un valor promedio. En el resto del texto, cada vez que nos
refiramos al valor de una propiedad en el seno del fluido, estaremos aludiendo precisamente a este valor
R. LOBO/2007
de la interfase (fuera de la capa límite de concentración). Por otra parte, debemos señalar
que estamos considerando que la transferencia de masa ocurre en solución diluida. El
tratamiento general del tema se basará en este supuesto, aunque más adelante
comentaremos sobre la transferencia de masa en soluciones concentradas y/o con altas tasas
de transferencia de masa.
Es necesario advertir que esta sencilla definición requiere de algunas precisiones. En primer
lugar, la definición dada en la Ec.(6.1-3) no es aceptada universalmente. Como ya
apuntamos arriba, es frecuente encontrar en la literatura especializada de transferencia de
masa varias definiciones distintas a la dada aquí. A diferencia del coeficiente de difusión,
que tiene siempre las dimensiones [L2/t] sin importar la forma de la ley de Fick, el
coeficiente de transferencia de masa tendrá dimensiones que dependerán de la forma de
expresar la diferencia de concentraciones y el flux. Esto se hace evidente en la Tabla 4.1-1.
Tabla 4.1-1. Algunas definiciones del coeficiente de transferencia de masa
Definición Unidades del coeficiente de

N A [ = ] mol cm ⋅ s ; nA [ = ] g cm ⋅ s transferencia de masa
2 2
N A = kc ∆CA (6.1-4) cm s
N A = kG ∆pA (6.1-5) 2
mol cm ⋅ s ⋅ atm
N A = k y ∆y A ; N A = k x ∆ x A (6.1-6) mol cm ⋅ s
2
nA = k ρ ∆ρ A (6.1-7) cm s
La definición dada por la Ec.(6.1-4) es común, ya sea en un líquido o de un gas; es

frecuente que este coeficiente se use como kL en líquidos. La Ec.(6.1-5) se emplea
exclusivamente en gases, donde la fuerza motriz es una diferencia de presiones parciales
del compuesto que se transfiere. Cuando la fuerza impulsora de la transferencia se expresa
en fracción mol, se emplean las Ecs.(6.1-6), dependiendo de si se trata de un líquido (kx) o
de un gas (ky). Aunque no de empleo común, no es raro encontrar en la literatura el uso de
la Ec.(6.1-7) donde la diferencia de concentración se expresa en unidades másicas. La
utilización de las variadas definiciones pudiera conducir a confusión; sin embargo, todas
ellas expresan lo mismo y, una vez que se cuenta con un valor del coeficiente expresado en
cualesquiera de sus modalidades, es posible convertirlo a otra forma diferente mediante una
simple conversión de unidades. Esto se ilustrará en el Ejemplo 6.1-2.
La diferencia de concentración que impulsa la transferencia de masa (∆CA u otra) debe

definirse sólo después de una cuidadosa reflexión sobre la situación concreta. En mi
experiencia como profesor, he encontrado que los estudiantes, incluyendo los de posgrado,
promedio. Las propiedades en el seno del fluido se designarán con el subíndice b; y las de la interfase con el
subíndice i. En la Fig. 6.1-1 se muestra una posible distribución de concentraciones promedio en la cercanía de
una interfase, así como las fluctuaciones. Contra lo que pudiera pensarse, la intensidad de la turbulencia es
mayor muy cerca de la interfase que lejos de ella [Klebanoff, PS (1955) NACA Report 1247]. En la Sec. 7.6 se
analizan con más detalle las propiedades promedio en flujo turbulento.
R. LOBO/2007
tienen dificultades para definir correctamente esta diferencia de concentración,

particularmente cuando se trata de emplear el concepto en situaciones prácticas que
involucran equipos de proceso. Sin embargo, la solución a este problema es sencilla, con tal
que prestemos la debida atención. Cualquier proceso de transferencia de masa es dinámico,
cambia constantemente en el tiempo y en el espacio. Elaboremos una regla general para
identificar la diferencia de concentración que emplearemos en este texto:
La diferencia de concentración a utilizar es aquella que prevalece entre la superficie

desde donde se tranfiere el compuesto y la del fluido hacia donde ocurre la
transferencia.
Tratemos ahora de analizar el proceso dinámico de transferencia de masa con dos ejemplos
que ilustran situaciones extremas de variación de la diferencia de concentración en el tiempo
y en el espacio.
A CAi
CAb ( t )
CAb ( t )
líq. con concentración
uniforme CA b ( t )
Figura 6.1-2. Disolución de un sólido en un líquido. La diferencia de concentración [CAi−CAb(t)] varía en el

tiempo y no varía en el espacio.
Consideremos primero un sistema donde la diferencia de concentración varía en el tiempo,

pero no en el espacio, el cual se ilustra en la Figura 6.1-2. En esta figura representamos el
proceso de disolución de múltiples partículas de un sólido A puro que es poco soluble en un
líquido B. Para no complicar el experimento, supongamos que durante el tiempo que dura el
proceso la cantidad disuelta es pequeña y que por ello, el tamaño (y el área) de las partículas
permanece constante. El líquido en el tanque se encuentra bien agitado, por lo que podemos
suponer que la concentración promedio del soluto en el líquido, CAb, es uniforme, esto es, no
varía en el espacio. Sin embargo, debido a que el sólido se está disolviendo conforme pasa el
tiempo, CAb aumenta con tiempo, esto es, CAb =CAb(t). Por otra parte, si ponemos nuestra
atención en la interfase sólido-líquido y hacemos uso de la condición de equilibrio en la
interfase, veremos que CAi= CAsat mientras exista sólido A sin disolverse. En resumen, CAi
permanece constante y CAb aumenta con el tiempo. Si aplicamos la regla general enunciada
arriba, veremos que la diferencia de concentración a utilizar es [CAi−CAb(t)]. Esta diferencia
de concentración es la que existe en todo momento para impulsar la transferencia de masa
R. LOBO/2007
entre el sólido y el líquido. Dado que CAi es constante y CAb(t) aumenta con el tiempo, la
diferencia de concentración [CAi−CAb(t)] disminuye con el tiempo. Naturalmente, el flux
de disolución N A i = kc ( CA i − CA b ) será mayor a tiempos cortos que a tiempos largos.
gas libre de soluto

seno del líquido
líquido puro CAb ( z )
pAb ( z )
CAi ( z )
empaque L pAi ( z ) G
inerte
CAb ( z )
gas con soluto seno del gas

pAb ( z )
líquido con soluto
Figura 6.1-3. Diferencias de concentración en la absorción de un soluto. Las diferencias de concentración en

el gas y en el líquido cambian con la posición a lo largo de la distancia.
Examinemos ahora un caso en el que la diferencia de concentración varía en el espacio, pero

no en el tiempo. El sistema se ilustra en la Figura 6.1-3 y consiste de una columna de
absorción. Esquemáticamente, este equipo es un recipiente cilíndrico alto que está lleno de un
sólido inerte llamado empaque o relleno (estos dispositivos se llaman también columnas
empacadas o columnas de relleno). Un gas, que contiene un compuesto A, fluye
ascendentemente a lo largo de la columna, y transfiere este compuesto a un líquido que fluye
descendentemente. Conforme el líquido desciende, la concentración del soluto aumenta, en
tanto la concentración de A en el gas se reduce conforme éste asciende. El propósito del
empaque es provocar un área grande de contacto gas-líquido y así aumentar la cantidad
transferida de soluto. Para simplificar nuestra discusión consideremos que la absorción ocurre
en estado estacionario y que en la dirección radial de la columna las concentraciones del
soluto en gas y en el líquido son uniformes. En esta situación ocurre un doble proceso de
transferencia de masa en serie del soluto: (a) del seno del gas a la interfase con el líquido,
donde se disuelve, y (b) de la interfase del lado del líquido hasta el seno del líquido. Por el
momento, concentremos nuestra atención en el líquido. Conforme el líquido desciende,
ambas CAi y CAb cambian; entonces, tenemos una situación en la que CAi= CAi(z) y
CAb=CAb(z), y la fuerza motriz (CAi−CAb) del proceso de transferencia varía con la altura
de la columna. Algo semejante ocurre si miramos el proceso desde el gas. Como
consecuencia de la variación de las concentraciones a lo largo de la columna, la cantidad
de masa de A transferida de una fase a la otra también varía con la altura de la columna.
En conclusión, debemos usar diferencias de concentración locales y tomar en cuenta su
variación con la distancia z.
R. LOBO/2007
Un parámetro esencial en los problemas de transferencia de masa convectiva es el área de

la superficie interfacial, o simplemente área interfacial. Esta es el área a través de la cual
ocurre la transferencia de masa. Su importancia radica en que de ella depende la cantidad
de masa transferida: entre mayor sea el área interfacial mayor será la cantidad de masa que
se tranfiere, aunque el flux permanezca constante:
( desde la interfase )
Flujo molar de A transferido = W = N × área interfacial = N A
Ai Ai Ai i (6.1-8)
El área interfacial es un parámetro que el ingeniero puede controlar mediante una selección
juiciosa de las condiciones de operación, principalmente de los flujos de gas y líquido. Por
ello, es importante identificar y cuantificar el área de transferencia de masa. En algunos
casos, la identificación y cuantificación de esta área son relativamente sencillas porque son
evidentes de la geometría del problema. Por ejemplo, en el sistema ilustrado en la Figura
6.1-2 es claro que el área de transferencia de masa es la suma del área de las partículas
sólidas, ya que desde su superficie ocurre la transferencia. Sin embargo, existen muchos
sistemas en los que el área interfacial no es fácil de identificar y cuantificar. Pensemos en el
sistema de la Figura 6.1-3, donde es prácticamente imposible establecer el área de
transferencia entre el líquido y el gas, ya que éste burbujea caóticamente conforme asciende
por los huecos que dejan libres el empaque y el líquido que desciende. Sabemos que el área
de transferencia de masa es el área de contacto gas-líquido, pero en la situación descrita no
es posible saber cuál es el valor de esta área (es tentador, pero erróneo, pensar que el área
interfacial es el área del empaque sólido; ¿por qué es erróneo?). En casos como el de este
ejemplo, es común introducir un parámetro adicional, denominado área específica de
contacto gas-líquido, que se representa por a y se define por:
dAi área interfacial

a= = [ =] distancia −1 (6.1-9)
dV volumen del lecho empacado
Con ayuda de este parámetro1, podemos definir un flujo molar de A por unidad de volumen
de lecho:
( )
N A i a = ( kc a ) CA i − CA b [ = ]
moles de A transferidas
volumen de lecho ⋅ tiempo
(6.1-10)
1 La definición dada por la Ec,(6.1-9) se plantea en términos de los diferenciales de área interfacial y de
volumen del lecho, porque a puede variar con la altura del lecho fijo debido, entre otras causas, a variaciones
grandes en los flujos del gas y del líquido. Consideremos que si la cantidad de masa del soluto A transferida
del gas al líquido es grande, el flujo másico total de gas se reducirá apreciablemente conforme asciende en la
columna y el flujo de líquido aumentará conforme desciende. El parámetro a será aproximadamente constante
sólo cuando las cantidades de soluto transferidas sean pequeñas.
R. LOBO/2007
2
longitud longitud
kc a [ = ] ⋅ 3
= tiempo
−1
(6.1-11)
tiempo longitud
donde el producto (kca) recibe el nombre de coeficiente de transferencia de masa

volumétrico. En vista de que es difícil medir NAi y a por separado, se acostumbra medir el
producto el (NAia), lo cual es relativamente fácil. El área específica es un parámetro que
puede obtenerse experimentalmente (deseable, aunque no sencillo) o se puede estimar de
correlaciones que veremos más adelante. El área específica a depende de los flujos de gas y
líquido, de las propiedades de los fluidos y de las características físicas del empaque.
Podemos resumir nuestra discusión sobre la definición del coeficiente de transferencia de

masa en los siguientes puntos: (1) no definir un coeficiente de transferencia de masa, en tanto
no se defina el área interfacial y la diferencia de concentración a utilizar; (2) identificar la
interfase desde la cual ocurre la transferencia y la fase hacia la cual sucede; (3) identificar y
cuantificar el área interfacial a través de la cual ocurre la transferencia de masa; (4) identificar
la diferencia de concentración que impulsa la transferencia, prefiriendo las locales; (5) definir
el coeficiente de transferencia de masa de acuerdo a los puntos (2) a (5). como veremos, el
uso del modelo del coeficiente de transferencia de masa en conjunción con los balances de
masa son una poderosa herramienta para resolver importantes problemas prácticos.
Ejemplo 6.1-1. Cálculo de diferentes coeficientes de transferencia de masa
Un líquido A se evapora a 25 ºC y 1 atm de presión dentro de un tubo de 0.3 cm de

diámetro interior y cuyo extremo inferior se encuentra sellado. Por encima del extremo
abierto del tubo fluye una mezcla gaseosa A-B, donde la presión parcial de A es 0.01 atm.
Una vez que se ha estabilizado el proceso, se determina que la rapidez de evaporación del
líquido A es 3.92×10−5 mol/h. La presión de vapor del líquido A a la temperatura del
experimento es 0.15 atm.
Determinar los coeficientes de transferencia de masa kc, ky, y kG, suponiendo que la mezcla
gaseosa es diluida.
Solución. El sistema se encuentra representado en la Figura 4.4-1. Las definiciones de los

coeficientes de transferencia de masa se pueden tomar de la Tabla 6.1-1,
N A = kc ( CA1 − CA2 ) ; N A = k y ( yA1 − yA2 ) ; N A = kG ( pA1 − pA2 ) (6.1-12)
El flux total NA se puede obtener de dividir la rapidez de evaporación por el área de

transferencia de masa, que es la misma que el área seccional del tubo,
R. LOBO/2007
−5
WA WA 3.82 × 10 mol hr 1 hr −7
NA = = = 2
= 1.50 × 10 mol cm ⋅ s
Ai π dt2 4 2
π (0.3 cm) 4 3600 s
Las diferencias de concentración son las que prevalecen entre la interfase gas/líquido y el
extremo abierto del tubo, a saber
pA1 = pAsat = 0.15 atm ; pA2 = 0.01 atm
pAsat 0.15atm
CA1 = CAsat = = 3
−6
= 6.13 × 10 mol cm
3
RT 82.05 cm ⋅ atm mol ⋅ K × 298 K
pA2 0.01 −6
CA2 = = = 0.41 × 10 mol cm
3
RT 82.05 × 298
yA1 = pA1 P = 0.15 atm 1 atm = 0.15

yA2 = pA2 P = 0.01 atm 1 atm = 0.01
Sustituyendo los valores del flux y las concentraciones en las definiciones del coeficiente
de transferencia de masa,
−7 2
NA 1.50 × 10 mol cm ⋅ s −2
kc = = = 2.62 × 10 cm/s
CA1 − CA2 ( 6.13 − 0.41) × 10 −6 mol cm3
−7 2
NA 1.50 × 10 mol cm ⋅ s −6 2
ky = = = 1.07 × 10 mol cm ⋅ s
yA1 − yA2 ( 0.15-0.01)
−7 2
NA 1.50 × 10 mol cm ⋅ s −6
kG = = 2
= 1.07 × 10 mol cm ⋅ s ⋅ atm
pA1 − pA2 ( 0.15-0.01) atm
Ejemplo 6.1-2. Conversión entre diferentes coeficientes de transferencia de

masa
Se cuenta con valores del coeficiente de transferencia de masa de un gas A cuyo valor es
kc=0.093 cm/s. Asímismo, el mismo soluto se transfiere en agua con un coeficiente de
transferencia de masa kc=0.13 cm/s. La presión es de 1 atm y la temperatura 25 °C.
Calcular kG, ky, kL y kx.
R. LOBO/2007
Solución. Suponemos que la transferencia del gas se da del seno del gas a la interfase. Si
consideramos que el flux transferido es el mismo, independientemente del tipo de
coeficiente de transferencia de masa y de la fuerza motriz empleada:
( ) ( )
N A i = kc CA b − CA i = kG pA b − pA i = k y ( yAb − yA i ) (6.1-13)
De esta relación podemos obtener:
CA b − CA i kc pA b − pA i k
kG = kc = = c
pA b − pA i RT pA b − pA i RT
donde hemos usado el hecho de que para un gas CA = pA RT . Sustituyendo valores:
0.093 cm/s −6
kG = = 3.8 × 10
2
mol cm ⋅ atm ⋅ s
(82.05 atm ⋅ cm 3
/ mol K ) ( 298 K )
De la Ec.(6.1-13) también podemos obtener:
k y = kG
( pAi − pAb ) = k P ( yAi − yAb ) = k
GP
( yAi − yAb ) G
( yAi − yAb )
Sustituyendo valores
( )
k y = 3.8 × 10 −6 mol cm 2 ⋅ atm ⋅ s (1 atm ) = 3.8 × 10−6 mol cm 2 ⋅ s
En cuanto a los coeficientes de transferencia de masa del líquido, podemos tratarlos de

modo semejante a como lo hicimos con el gas sabiendo que el flux es el mismo:
( ) ( )
N A i = kc CA i − CA b = k L CA i − CA b = k x ( xAi − xAb ) (6.1-14)
En consecuencia,
kc = k L = 0.13 cm s
y
CA i − CA b xAi − xAb
kx = kL = k L CL = k L CL
xAi − xAb xAi − xAb
R. LOBO/2007
Aquí hemos hecho uso de CA = CL xA , donde CL es la concentración molar total de la

solución líquida. Ahora se nos presenta un pequeño problema, ya que no conocemos CL; sin
embargo, si consideramos que la solución líquida es diluida, podremos suponer que la
concentración molar total de la solución es aproximadamente igual a la del solvente, en este
caso agua, esto es
ρL ρ H 2O g 1 mol
CL = ≅ CH 2O = = 1.0 3
= 0.056 mol cm
3
ML M H 2O cm 18 g
Entonces, sustituyendo valores:
k x = ( 0.13 cm s ) ( 0.056 mol cm 3 ) = 7.28 × 10 −3 mol cm 2 ⋅ s
Como conclusión general de este ejemplo podemos establecer las relaciones entre algunas
de las definiciones del coeficiente de transferencia de masa:
RT k y
para gases: kc = kG RT = k y =
P CG
(6.1-15)
k
para líquidos: kc = k L = x
CL
6.2 Algunas expresiones teóricas del coeficiente de transferencia de masa

En la sección anterior introducimos una definición empírica del coeficiente de transferencia
de masa como una forma de calcular el flux de masa en una interfase. En esta sección
veremos que esta manera de cálculo del flux tiene equivalentes precisos con el flux que
puede calcularse mediante el uso de los balances de masa y la ley de Fick. Esta
equivalencia ocurre cuando el transporte de masa ocurre únicamente por difusión o por
difusión en flujo laminar. Ello permite encontrar expresiones teóricas del coeficiente de
transferencia de masa que ilustran las variables que afectan a este coeficiente y que son la
base de algunos de los modelos de la transferencia de masa convectiva que se estudiarán en
el siguiente capítulo.
En esta sección derivaremos expresiones del coeficiente de transferencia de masa para las
siguientes situaciones que ya estudiamos en otras secciones del texto: (a) difusión en una
película estancada en solución diluida (Sec. 2.1) y contradifusión equimolar (Sec. 4.5); (b)
difusión y convección en una película estancada en soluciones concentradas (Sec. 4.4); (c)
R. LOBO/2007
difusión en un medio seminfinito (Sec. 3.2); (d) difusión y convección forzada en flujo
laminar en la absorción de un soluto (Sec. 4.6) y; (e) difusión y convección forzada en flujo
laminar en la disolución de un sólido (Sec. 4.7). Lo que haremos será comparar las
expresiones integradas del flux para cada uno de las situaciones concretas con la definición
del coeficiente de transferencia de masa dada por la Ec.(6.1-1) y estableceremos la relación
pertinente. La discusión estará limitada a soluciones binarias isotérmicas.
Difusión en solución diluida y contradifusión equimolar en soluciones concentradas.

En ambos casos, sabemos que el flux total es igual al flux difusivo
N A,z = J A,z (6.2-1)
La expresión del flux difusivo JA para difusión en una película en solución diluida está dada
por la Ec.(2.1-14) y el de NA para contradifusión equimolar está dada por la Ec. (4.5-4) En
ambos casos, el flux es constante:
N A,z z =0 = J A,z z =0 =
DAB
Lz
( )
CA0 − CA L = N A i (6.2-2)
Si igualamos la Ec.(6.2-2) con la Ec.(6.1-1), podemos hacer las siguientes equivalencias
( CAi − CAb ) = ( CA0 − CA L ) (6.2-3)
DAB
kc = (6.2-4)
Lz
Como podemos apreciar en la Ec.(6.2-4), el coeficiente de transferencia de masa es

proporcional al coeficiente de difusión e inversamente proporcional a la longitud de la
trayectoria de difusión. Aunque sencillo, este resultado es muy importante, ya que
proporciona los elementos teóricos de la llamada teoría de la película de la transferencia de
masa convectiva que se abordará en el siguiente capítulo.
Difusión y convección en una película estancada en soluciones concentradas. La

situación física es la representada por la Figura 4.4-1. El flux total de A que se evapora de
la interfase líquido-gas es la dada por Ec.(4.4-10) y es constante para toda z
C DAB ⎛ 1 − yAL ⎞ C DAB ⎛ yBL ⎞

N A,0 = ln ⎜ ⎟= ln ⎜ ⎟ = N Ai (6.2-5)
Lz ⎝ 1 − yA0 ⎠ Lz ⎝ yB 0 ⎠
La Ec. (6.1-1) no es directamente comparable con la Ec.(6.2-5) ya que en ésta la diferencia

de concentración no es lineal. Sin embargo, si multiplicamos y dividimos la Ec.(6.2-5) por
R. LOBO/2007
⎡⎣(1 − yAL ) − (1 − yA 0 ) ⎤⎦
podemos obtener
C DAB (1 − yAL ) − (1 − yA0 ) ⎛ 1 − yAL ⎞

NA = ln ⎜ ⎟ (6.2-6)
Lz (1 − yAL ) − (1 − yA0 ) ⎝ 1 − yA0 ⎠
En la Ec.(6.2-6) podemos reconocer una diferencia de concentración media logarítmica,

que definiremos por
(1 − yAL ) − (1 − yA0 ) = yBL − yB0

(1 − yA )ln = ( yB )ln = (6.2-7)
⎛ 1 − yAL ⎞ ⎛y ⎞
ln ⎜ ⎟ ln ⎜ BL ⎟
⎝ 1 − yA0 ⎠ ⎝ yB0 ⎠
Entonces, usando la Ec. (6.2-7) podemos escribir la Ec. (6.2-6) de la siguiente manera:
C DAB
N A,0 = (
Lz ( yB )ln
yA0 − yA L = ) DAB
Lz ( yB )ln
(
CA0 − CA L ) (6.2-8)
Ahora sí podemos comparar las Ecs. (6.1-1) y (6.2-8), y es entonces evidente que
( CAi − CAb ) = ( CA0 − CA L ) (6.2-9)
DAB
kc = (6.2-10)
Lz ( yB )ln
La Ec.(6.2-10) es la relación buscada. Es importante notar que en este caso el coeficiente de

transferencia de masa depende de la diferencia de concentración media logarítmica
definida por la Ec.(6.2-7) debido a que la solución es concentrada y/o existen altas tasas de
transferencia de masa, de modo que el término yA ( N A + N B ) del flux total de A no puede
despreciarse. Ello quiere decir que si comparamos este coeficiente de transferencia de masa
con el obtenido para solución diluida en la Ec. (6.2-4), éste será menor que aquél. El
resultado expresado por la Ec.(6.2-10) es importante porque proporciona las bases teóricas
para calcular el coeficiente de transferencia de masa en situaciones generales, a partir del
coeficiente de transferencia de masa que se calcula para soluciones diluidas. En la Sec. 7.5
continuaremos el análisis de este efecto.
R. LOBO/2007
Difusión en un medio seminfinito en solución diluida (teoría de penetración). En este

caso, la situación física está representada en la Figura 3.2-2, y el flux en la interfase es la
función del tiempo dada por la Ec.(3.2-13)
J A, z z =0 =
DAB
πt
(
CA0 − CA ∞ = N A i ) (6.2-11)
Igualando esta ecuación con la Ec.(6.1-1), es evidente que
( CAi − CAb ) = ( CA0 − CA ∞ ) (6.2-12)
DAB
kc = (6.2-13)
πt
Observemos dos cosas. En primer lugar, ahora el coeficiente de transferencia de masa es

una función del tiempo, porque la transferencia de masa se encuentra en régimen transitorio
y no en estado estacionario como en los casos anteriores, esto es, es un coeficiente
instantáneo. En segundo lugar, notemos que cuando el tiempo es igual a cero, el valor del
coeficiente es infinito.
En este tipo de casos, con frecuencia es conveniente definir un flux y un coeficiente de

transferencia de masa promedio en el tiempo que dura la transferencia. El flux promedio
durante el tiempo tc está dado por Ec. (3.2-14)
( J A, ) z=0 = t1 ∫0 DAB
( )
tc
J A, z z =0 dt =2 CA0 − CA∞ = N A i (6.2-14)
z
c π tc
De esta expresión es evidente que el coeficiente de transferencia de masa promedio en el

tiempo de contacto tc es
DAB
kc = 2 (6.2-15)
π tc
La Ec.(6.2-15) es la relación buscada. Observemos que el coeficiente promedio tiene un

valor que es el doble del coeficiente instantáneo evaluado al mismo tiempo y es
proporcional a √DAB. Los resultados expresados en las Ecs.(6.2-13) y (6.2-15) son
importantes porque son la base de la teoría de penetración de la transferencia de masa
convectiva en las interfases fluido-fluido que se expondrá en el Capítulo 7, que a su vez son
de gran utilidad en las aplicaciones que involucran sistemas gas-líquido y líquido-líquido,
que se estudiarán en el Capítulo 7.
R. LOBO/2007
Difusión y convección forzada en flujo laminar en la absorción de un soluto (teoría de

penetración). En este caso, la situación física está ilustrada en la Figura (4.6-1) La
característica física importante del problema es que el soluto transferido desde la interfase
gas-líquido experimenta una velocidad constante igual a vmax en dirección perpendicular a
la dirección de la transferencia. El flux en la interfase está dado por la Ec.(4.6-13).
N A,0 =
DAB vmax
πy
(
CAsat − 0 = N A i) (6.2-16)
Si comparamos esta ecuación con la (6.1-1), será evidente que
( CAi − CAb ) = ( CAsat − 0 ) (6.2-17)
DAB vmax
kc = (6.2-18)
πy
Notemos que en este caso el coeficiente de transferencia de masa es una función de la

posición a lo largo de la película descendente, esto es, se trata de un coeficiente de
transferencia de masa local; observemos también que en y=0 el coeficiente tiene un valor
infinito y su valor respecto a longitud de la película líquida y siempre disminuye hasta
llegar a Ly.
Podemos escribir el coeficiente de transferencia de masa local definido por Ec.(6.2-18) en

términos de números adimensionales. Si multiplicamos ambos lados de la ecuación por
(y/DAB), multiplicamos y dividimos el cociente dentro del radical por la viscosidad
cinemática y sustituimos la velocidad máxima de la película líquida en términos de su
velocidad promedio 〈vy〉=(2/3)vmax,
( )
12 12
⎡ 2 ν⎤ ⎛ 3 y vy
kc y 1 ⎢ y 2 DAB 3 v y ν ⎞
= Sh y = ⋅ ⋅ ⎥ =⎜ ⋅ ⋅ ⎟
DAB π ⎢ DAB
2
y ν⎥ ⎜ 2π ν DAB ⎟
⎣ ⎦ ⎝ ⎠ (6.2-19)
kc y
( ) ( )
12 12
Sh y = = 0.691 Re y Sc = 0.691 Pem y
DAB
donde Shy es el número de Sherwood local, lo cual queda indicado con el subíndice y que
indica la distancia variable a lo largo de la película líquida. Notemos que el número de
Reynolds y el número de Peclet están basados en la velocidad promedio de la película
descendente y que también son de carácter local, es decir, varían con la distancia y.
R. LOBO/2007
A menudo es deseable obtener un valor promedio del coeficiente de transferencia de masa

sobre toda la distancia Ly. El problema se reduce a encontrar el flux promedio que penetra a
la película descendente sobre toda la superficie LxLy y comparar términos. El flux promedio
está dado por la Ec.(4.6-14)
N A,0 = 2
DAB vmax
π Ly
(
CAsat − 0 = N A i ) (6.2-20)
(22)
Es claro entonces que el coeficiente de transferencia de masa promedio sobre la distancia Ly
será
DAB vmax
kc = 2 (6.2-21)
π Ly
Observemos que el coeficiente promedio tiene un valor que es el doble del local evaluado a
la misma distancia [Ec.(6.2-18)]. El coeficiente de transferencia de masa promedio dado en
la Ec.(6.2-21) la podemos escribir en términos de números adimensionales, de modo
semejante a como lo acabamos de hacer con el coeficiente local. Si multiplicamos ambos
lados de la ecuación por (Lz/DAB), multiplicamos y dividimos el cociente dentro del radical
por la viscosidad cinemática y sustituimos la velocidad máxima de la película líquida en
términos de su velocidad promedio 〈vy〉=(2/3)vmax, podemos obtener:
( )
12 12
⎡ 2 ν⎤ ⎛ 6 Lz v y ⎞
kc Lz 2 ⎢ L2z DAB 3 v y ν Lz
= ⋅ ⋅ ⎥ =⎜ ⋅ ⋅ ⋅ ⎟
DAB π ⎢ DAB
2
Ly ν⎥ ⎜π ν DAB Ly ⎟
⎣ ⎦ ⎝ ⎠ (6.2-22)
( ) (L ) ( ) (L )
12 12 12 12
Sh Lz = 1.382 Re Lz Sc z Ly = 1.382 Pem Lz z Ly
donde Sh Lz es el número de Sherwood promedio. Es importante notar que la distancia

característica con la que se calculan los números de Sherwood y de Reynolds es la misma
Lz. Se ha demostrado que esta "correlación" representa muy bien los datos experimentales,
siempre y cuando el tiempos de contacto gas-líquido Ly v y sean pequeños y que el soluto
no penetre apreciablemente en la película (ver Ejemplo 4.6-1), esto es,
Ly vy tiempo de convección
= ≤ 0.1 (6.2-23)
L2z DAB tiempo de difusión
Podemos apreciar que la Ec.(6.2-21) es muy semejante a la (6.2-15). De hecho, los casos de
difusión en un medio seminfinito y el de difusión en una película descendente para tiempos
R. LOBO/2007
de contacto pequeños constituyen el mismo problema; lo que varía es el punto de vista del
observador. Si pudiéramos viajar durante un tiempo vmax /Ly sobre un paquete de fluido que
desciende a una velocidad vmax en la película descendente, observaríamos que la difusión
ocurre hacia un medio estacionario seminfinito. Si ambos problemas son equivalentes,
conviene apuntar que la imagen física de ambos constituye la base de la llamada teoría de
penetración del coeficiente de transferencia de masa que se verá en el Capítulo 7.
Difusión y convección forzada en flujo laminar en la disolución de un sólido. La

situación física es la misma que la ilustrada en la Figura 4.7-1. La característica más
importante de la situación física es que el soluto experimenta un gradiente de velocidad en
dirección normal a la de la transferencia. El flux del soluto A que se disuelve en la interfase
líquido-sólido y se transfiere a la película líquida laminar está dado por la Ec.(4.7-10)
13
D ⎛ a ⎞
N A,0 = AB ⎜ ⎟ (C sat
)
− 0 = NAi (6.2-24)
()
Γ 43 ⎝ 9DAB y ⎠
A
donde Γ(4/3)=0.8930 y a se define por
ρ gLz 2vmax 3 v y
a= = = (6.2-25)
µ Lz Lz
El flux total de A puede expresarse en términos de la velocidad promedio 〈vz〉 sustituyendo la

Ec.(6.2-25) en la (6.2-24):
13
D ⎛ vy ⎞
N A,0 = AB ⎜ ⎟ (C sat
)
− 0 = N Ai (6.2-26)
()
Γ 3 ⎝ AB z ⎟⎠
4 ⎜ 3D L y
A
Entonces, si comparamos las Ecs.(6.1-1) con la (6.2-26), será evidente que:
( CAi − CAb ) = ( CAsat − 0 ) (6.2-27)
13
D ⎛ vy ⎞
kc = AB ⎜ ⎟ (6.2-28)
()
Γ 43 ⎜⎝ 3DAB Lz y ⎠⎟
En esta ecuación es claro que kc es un coeficiente local, ya que varía con la distancia
vertical y. Podemos escribir esta relación en términos adimensionales, de modo semejante a
como lo hicimos para llegar a la Ec.(6.2-19). Si multiplicamos ambos lados de la
R. LOBO/2007
Ec.(6.2-28) por (y/DAB) y multiplicamos y dividimos el cociente dentro del radical por la
viscosidad cinemática, y rearreglamos
13 13
kc y 1 ⎡ y3 vy ν⎤ 1 ⎛ 1 y vy ν y ⎞
= Sh y = ⎛ ⎞ ⎢ ⋅ ⎥ = ⎛ ⎞⎜ ⋅ ⋅ ⋅ ⎟
DAB Γ ⎜⎝ 43 ⎟⎠ ⎢⎣ 3 DAB yLz ν ⎥⎦ 4 ⎜3
Γ ⎜⎝ 3 ⎟⎠ ⎝ ν DAB Lz ⎟
⎠ (6.2-29)
( ) ( y Lz )1 3 = 0.776 ( Pem y ) ( y Lz )1 3
13 13
Sh y = 0.776 Re y Sc
Podemos calcular también un coeficiente de transferencia de masa promedio sobre la

superficie LxLy usando la expresión del flux promedio dado en la Ec.(4.7-11). Si esta
ecuación la escribimos en términos de la velocidad promedio vz , podemos ver que:
13
2DAB ⎛ vy ⎞
kc = ⎜ ⎟ (6.2-30)
()
Γ 73 ⎜⎝ 3DAB Ly Lz ⎟
⎠
donde Γ(7/3)=1.1906. También podemos escribir esta relación en términos adimensionales,

de modo semejante a como lo hicimos para llegar a la Ec.(6.2-22). Si multiplicamos ambos
lados de la Ec.(6.2-30) por (Lz/DAB) y multiplicamos y dividimos el cociente dentro del
radical por la viscosidad cinemática, podemos llegar a:
13 13
kc Lz 2 ⎛ L3z vy ν⎞ ⎛ v y Lz ν L ⎞
= ⎜ ⋅ ⋅ ⎟ = 1.165 ⎜ ⋅ ⋅ z ⎟
()
DAB Γ 7 ⎜ 3 DAB Ly Lz ν ⎟
3 ⎝ ⎠
⎜ ν
⎝
DAB Ly ⎟
⎠ (6.2-31)
( ) (L ) ( ) (L )
13 13 13 13
Sh Lz = 1.165 Re Lz Sc z Ly = 1.165 Pe Lz z Ly
donde Sh Lz es el número de Sherwood promedio. Éste y el número de Reynolds están

basados en la misma distancia característica Lz. Esta correlación también reproduce bien los
datos experimentales si se cumple el criterio expresado en la Ec.(6.2-23), esto es, si se
emplean tiempos de contacto cortos. Los resultados alcanzados son importantes porque,
aunque idealizados, nos ayudan a comprender los efectos de un campo de velocidad sobre
la transferencia de masa convectiva entre una superficie y un fluido en movimiento.
Es útil comparar los exponentes en las correlaciones resultantes de las Ecs.(6.2-22) y

(6.2-31). En la primera el exponente es 1/2 mientras que en la segunda el exponente es 1/3.
En ambos casos el flujo es laminar, pero en el primero la velocidad a la que está expuesto el
soluto que se transfiere es constante (vmax), con cero gradientes, en tanto que en el segundo
el soluto experimenta una velocidad lineal, con un gradiente de velocidad constante.
R. LOBO/2007
Conviene recapitular dos cuestiones que aparecen en los dos casos de transferencia de masa
por convección forzada en la película laminar descendente, ya que tienen implicaciones
generales. (1) A menudo una de las dimensiones que componen el área de transferencia se
mantiene constante (Lx en los casos de la película laminar descendente), por lo que el
coeficiente local aparece sólo como función de la dimensión variable (y en los casos
mencionados). (2) Aunque hemos calculado los coeficientes de transferencia de masa
promedio a partir de la expresión del flux promedio en el área de transferencia, el mismo
resultado se puede obtener promediando el coeficiente de transferencia de masa local sobre
el área de la película líquida expuesta a la transferencia. En general, el coeficiente de
transferencia de masa promedio se define por
1
kc =
Ai ∫ k dA
Ai
c i (6.2-32)
donde Ai es el área interfacial desde donde ocurre la transferencia. Las mismas

observaciones aplican para los números de Sherwood local y promedio.
Para concluir esta sección detengámonos a reflexionar por un momento en los resultados.
En cada uno de los incisos abordamos un problema típico de difusión y convección en un
medio estacionario o en flujo laminar. En todos los casos encontramos una equivalencia
entre el flux promedio calculado a partir de los modelos de la ley de Fick y del coeficiente
de transferencia de masa. Esta equivalencia no debe sorprendernos, ya que la realidad que
tratan de describir ambos modelos es la misma. Sin embargo, reiteremos que esta
equivalencia se refiere únicamente al cálculo del flux promedio, pues con el modelo del
coeficiente de transferencia de masa no es posible conocer el perfil de concentraciones en la
zona de transferencia, sino solamente la diferencia de concentración entre los dos planos
limítrofes de dicha zona.
Por otra parte, es importante destacar que en las situaciones que acabamos de analizar, no
ha estado presente el mecanismo de convección forzada en flujo turbulento, que es común
en muchos problemas de transferencia de masa en interfases. En la Sec. 6.3 explicaremos
cómo opera el mecanismo de transferencia de masa por convección forzada en flujo
turbulento y veremos algunos ejemplos sobre cómo calcular el valor del coeficiente de
transferencia de masa en estas situaciones.
Ejemplo 6.2-1. Absorción de CO2 en una película laminar descendente
Se tiene una superficie vertical de 100 cm de largo y 62.8 cm de ancho, la cual se irriga
uniformemente con un flujo másico de agua de 78 g/s. El agua (A) que desciende forma
una película laminar de 3.24×10−2 cm de espesor. Paralelo a esta película líquida fluye CO2
puro (B), el cual se absorbe en ella. Se cuenta con los siguientes datos a 25 °C: la densidad
R. LOBO/2007
de la solución CO2-agua que se forma es 1 g/cm3; su viscosidad es 0.89 cp; el coeficiente de

difusión del CO2 en agua es 2×10−5 cm2/s; y la solubilidad del CO2 en agua es de 3.37×10−4
mol/cm3.
Calcular:
(a) el coeficiente de transferencia de masa local;
(b) el coeficiente de transferencia de masa promedio;
(c) la concentración promedio; y
(d) el flujo molar total de CO2 absorbido, a los 50 y 100 cm de la longitud vertical de la
película.
Solución. (a) En primer lugar, debemos verificar si las condiciones bajo las cuales la teoría
desarrollada en la Sec. 4.6 se cumplen. Los criterios de aplicación para esta teoría están
dados por las Ecs.(4.6-28) y (4.6-41):
vmax Lz
DAB
>> 1 ; (
0 (Y ) ≤ 0 103 cm ) (6.2-33)
donde Y es la longitud vertical de la película. Evidentemente, el segundo criterio se cumple,

ya que la longitud vertical de la película líquida en nuestro caso es de 100 cm. Debemos
ahora verificar si el primer criterio se cumple. Para ello, calculemos primero vmax, sabiendo
que en un flujo laminar, como es el de nuestro problema, la velocidad promedio se
relaciona con la velocidad máxima mediante vmax=(3/2)〈vy〉. La velocidad promedio 〈vy〉 se
relaciona con el flujo másico wy mediante
wy = ρ v y AS (6.2-34)
donde AS es el área transversal al flujo. Entonces,
wy wy 78 g/s
vy = = = = 38.3 cm/s
ρ AS ρ Lx Lz (1 g/cm )( 62.8 ⋅ 3.24 × 10
3 −2
cm
2
) (6.2-35)
vmax =
3
2
vy =
3
2
( 38.3 cm/s ) = 57.5 cm/s
Veamos ahora si se cumple el primer criterio:
vmax Lz ( 57.5 cm/s ) ( 3.24 × 10 cm )

−2
= −5 2
4
= 9.32 × 10 >> 1 (6.2-36)
DAB 2 × 10 cm /s
Es claro entonces que sí podemos usar la teoría desarrollada en la Sec. 4.6 y sus
consecuencias en el cálculo del coeficiente de transferencia de masa, el flux y el flujo
molares expuestos en esta sección.
R. LOBO/2007
El coeficiente de transferencia de masa local está dado por la Ec.(6.2-18):
DAB vmax
kc =
πy
Si evaluamos esta expresión en y=50 y 100 cm de longitud vertical de la película

obtenemos el valor del coeficiente de transferencia de masa local a estas distancias:
kc =
DAB vmax
=
( 2 × 10 −5 2
cm /s ) ( 57.5 cm/s ) = 2.71× 10 −3
cm/s (6.2-37)
π ( 50 cm )
50 cm
π y 50 cm
kc =
DAB vmax
=
( 2 × 10 −5
cm /s
2
) ( 57.5 cm/s ) = 1.91× 10 −3
cm/s (6.2-38)
π (100 cm )
100 cm
π y 100 cm
(b) El coeficiente de transferencia de masa promedio para el problema que estamos tratando
está dado por la Ec.(6.2-21):
DAB vmax
kc = 2
π Ly
Si evaluamos esta expresión para una película de longitud Ly=50 y 100 cm de longitud
vertical de la película obtenemos:
( 2 × 10 −5 2
cm /s ) ( 57.5 cm/s ) = 5.42 × 10
( k c )0
50 cm −3
= 2⋅ cm/s (6.2-39)
π ( 50 cm )
( 2 × 10 −5 2
cm /s ) ( 57.5 cm/s ) = 3.82 × 10
( kc )0
100 cm −3
= 2⋅ cm/s (6.2-40)
π (100 cm )
En estas ecuaciones se indica el rango de distancias sobre las cuales se está obteniendo el
coeficiente de transferencia de masa promedio. Debemos notar también la diferencia de
magnitud entre los coeficientes de transferencia de masa locales y los promedio. La
magnitud de estos últimos es dos veces la de los coeficientes locales a la misma distancia.
(c) La concentración promedio se puede calcular de la Ec.(4.6-49)
R. LOBO/2007
CAsat 4DAB Ly
CA b y = Ly =
Lz π vmax
Si sustituimos los valores correspondientes en esta ecuación:
(3.37 × 10 mol/cm ) −4 3
(
4 ⋅ 2 × 10
−5 2
cm / s ) (50 cm )
CA b L =50 cm =
y
(3.24 × 10 cm ) −2
π ( 57.5 cm/s )
−5 3
(6.2-41)
= 4.89 × 10 mol/cm
CA b =
(3.37 × 10 mol/cm )−4 3
(
4 ⋅ 2 × 10
−5 2
)
cm / s (100 cm )
L y =100 cm
(3.24 × 10 cm ) −2
π ( 57.5 cm/s ) (6.2-42)
−5 3
= 6.92 × 10 mol/cm
De estos resultados podemos percatarnos que la concentración ha aumentado más en los

primeros 50 cm que en los siguientes 50 cm. Ello refleja que la absorción es más rápida en
los primeros 50 cm que en los segundos, dada la diferente magnitud de los coeficientes de
transferencia de masa promedio en cada zona de la película. Esta misma tendencia debe
observarse en el cálculo del flujo molar. Veamos.
(d) El flujo molar absorbido lo podemos calcular de la Ec.(4.6-15):
WA,0 = Lx Ly
4DAB vmax
π Ly
CAsat − 0 ( )
La sustitución de los valores correspondientes en esta ecuación produce:
(
4 2 × 10
−5 2
)
cm / s ( 57.5 cm/s )
(WA,0 )0 = ( 62.8 × 50 cm ) (3.37 × 10 )=
50 cm 2 −4 3
mol/cm
π ( 50 cm )
−3
(6.2-43)
= 5.73 × 10 mol/s
(
4 2 × 10
−5 2
)
cm / s ( 57.5 cm/s )
(WA,0 )0 = ( 62.8 × 100 cm ) (3.37 × 10 )=
100 cm 2 −4 3
mol/cm
π (100 cm )
−3
(6.2-44)
= 8.10 × 10 mol/s
Las magnitudes de los flujos molares guardan la misma proporción que las magnitudes de
las concentraciones promedio, que a su vez guardan la misma proporción inversa a las
magnitudes de los coeficientes de transferencia de masa locales y promedios, esto es,
R. LOBO/2007
(WA,0 )0 CA b L =50 cm ( kc )
50 cm 100 cm
y kc 100 cm
= = 0
= = 0.7 (6.2-45)
(WA,0 )0 CA b L =100 cm ( kc )
100 cm 50 cm
kc 50 cm
y 0
6.3 Transferencia de masa por convección forzada en flujo turbulento

Con el propósito de ubicarnos en nuestro estudio de la transferencia de masa interfacial, es
conveniente iniciar esta sección con una síntesis de lo que hemos avanzado hasta el
momento. En la Sec. 6.1 introducimos una definición empírica del coeficiente de
transferencia de masa, como una manera "sencilla" de calcular el flux de masa en las
interfases. Ahí señalamos que su definición no es única, ya que existen diferentes tipos de
coeficiente, dependiendo de la forma de expresar la fuerza motriz de la transferencia
(∆CA). Indicamos también los requisitos de definir claramente la diferencia de
concentración y el área a través de la cual ocurre la transferencia. Dejamos establecido
también que debemos preferir el uso de coeficientes de transferencia de masa y diferencias
de concentración locales porque representan mejor el proceso de transferencia de masa.
En la Sec. 6.2 recapitulamos una serie de situaciones de transferencia de masa en las que
encontramos una equivalencia para el cálculo del flux entre el uso de la ley de Fick y el del
coeficiente de transferencia de masa. Lo característico de estos casos fue que el medio
donde ocurre la transferencia de masa interfacial era estacionario (v=0), o donde el flujo era
laminar1. En dichas situaciones mostramos que era posible encontrar una expresión teórica
para calcular el coeficiente de transferencia de masa; de hecho, cuando el flujo es laminar
siempre es posible calcular teóricamente el coeficiente de transferencia de masa. Sin
embargo, dichas predicciones teóricas pueden desviarse de los valores experimentales y
deben ser corregidas para una mayor precisión.
Ahora es necesario que avancemos nuestro estudio de la transferencia de masa interfacial

cuando el flujo es turbulento. ¿Cómo ocurre la transferencia de masa en flujo turbulento?
¿De qué depende el coeficiente de transferencia de masa? ¿Cómo obtener valores del
coeficiente de transferencia de masa para situaciones prácticas? ¿Cómo se expresa este
coeficiente en términos del número de Sherwood cuando el flujo es turbulento? A
responder estas preguntas dedicaremos el resto de este capítulo.
En esta sección haremos una descripción cualitativa de la transferencia de masa en flujo

turbulento en una interfase sólido-fluido desde el punto de vista de su mecanismo2. Hay que
señalar, sin embargo, que en muchas de las interfases de interés práctico, principalmente las
fluido-fluido los fenómenos son aún más complejas. El marco conceptual que usaremos
proviene de la teoría de la capa límite, pero en estos momentos nos abstendremos de su
descripción matemática para concentrarnos en los conceptos. Como ganancia adicional de
1 En flujo laminar, el mecanismo de transferencia de masa perpendicular al flujo es la difusión.

2 El estudiante interesado puede consultar mayores detalles en la sección correspondiente en el Capítulo 7.
R. LOBO/2007
esta descripción nos será posible racionalizar la definición del coeficiente de transferencia
de masa en flujo turbulento, esto es, darle un sentido “teórico”.
Consideremos el flujo sobre la placa plana horizontal hecha de un compuesto A poco

soluble que se muestra en la Figura 6.3-1. La solución que el soluto A forma con el fluido
es diluida. El fluido en movimiento pasa sobre la superficie sólida en condiciones de flujo
turbulento, donde la velocidad promedio en la corriente libre es vb,z. La velocidad de
pequeñas partículas o "paquetes" de fluido no puede describirse de una manera clara como
en flujo laminar. En flujo turbulento, no existen las líneas de corriente ordenadas que
caracterizan al flujo laminar, sino que pequeñas porciones o paquetes de fluido se mueven
rápidamente casi al azar en las tres direcciones, formando pequeños remolinos (eddies, en
inglés) Aunque macroscópicamente un flujo turbulento puede parecer estable y
unidireccional, a nivel microscópico es inherentemente inestable y tridimensional, ya que
existen pequeñas fluctuaciones en el valor de la velocidad promedio1 (ver Sec. 7.6).
Cuando el soluto se disuelve desde la superficie sólida y se transfiere al fluido, existe una
concentración comparativamente alta de dicho soluto en el fluido en inmediato contacto
con la superficie, misma que va disminuyendo conforme aumenta la distancia desde la
superficie. Si pudiéramos analizar pequeñas porciones de fluido, resultaría que muestras
adyacentes no tendrían siempre la misma concentración. Ello se debe a que los pequeños
remolinos de fluido están acarreando continua y rápidamente el soluto hacia el seno del
fluido. Esta acción de los remolinos da como resultado una transferencia de masa rápida en
comparación con la que ocurriría por difusión molecular.
corriente libre vb. z

CAb vb. z
remolinos
o "eddies"
CAb vb. z
vb. z
CAb
zona de capa límite
CA ( y, z ) transición turbulenta
y N Ai N Ai N Ai
CAi v ( y, z )
↑↑↑ ↑↑↑ ↑↑↑
z
capa límite laminar subcapa laminar
1 En flujo turbulento todas las propiedades fluctúan alrededor de un valor promedio. De hecho, en flujo
turbulento vb,z es un doble promedio. Es un promedio de la velocidad sobre el área perpendicular al flujo y es
un promedio en el tiempo (ver Sec. 7.6). De hecho deberíamos escribir vb, z y no vb, z . Sin embargo, por
simplicidad en la notación mantendremos esta última.
R. LOBO/2007
Figura 6.3-1. Transferencia de masa por convección forzada en flujo turbulento en la capa límite.
Debido a que el sólido frena al fluido que está en inmediato contacto con él, en la vecindad
de la interfase se forma una película donde el fluido fluye en régimen laminar. Esta zona
de denomina la capa límite laminar y en ella el gradiente de velocidad es grande. En la
capa límite laminar no existen remolinos, por lo que el proceso de transferencia de masa en
dirección normal al flujo se da únicamente por difusión molecular. Como el transporte de
masa por difusión es lento, el gradiente de concentración en la zona laminar es grande, ya
que la concentración cambia de CAi hasta CAb sobre una distancia muy pequeña. Cabe
señalar que en la Figura 6.3-1 se muestra el mismo grosor para la capa límite de velocidad
(o hidrodinámica) y para la capa límite de concentración, aunque por lo general la capa
límite de concentración es más delgada que la hidrodinámica; asimismo, las diferentes
zonas aparecen claramente delimitadas, aunque en realidad no hay una frontera precisa
entre ellas.
Adyacente a la capa límite laminar, existe una zona que se denomina de transición, donde
comienzan a aparecer los remolinos ocasionados por el flujo turbulento. A la acción de
acarreo de masa del soluto A por los remolinos se le llama difusión turbulenta1. La
transferencia de masa de A se incrementa y se da por la suma de la difusión molecular y la
difusión turbulenta. En la región de transición no existen cambios abruptos, sino graduales,
de transporte molecular de masa a transporte turbulento. Los gradientes de concentración
no son tan marcados como en la película laminar.
Contigua a la región de transición, se encuentra la región turbulenta, donde la mayor parte

del transporte de masa se debe al acarreo del soluto A por los remolinos, esto es, a la
difusión turbulenta. Se ha demostrado que la intensidad de la turbulencia es más alta en las
regiones cercanas a la interfase que en la corriente libre del fluido y que la acción de los
remolinos puede llega hasta ella. En la región turbulenta, la difusión molecular continúa
presente, pero contribuye poco al transporte total de masa. Los cambios de concentración en
esta zona son pequeños, ya que los remolinos tienden a uniformar la concentración y a
disminuir los gradientes que pudieran existir.
Esta esquemática descripción del proceso de transferencia de masa en una interfase bien
definida y con un flujo relativamente bien caracterizado da una idea sobre la complejidad
del estudio de la transferencia de masa en las interfases. Se presenta entonces un nuevo
problema: ¿cómo cuantificar la difusión turbulenta? Debido a que la teoría de transporte
turbulento no se encuentra aún suficientemente desarrollada2, se ha intentado cuantificar la
difusión turbulenta por expresiones análogas a las de difusión molecular, esto es
1 La difusión turbulenta es un mecanismo de transporte de masa de origen convectivo y diferente de la

difusión molecular.
2 Diversos expertos consideran al flujo turbulento como el principal problema no resuelto de la física clásica.
R. LOBO/2007
(t ) () t dCA
J A, y = − DAB (6.3-1)
dy
En esta ecuación, J A,(t )y representa el flux por difusión turbulenta en la dirección y normal al
flujo. DAB( t ) es la difusividad turbulenta o difusividad de remolino (eddy diffusivity en
inglés), y tiene las mismas unidades que el coeficiente de difusión molecular; sin embargo,
a diferencia de éste, DAB ( t ) no es una propiedad de la mezcla, sino que es una variable que
depende principalmente de las características del flujo y de la geometría del sistema. CA es

la concentración promedio en el tiempo1 cuyo significado preciso debe consultar el
estudiante en la Sec. 7.6.
El valor de DAB( t ) es cero en la superficie sólida y se incrementa conforme aumenta la
distancia desde ésta. Si suponemos que los fluxes por difusión molecular y turbulenta son
aditivos, y si usamos un valor promedio de la difusividad turbulenta, DAB ( t ) , es posible
integrar una expresión para el flux difusivo total desde la interfase
()
( )( )
t
D + DAB
( total )
J Ai = AB
Ly
( )
CA i − CA b [ = ] flux difusivo + flux difusivo
molecular turbulento
(6.3-2)
donde Ly es una distancia indefinida sobre la cual ocurre el proceso de transferencia de

masa. En adelante, omitiremos la barra sobre la difusividad turbulenta y entenderemos que
siempre usaremos su valor promedio. Es necesario recalcar que la Ec.(6.3-2) es una
racionalización que hacemos para tomar en cuenta el transporte de masa debido al flujo
turbulento, transporte que frecuentemente se describe cuantitativamente con una ecuación
semejante a la de la difusión molecular.
Para darle un sentido "teórico" al coeficiente de transferencia de masa en flujo turbulento

podemos comparar esta expresión del flux difusivo total con una en la cual el flux de masa
interfacial se exprese en términos del coeficiente de transferencia de masa2:
() t
DAB + DAB
( total )
J Ai =
Ly
( ) (
CA i − CA b = N A i = kc CA i − CA b ) (6.3-3)
1 Nuevamente hay que recordar que en flujo turbulento todas las propiedades fluctúan alrededor de un valor
promedio. CA que es el valor promedio de la concentración que se usa en esta sección. De hecho, en flujo
turbulento deberíamos escribir como CAb la concentración promedio CAb . Sin embargo, por simplicidad en
la notación mantendremos esta última.
2 En vista de que estamos considerando que la velocidad del fluido en la superficie del sólido es cero,
podemos hacer
J A( i ) = N A i
total
R. LOBO/2007
de donde es evidente que
() t
DAB + DAB
kc = (6.3-4)
Ly
Una vez derivada la Ec.(6.3-4), pudiéramos pensar que hemos obtenido algo concreto sobre
(t )
el coeficiente de transferencia de masa en flujo turbulento. Esto no es así, porque Ly y DAB
son desconocidos. Sin embargo, sí hemos avanzado en formarnos una idea de lo que sucede
cuando la difusión y la convección forzada turbulenta actúan simultáneamente como
mecanismos de transporte de masa en la vecindad de las interfases. También hemos podido
percatarnos de las dificultades para obtener expresiones rigurosas para describir dichos
procesos. De hecho, prácticamente todas las ecuaciones existentes para calcular el
coeficiente de transferencia de masa en situaciones en las que el flujo es turbulento se han
obtenido con base en correlaciones empíricas de mediciones experimentales. Esto marca
una diferencia importante con respecto a la transferencia de masa en flujo laminar.
Recordemos además que la discusión precedente se basó en que la solución hacia la cual
ocurría la transferencia de masa era diluida.
Cuando la solución que se forma en la vecindad de la interfase es concentrada (CAi grande)

y/o existen altas tasas de transferencia de masa (NAi+NBi grandes), debemos tomar en
cuenta el término convectivo del flux total de masa. La definición del flux total en términos
del coeficiente de transferencia de masa que incluye este efecto convectivo es1:
i (
N A i = J A(total) + CA i C )( N Ai + N Bi )= kc ( CAi − CAb ) + ( CAi C )( N Ai + N Bi ) (6.3-5)
de donde podemos despejar el coeficiente de transferencia de masa
kc =
(
N A i − CA i C )( NAi + NBi ) (6.3-6)
CA i − CA b
Ciertamente, esta expresión es más compleja que la que obtendríamos de la Ec.(6.1-3). La

diferencia está en el término ( CA i C )( N A i + N Bi ) . Si recordamos el análisis realizado en la Sec.
4.1, notaremos que este término se refiere a la convección que es importante cuando las
soluciones son concentradas y/o la rapidez de la transferencia de masa son altas. Cuando las
soluciones son diluidas puede suponerse razonablemente que la velocidad en la interfase es
cero. Sin embargo, cuando las soluciones son concentradas y/o cuando existen altas tasas de
transferencia de masa la velocidad en la interfase ya no puede suponerse igual a cero. Cuando
1 Bird et al. op. cit.
R. LOBO/2007
debemos considerar el término convectivo, es común definir el coeficiente de transferencia de

masa como1:
kc
NAi =
ϕN
( CAi − CAb ) (6.3-7)
donde
⎛ NAi CA b ⎞ ⎛ N A i CA i ⎞
⎜⎜ − ⎟⎟ − ⎜⎜ − ⎟⎟
N A i + N Bi C ⎠ ⎝ N A i + N Bi C
ϕN = ⎝ ⎠ (6.3-8)
⎛ NAi CA b ⎞
⎜ − ⎟
NAi ⎜ N A i + N Bi C ⎟
ln
N A i + N Bi ⎜ NAi CA i ⎟
⎜⎜ − ⎟
⎝ N A i + N Bi C ⎟⎠
es una especie de factor de corrección que debe aplicarse al coeficiente de transferencia de

masa en solución diluida. Cuando la transferencia se da en una situación de contradifusión
equimolar ϕN=1. En otros casos, el factor de corrección puede calcularse de diversas teorías; de
ellas, la más sencilla es la teoría de la película, que será la que utilicemos para calcular el factor
de corrección en el Capítulo 7.
Así, nuestro enfoque para soluciones los problemas donde existan altas concentraciones y/o
altas tasas de transferencia de masa consistirá en seguir utilizando la definición
(
N A i = k c CA i − CA b )
y, cuando el término convectivo ( CA i C )( N A i + N Bi ) no pueda ser ignorado, utilizaremos
los coeficientes de transferencia de masa obtenidos para solución diluida, haciéndoles una
corrección. El examen cuidadoso de las Ecs.(6.3-7) y (6.3-8) nos lleva a obtener una
conclusión muy importante: el coeficiente de transferencia de masa es función de la
rapidez de la transferencia de masa y de la magnitud de la concentraciones involucradas.
El origen de la funcionalidad de kc con la rapidez de la transferencia de masa reside en la

alteración de los perfiles de velocidad y concentración en la vecindad de la interfase causados
por el flujo del compuesto o compuestos que se transfieren a través de ella. En la situación
expuesta en la Figura 6.3-1 consideramos que el flux NAi que sale de la interfase y que es
perpendicular a la dirección principal del flujo, no afecta al perfil de velocidades ni al perfil
de concentraciones. Sin embargo, esto no siempre ocurre así. Por fortuna, este efecto es
pequeño para muchos casos de interés práctico, de modo que podemos utilizar la Ec.(6.1-1)
1 Bedingfield, CH y Drew, TB. (1950). Ind. Eng. Chem. 42, 1164.
R. LOBO/2007
como punto de partida. Una vez entendido el problema y si el caso lo amerita, podemos
considerar los efectos señalados y hacer las correcciones necesarias.
Ejemplo 6.3-1. Más sobre la definición del coeficiente de transferencia de masa
A menudo aparecen en la literatura especializada expresiones para calcular el flux o el flujo

de masa en algún equipo en términos de distintos coeficientes de transferencia de masa y
diferencias de concentraciones. Por ejemplo, pueden encontrarse expresiones como
⎡ ∆CA,ent − ∆CA,sal ⎤
N A = kln ⎢ ⎥ (6.3-9)
⎢⎣ ln ( ∆CA,ent ∆CA,sal ) ⎥⎦
donde los subíndices ent y sal indican las condiciones a la salida y a la entrada del equipo,
respectivamente. Por otra parte, en este texto nosotros hemos adoptado un coeficiente de
transferencia de masa local y una diferencia de concentración también local,
N A i = kc ( CA i − CA b ) porque representan mejor el proceso real de transferencia de masa que
está ocurriendo. En la Sec. 6.1 señalamos que un coeficiente de transferencia de masa no
queda definido en tanto no se defina un área de transferencia y una diferencia de
concentración. El problema consiste en usar este señalamiento para encontrar la relación
entre los coeficientes de transferencia de masa kc y kln, así como el significado de la
diferencia de concentración media logarítmica arriba indicada.
Solución. Consideremos el flujo de un fluido (que puede ser turbulento o laminar) dentro
de un tubo cuya pared es ligeramente soluble en el fluido, como se muestra en la Figura
6.3-2. La velocidad promedio del fluido es vb,,z y la concentración promedio del soluto a la
entrada del tubo es CAb(0)=CA0.
pared poco
∆z soluble
CAi N Ai
CA0 CAb ( Lz )
N A,z z CAb ( z ) N A,z z +∆z
vb , z z vb, z
R N Ai
Lz
Figura 6.3-2. Disolución desde la pared de un tubo. El balance de masa se realiza en términos de la
concentración promedio CAb.
R. LOBO/2007
En estado estacionario, el balance de masa sobre el soluto A en el elemento diferencial de

volumen 2πR∆z es:
⎛ flujo molar de A que entra en z − ⎞ ⎛ flujo molar de A transferido ⎞

⎜ ⎟+⎜ ⎟=0
⎝ flujo molar de A que sale en z + ∆z ⎠ ⎝ desde la superficie 2π R∆z ⎠
En términos matemáticos:
(
π R 2 N A,z z − N A,z z +∆z )+ N A i ∆Ai =0 (6.3-10)
donde NAi es el flux que se transfiere desde la superficie del sólido y ∆Ai es el área
interfacial de transferencia contenida en el elemento diferencial de volumen. Definamos las
cantidades que componen este balance:
( total )
N A, z = J A,z + CA b vb, z ≅ CA b vb, z ; (
N A i = k c CA i − CA b )
(6.3-11)
∆Ai = 2π R∆z
Debemos notar que hemos supuesto que el flux difusivo total en dirección z es
insignificante en relación con el convectivo; estamos usando vb,,z para representar la
velocidad promedio del fluido, la cual consideraremos constante. Si sustituimos estas
relaciones en la Ec.(6.3-10), dividimos por 2πR∆z, tomamos el límite cuando ∆z→0 y
rearreglamos, llegamos a:
d CA b 2
−
dz
+
Rvb, z
(
k c CA i − CA b = 0 ) (6.3-12)
sujeta a la condición de frontera:
z=0; CA b = CA b ( 0 ) = CA0 (6.3-13)
Examinemos por un momento la Ec.(6.3-12). Arriba dijimos que se puede definir un

coeficiente de transferencia de masa sólo después de definir un área de transferencia y una
diferencia de concentración. En nuestro caso, para realizar el balance de masa hemos
escogido una diferencia de concentración local ∆CA=CAi−CAb y un área interfacial
∆Ai=2πR∆z, lo cual significa que kc es el valor promedio en el área ∆Ai; al hacer ∆Ai→0, kc
se torna en el coeficiente de transferencia de masa local. Esto lo podemos apreciar también
si escribimos la Ec.(6.3-12) como
Rvb, z 1 dCA b
kc = = kc ( z ) (6.3-14)
2 CA i − CA b d z
R. LOBO/2007
donde es evidente que kc varía1 con z.
Retomemos ahora el desarrollo de nuestro ejemplo. Definamos por conveniencia:
∆C A = C A i − C A b (6.3-15)
Podemos sustituir esta definición en la Ec.(6.3-12), para luego separar variables
d ( ∆CA ) 2 Lz 1
=− ⋅ kc dz (6.3-16)
∆C A Rvb, z Lz
Por razones que serán evidentes abajo, hemos multiplicado y dividido el lado derecho de
esta ecuación por Lz. Si integramos esta ecuación entre la entrada (z=0) y la salida (z=Lz)
del tubo no es difícil llegar a:
⎛ ∆CA,sal ⎞ 2 Lz 1 Lz 2 Lz
ln ⎜
⎜ ∆C
⎝ A,ent
⎟⎟ = −
⎠ Rvb , z
⋅
Lz 0∫

kc dz = −
Rvb , z
⋅ kc (6.3-17)
kc
En el término medio de esta ecuación aparece la definición del coeficiente de transferencia

de masa promedio sobre toda la superficie interna del tubo 2πRLz. De esta ecuación
podemos obtener
Rvb, z ⎛ ∆C ⎞ v ⎛D ⎞ ⎛ ∆CA,ent ⎞
kc = ln ⎜ A,ent ⎟⎟ = b. z ⎜ ⎟ ln ⎜⎜ ⎟⎟ (6.3-18)
⎜ ∆C
2 Lz ⎝ A,sal ⎠ 4 ⎝ Lz ⎠ ⎝ ∆CA,sal ⎠
donde
∆CA,ent CA i − C A b ( 0 )
= (6.3-19)
∆CA,sal CA i − CA b ( Lz )
y D es el diámetro del tubo. Es muy importante notar que el coeficiente promedio es

función de la concentración, la velocidad promedio y de la razón de diámetro a longitud del
tubo (D/Lz).
1 Sin embargo, debemos señalar que esta variación de kc con la distancia se da en la zona de entrada del tubo
y hasta una cierta distancia, en la cual ni el flujo ni el perfil de concentraciones están completamente
desarrollados. Después de una longitud lo suficientemente larga kc adquiere un valor constante (y sigue siendo
un coeficiente local). Cuando el flujo es laminar, en la zona totalmente desarrollada kc=3.66(DAB/D), donde D
es el diámetro del tubo. Más adelante se discutirán más los valores de kc en flujo turbulento.
R. LOBO/2007
Naturalmente, podemos integrar la Ec.(6.3-16) desde la entrada del tubo hasta alguna
posición axial z cualquiera. En este caso obtenemos después de rearreglar:
CA i − CA b ( z ) ⎛ 4k z ⎞
= exp ⎜ − c ⎟ (6.3-20)
CA i − CA b ( 0 ) ⎜ v D⎟
⎝ b, z ⎠
donde ahora kc es el coeficiente de transferencia de masa promedio desde la entrada del

tubo hasta la distancia axial z. Notemos que aunque CAb siempre crece a lo largo del tubo
(ver Figura 6.3-3), el cociente de la concentración en el lado izquierdo de esta ecuación
vale la unidad en z=0 y decrece exponencialmente a lo largo del tubo.
CAb CAi
∆CA,sal
( C A i − CA b )
∆CA,ent
CAb ( z )
0 z Lz
Figura 6.3-3. Perfil de concentraciones en la transferencia de un soluto desde la pared interna de un tubo.
Ahora bien, si hacemos un balance de masa macroscópico (entre la entrada y la salida del
tubo) y lo relacionamos con la rapidez de transferencia de masa en todo el tubo podemos
escribir:
⎛ flujo molar de A que sale del tubo ⎞ ⎛ flujo molar promedio transferido ⎞
⎜ −flujo molar de A que entra al tubo ⎟ = ⎜ desde toda la superficie del tubo ⎟
⎝ ⎠ ⎝ ⎠
(π R v ) ⎡⎣C
2
b, z Ab ( Lz ) − CA b ( 0 ) ⎤⎦ = WA i = kc ( 2π RLz )( ∆CA )m (6.3-21)
En esta ecuación (πR2vb,z) es el flujo volumétrico a lo largo del tubo, WAi es el flujo molar
promedio transferido desde la toda el área de transferencia de masa Ai=2πRLz. Lo que no
está definido en esta ecuación es la diferencia de concentración (∆CA)m adecuada para usar
con el coeficiente de transferencia de masa promedio kc definido en la Ec.(6.3-17).
Entonces, si sustituimos kc de la Ec.(6.3-18) en la (6.3-21) y rearreglamos, obtendremos la
definición de (∆CA)m:
R. LOBO/2007
∆CA,ent ∆CA,sal

CA b ( Lz ) − CA b ( 0 ) [CA i − CA b ( 0 )] − [CA i − CA b ( Lz )]
( ∆C A ) m = = (6.3-22)
⎛ ∆CA,ent ⎞ ⎛ ∆CA,ent ⎞
ln ⎜ ⎟⎟ ln ⎜ ⎟⎟
⎜ ∆CA,sal ⎜ ∆CA,sal
⎝ ⎠ ⎝ ⎠
Para rearreglar el numerador de esta ecuación hemos sumado y restado la concentración

interfacial CAi. Es entonces claro que la diferencia de concentración (∆CA)m es la clásica
definición de la diferencia de concentración media logarítmica:
∆CA,ent − ∆CA,sal [CA i − CA b ( 0 )] − [CA i − CA b ( Lz )]

( ∆CA )m = ( ∆CA )ln = = (6.3-23)
⎛ ∆CA,ent ⎞ ⎡ CA i − C A b ( 0 ) ⎤
ln ⎜ ⎟⎟ ln ⎢ ⎥
⎜
⎝ ∆CA,sal ⎠ ⎢⎣ CA i − CA b ( Lz ) ⎦⎥
que es la misma que aparece en la Ec.(6.3-9). Esta diferencia de concentración también

puede escribirse como:
∆CA,sal − ∆CA,ent [CA i − CA b ( Lz )] − [CA i − CA b ( 0 )]

( ∆CA )ln = = (6.3-24)
⎛ ∆CA,sal ⎞ ⎡ CA i − CA b ( Lz ) ⎤
ln ⎜ ⎟⎟ ln ⎢ ⎥
⎣⎢ CA i − CA b ( 0 ) ⎦⎥
⎜
⎝ ∆CA,ent ⎠
Resumamos lo que hemos encontrado en nuestro análisis. En primer lugar, el denominado

coeficiente de transferencia de masa logarítmico kln no es otra cosa que el coeficiente de
transferencia de masa promedio kc , mismo que se obtiene de promediar el coeficiente de
transferencia de masa local a lo largo del área interfacial del equipo. La diferencia de
concentración media logarítmica surge de hacer un balance de masa macroscópico (sobre
todo el equipo), e igualarlo al flujo promedio transferido en todo el tubo, y resulta ser la
diferencia de concentración adecuada a usar con el coeficiente promedio kc .
¿Por qué no usar directamente el coeficiente promedio y la diferencia de concentración

media logarítmica? ¿Por qué pasar por el balance de masa diferencial de masa usando un
coeficiente y una diferencia de concentración local? En mi experiencia como profesor he
encontrado que el estudiante, especialmente el principiante, debe recordar constantemente
que la transferencia de masa en la realidad es un proceso local. No puedo sino recomendar
al estudiante que siempre realice los balances diferenciales y use coeficientes locales. Es en
el proceso de integración de dichos balances que surgen los coeficientes promedio y las
diferencias de concentración adecuadas a ellos. Esta recomendación quizá implique un
poco más de trabajo, pero creo que las probabilidades de cometer errores conceptuales (y
de cálculo) se reducirán apreciablemente.
R. LOBO/2007
Ejemplo 6.3-2. Determinación experimental del coeficiente de transferencia de

masa en un lecho empacado
Se tiene un pequeño lecho cilíndrico de 60 cm de longitud, empacado de esferas de naftaleno

(A) de 0.2 cm de diámetro. El área específica del lecho de las esferas empacadas es
aproximadamente de 25 cm2 por cada cm3 de volumen del lecho y se representa por a. A este
lecho se alimenta aire (B) puro y seco a 1 atm y 75 ºC, a una velocidad superficial1 de 2 m/s.
La presión parcial del naftaleno al final del tubo es de 2 mm de Hg. La presión de vapor del
naftaleno a la temperatura del aire es 5 mm de Hg. Puede suponerse que durante el tiempo
que dura el experimento las esferas de naftaleno no cambian su tamaño.
Calcula el coeficiente de transferencia de masa para la evaporación del naftaleno.
aire seco T = 75 °C
CA b = 0 P = 1 atm
z N A,z
z
z + ∆z
60 cm
A N Ai
sólo el área no sombreada
está disponible al flujo
aire + naftaleno
Figura 6.3-4. Representación del sistema de la transferencia de masa en el lecho empacado del Ejemplo 6.3-2.
Solución. El sistema se representa en la Figura 6.3-4. En primer lugar, y aunque parezca

obvio, hay que decir que el gas es la fase de interés. Éste fluye a lo largo del lecho
empacado a una velocidad que debe definirse claramente, ya que de esta definición
dependerá qué área transversal al flujo sea la correcta. Si vamos a usar la velocidad
superficial para cuantificar el flujo del gas, debemos ser consistentes y usar como área
1 La velocidad superficial es una forma común de expresar la velocidad de un fluido en lechos empacados. La
velocidad superficial es la que tendría el fluido suponiendo que el tubo está vacío, esto es, suponiendo que
toda el área seccional del tubo está disponible al flujo. Otra velocidad común es la intersticial, que se define
como la velocidad superficial dividida por la fracción volumétrica de huecos en el lecho (fracción vacía); en
esta definición sólo una fracción del área seccional del tubo está realmente disponible al flujo del fluido.
Como se verá en este ejemplo, hay que definir claramente con que velocidad se trabaja.
R. LOBO/2007
transversal al flujo el área seccional del tubo vacío; si usáramos la velocidad intersticial,
deberíamos usar como área transversal al flujo el área transversal al flujo de los huecos que
deja libre el sólido (ver el dibujo central de la Figura 6.4-3). En nuestro caso, el dato que
tenemos es de la velocidad superficial y, en consecuencia, consideraremos que el área
transversal al flujo es el área seccional del tubo vacío.
Observamos que el naftaleno se transfiere desde la superficie de las esferas hacia el seno
del aire. El área total de transferencia de masa será aproximadamente el área de todas las
esferas que componen el lecho empacado. En nuestro caso, este dato se expresa
normalizado por el volumen del lecho como el área específica a=25 cm2/cm3.
La concentración del naftaleno en el aire irá aumentando a lo largo del lecho en la dirección
del flujo. La diferencia de concentración que rige el proceso de transferencia de masa varía
entre los valores de la concentración en la superficie de las esferas (que es constante e igual
a la de saturación) y la que prevalece en el seno del aire; en otras palabras, usaremos una
diferencia de concentración local.
Habiendo aclarado lo anterior, procedemos a hacer un balance de masa de A sobre un

elemento diferencial de volumen1 de la fase gaseosa del lecho empacado, que podemos
escribir como
⎛ flujo molar de A que entra en z − ⎞ ⎛ flujo molar de A transferido desde ⎞

⎜ ⎟+⎜ ⎟=0
⎝ flujo molar de A que sale en z + ∆z ⎠ ⎝ la superficie interfacial contenida en ∆V ⎠
(
π R 2 N A, z z − N A, z z +∆ z )+ N A i ∆Ai =0 (6.3-25)
donde NA,z es el flux total de naftaleno en la fase gaseosa en dirección z y NAi es el flux de
naftaleno que se transfiere desde la superficie de las esferas hacia el seno del gas. El
elemento del área interfacial lo obtenemos de multiplicar el área específica del sólido por el
volumen del lecho:
( )
∆ Ai = ( a ) π R 2 ∆z [ = ]
área interfacial
volumen de lecho
⋅ volumen de lecho = área interfacial (6.3-26)
Esta área multiplicada por el flux NAi será el flujo molar de A que se transfiere desde la
superficie interfacial contenida en ∆V.
En flujo turbulento, el flux total incluye las difusiones molecular y turbulenta, y un término
convectivo. Para nuestro ejemplo podemos escribir
1 Como estamos usando la velocidad superficial, tenemos que considerar que el tubo está vacío y entonces el
elemento diferencial de volumen es ∆V=πR2∆z.
R. LOBO/2007
(total)
N A, z = J A, z + CA b vb, z ≅ CA b vb, z (6.3-27)
donde hemos supuesto que flux difusivo total es pequeño comparado con el convectivo. El
flux interfacial está dado por
(
N A i = k c CA i − CA b ) (6.3-28)
Podemos sustituir las Ecs. (6.3-26) a (6.3-28) en la (6.3-25) y obtener:
π R 2 ⎡⎣( CA b vb, z ) z − ( C Ab vb, z ) ⎤

z +∆z ⎦ + kc ( CAi − CAb ) ( a ⋅ π R2∆z ) = 0 (6.3-29)
Si dividimos esta ecuación por πR2∆z, tomamos lím ∆z→0 y rearreglamos, podemos
obtener
d CA b kc a
dz
=
vb, z
(
CA i − CA b ) (6.3-30)
Esta ecuación se encuentra sujeta a una condición de frontera a la entrada del lecho,
z=0; CA b = 0 (6.3-31)
la cual indica que el aire entra al lecho sin naftaleno. Si integramos la Ec.(6.3-30),
CA b ⎛ ka ⎞
= 1 − exp ⎜ − c z ⎟ (6.3-32)
CA i ⎜ v ⎟
⎝ b, z ⎠
donde kc es el coeficiente de transferencia de masa promedio para la evaporación del

naftaleno entre la entrada y la distancia axial z. De esta ecuación podemos despejar el
coeficiente para obtener,
⎛v ⎞ ⎛ CA b ⎞ ⎛ vb, z ⎞ ⎛ pAb ⎞
kc = − ⎜ b , z ⎟ ln ⎜⎜ 1 − ⎟ = −⎜ ln
⎟ ⎜ ⎜ 1 − ⎟ (6.3-33)
⎝ az ⎠ ⎝ CA i ⎟⎠ ⎝ az ⎠ ⎝ pA i ⎟⎠
Si evaluamos esta expresión a la salida del lecho (z=60 cm), el coeficiente de transferencia
de masa promedio sobre todo el lecho es:
200 cm/s ⎛ 2 mm Hg ⎞ −2
kc = − ln ⎜ 1 − ⎟ = 6.81 × 10 cm s (6.3-34)
( 25 cm 2
/ cm
3
) ⋅ ( 60 cm ) ⎝ 5 mm Hg ⎠
R. LOBO/2007
Este es el resultado buscado. El estudiante ya habrá notado las semejanzas de este ejemplo
con el anterior. La metodología seguida es en esencia la misma, pero en el presente ejemplo
hemos introducido algunas variantes debido a la naturaleza del lecho empacado.
Ejemplo 6.3-3. Determinación del coeficiente de transferencia de masa en la

disolución de un sólido
Se desea estudiar la disolución de un sólido puro B. Para ello, se coloca una solución 0.1 M
de un ácido (A) en un recipiente y se agita. A un tiempo determinado (t=0) se añade una
determinada cantidad del sólido en forma de pequeñas partículas, que quedan suspendidas
en la solución. La reacción de disolución es tan rápida que puede considerarse que el ácido
reacciona instantáneamente en la superficie del sólido. Puede suponerse también que
durante el tiempo que dura el experimento el tamaño de las partículas del sólido no cambia
apreciablemente. El avance de la disolución se cuantifica mediante el análisis de la
concentración del ácido en el seno de la solución. De este análisis se obtienen los datos de
la Tabla 6.3-1.
Tabla 6.3-1. Datos experimentales del consumo de ácido en la disolución de un sólido
t CA t CA t CA t CA
(min) (M) (min) (M) (min) (M) (min) (M)
0 0.1 4 0.040 8 0.015 12 0.005
1 0.075 5 0.032 9 0.011
2 0.062 6 0.024 10 0.008
3 0.050 7 0.018 11 0.006
Determinar el coeficiente de transferencia de masa del ácido hacia las partículas del sólido.
Solución. El sistema es muy semejante al ilustrado en la Figura 6.1-2, con la diferencia que
ahora el compuesto de interés es el ácido. El área de transferencia de masa es el área total
de las partículas sólidas que se están disolviendo, aunque no conocemos su valor. Podemos
considerar que la solución está bien agitada y que por ello la concentración del ácido es
uniforme en todo el tanque. La diferencia de concentración del ácido es la existente entre el
seno de la solución y la superficie del sólido; en vista de que la reacción es instantánea, esta
última concentración será cero. Entonces, la diferencia de concentración irá disminuyendo
con el tiempo.
Ya dijimos que la concentración es uniforme en el espacio, por lo que no hay balance de

masa diferencial que realizar, sino un balance sobre el ácido en la solución líquida. Este
balance puede enunciarse como:
( )( )( )(
entrada de A − salida de A + generación ( − consumo) de = acumulación de A
a la solución de la solución A en la solución en la solución )
R. LOBO/2007
A la solución no entra A, por lo que el primer término de esta ecuación es cero. El

compuesto A sale de la solución a través de la superficie interfacial con el sólido, por lo
que este término es igual al flux interfacial de A multiplicado por el área interfacial.
Notemos con cuidado que el ácido no se está consumiendo por reacción en la solución, sino
en la interfase líquido-sólido, por lo que el término de generación del balance es cero.
Finalmente, como resultado de la salida de A en la interfase líquido-sólido, se producirá un
cambio en el número de moles del ácido por unidad de tiempo en la solución. Entonces el
balance de masa sobre A queda como:
(
d V CA b ) = −k
− N A i Ai =
dt
c ( CAb − CAi ) Ai = −kc Ai CAb (6.3-35)
donde NAi es el flux interfacial del ácido que llega a la superficie de las partículas; Ai es el
área del total de las partículas; y V es el volumen de la solución. El término (VCAb)
representa el total de moles de A presente en la solución en cualquier instante, por lo que su
derivada con respecto al tiempo es el cambio en el tiempo de los moles totales de ácido en
la solución. En vista de que el volumen de la solución es constante, podemos sacarlo de la
derivada para obtener:
dCA b ⎛ Ai ⎞
− =⎜ ⎟ kc C Ab = kc a CA b (6.3-36)
dt ⎝V ⎠
donde a=Ai/V es el área superficial interfacial del total de las partículas por unidad de
volumen de solución. La condición inicial a la que está sujeta la Ec.(6.3-36) es:
t=0; CA b = CA0 = 0.1 M (6.3-37)
Podemos integrar la Ec.(6.3-36) sujeta a la (6.3-37) y obtener:
⎛ CA b ⎞
− ln ⎜ ⎟ = ( kc a ) t (6.3-38)
C
⎝ A0 ⎠
donde ( kc a) representa el coeficiente de transferencia volumétrico promedio entre el inicio de

la disolución y el tiempo t. De esta ecuación es evidente que si graficamos los datos
experimentales como −ln(CAb/CA0) vs. t, la pendiente de la línea resultante será el valor de
( kc a) . Es importante notar que no podemos obtener el valor del coeficiente de transferencia
de masa kc , sino únicamente el producto ( kc a) , ya que no conocemos el área interfacial ni el
volumen de la solución. Tangencialmente, debemos señalar que no es recomendable sustituir
un valor de CAb y el correspondiente tiempo en la (6.3-38) para obtener ( kc a) , ya que los
datos experimentales tienen error. El mejor procedimiento es ajustar la mejor línea a los datos
R. LOBO/2007
graficados según la Ec.(6.3-38), y entonces obtener ( kc a) de la pendiente. Esto se muestra en

la Figura 6.3-5.
De la pendiente de la línea de la Figura 6.2-5 obtenemos el resultado deseado:
kc a = 0.249 min −1 = 4.15 × 10−3 s −1 (6.3-39)
3.0
⎛ CA b ⎞
− ln ⎜ ⎟
⎝ CA0 ⎠
2.5
2.0
1.5
pendiente = 0.249 min −1
1.0
0.5
0.0
0 2 4 6 8 10 12
tiempo ( min )
Figura 6.3-5. Obtención del coeficiente de transferencia de masa volumétrico en la disolución de un sólido.
El proceso de disolución pudiera haberse seguido mediante la medición de la concentración

del sólido disuelto en la solución. El cálculo del coeficiente para la disolución del sólido de
lo hubiéramos hecho de manera análoga a la expuesta. En este caso, el área de transferencia
de masa seguiría siendo desconocida y la diferencia de concentración sería la del sólido
disuelto en la interfase sólido/solución (la de saturación) y la prevaleciente en el seno de la
solución. Finalmente conviene recordar nuevamente que al producto (kca) se le denomina
coeficiente de transferencia de masa volumétrico.
6.4 Análisis dimensional y correlaciones

Hasta aquí hemos obtenido un conocimiento básico de la transferencia de masa interfacial.
Este conocimiento será de utilidad para poder hacer uso correcto de la gran cantidad de
literatura existente sobre el tema. Muchas veces el ingeniero se ve confrontado con
problemas prácticos cuya solución requiere de conocer valores del coeficiente de
transferencia de masa. A lo largo del texto hemos planteado y resuelto varios problemas en
los que se ha involucrado el coeficiente de transferencia de masa. En esos casos, siempre
supusimos que era un parámetro conocido aunque, en realidad, es necesario determinarlo
R. LOBO/2007
para obtener la solución. Para el análisis y diseño de equipos de separación de compuestos,

tales como la absorción, la destilación, la extracción líquido-líquido y líquido-sólido, así
como los químicos y bioquímicos, es necesario contar con el conocimiento de valores del
coeficiente de transferencia de masa adecuados a la situación específica. Baste por el
momento decir que este parámetro es fundamental para conocer ya sea la rapidez con que
uno o varios compuestos pueden pasar de una fase a otra, ya sea para determinar los límites
de la rapidez con que es posible llevar a cabo una reacción química o bioquímica.
Por lo anterior, en ésta y la siguiente sección nos dedicaremos a estudiar cómo obtener
valores del coeficiente de transferencia de masa para describir los procesos de transferencia
de masa en diversas situaciones. Una vez que sepamos cómo calcular estos valores, podremos
utilizarlos para el diseño y análisis de equipos como los mencionados. En el Capítulo 8
daremos algunos ejemplos sencillos de cálculo de equipo donde se involucra el conocimiento
del coeficiente de transferencia de masa.
Podemos decir que existen tres métodos para la obtención de coeficientes de transferencia
de masa:
• El primero es la predicción teórica del coeficiente a partir de la teoría de la capa límite.

Este método, aunque ha producido algunos resultados de relevancia práctica, se
encuentra limitado por el tipo de problemas que se han podido resolver. No obstante, el
método tiene mucha relevancia conceptual y lo estudiaremos en la Sec. 7.4.
• El segundo método consiste en recurrir a datos experimentales, ya sean obtenidos por

nosotros o por otros personas, para calcular los coeficientes de transferencia de masa
mediante correlaciones1 de dichos datos. Los datos experimentales de una situación física
particular que se han obtenido para diferentes tipos de solutos, fluidos, velocidades y
geometrías, pueden correlacionarse en términos de números adimensionales. En la Sec.
5.5 hemos discutido algunos aspectos del análisis dimensional que nos ha llevado, por
ejemplo, a identificar a los números adimensionales Re y Sc como de gran de utilidad en
problemas de transferencia de masa. Ahora extenderemos nuestro estudio para incluir una
técnica de análisis dimensional empírica, pero que tiene la ventaja de proporcionar la
forma funcional de las ecuaciones para correlacionar los datos experimentales en términos
de números adimensionales. Este método lo estudiaremos en la presente sección. En la
siguiente sección presentaremos un número de correlaciones basadas en datos
experimentales para estimar el valor del coeficiente de transferencia de masa en diversas
situaciones.
• El tercer método de obtención del coeficiente de transferencia de masa consiste en

utilizar las llamadas analogías entre las transferencias de momentum, calor y masa. El
1 Una correlación es una ecuación empírica que "ajusta" bien muchos datos experimentales. Para nuestro
caso, la variable dependiente es generalmente un número adimensional que contiene al coeficiente de
transferencia de masa y las variables independientes son principalmente Re y Sc, aunque existen otras formas
de correlaciones.
R. LOBO/2007
fundamento de este método está en que las ecuaciones generales que describen los
transportes de masa, momentum y energía son análogas (ver Sec. 5.5) ya que, en ciertas
condiciones, estos tres transportes ocurren por mecanismos similares en el fluido y en la
vecindad de las interfases. Así, es posible utilizar una correlación obtenida para un
problema de transferencia de calor para elaborar otra correlación para obtener el
coeficiente de transferencia de masa para un problema análogo de transferencia de masa.
El método de las analogías para obtener el coeficiente de transferencia de masa lo
estudiaremos en la Sec. 6.6.
Aunque ya tenemos antecedentes sobre algunos de los números adimensionales relevantes

en transferencia de masa, es conveniente extender nuestra discusión aquí.
A. Números adimensionales importantes en transferencia de masa
Anteriormente hemos visto que los números adimensionales generalmente tienen una
interpretación física, ya que representan cocientes de propiedades o de rapideces de
procesos. Aunque resulte evidente, es de vital importancia comprender con toda precisión
a qué se refieren cada una de las variables, parámetros y propiedades que entran en la
definición de cada número adimensional, ya que no siempre significan lo mismo, sino una
situación física determinada. Ejemplifiquemos. Como ya señalamos, la descripción del
flujo en la vecindad de las interfases es uno de los problemas principales en la predicción
del coeficiente de transferencia de masa y el principal parámetro que caracteriza el flujo por
convección forzada es el número de Reynolds, Re
L0 ρ v0 fuerzas inerciales
Re L0 = [ = ] fuerzas viscosas (6.4-1)
µ
Observemos que ahora Re tiene un subíndice. L0 es una distancia característica que se

emplea para calcular la magnitud de Re; esta distancia característica es alguna dimensión
importante del sistema bajo consideración. Por ejemplo, para flujo alrededor de una esfera,
L0 es el diámetro de ésta; para flujo sobre una placa plana, L0 es la longitud de la placa; para
flujo dentro de un tubo cilíndrico, L0 puede ser el diámetro del tubo. En adelante, el número
de Reynolds siempre lo emplearemos con algún subíndice, que indicará la distancia
característica empleada. La velocidad característica v0 también puede referirse a varias
cosas: en un tubo vacío, v0 puede ser la velocidad másica promedio del fluido; en un lecho
empacado, v0 puede ser la velocidad superficial vsup; o bien puede ser la velocidad
intersticial vi=vsup/ε, donde ε es la fracción volumétrica hueca del lecho; en una película de
líquido que desciende en flujo laminar v0 puede ser la velocidad promedio o la velocidad
máxima del líquido que desciende. Como podemos ver, cada vez que usemos un número
adimensional debemos ser precisos en cuanto el significado de los elementos que lo
componen.
R. LOBO/2007
Cuando existen diferencias de densidad apreciables dentro de una fase dada que dan origen
a la convección natural, el número de Grashof es importante y se define por
gL30 ( ∆ρ ρ ) fuerzas boyantes

Gr = [ = ] fuerzas viscosas (6.4-2)
(µ ρ ) 2
Otro número adimensional importante en transferencia de masa ya conocido por nosotros es

el número Schmidt, que se define por
ν µ rapidez del transporte difusivo de momentum

Sc = = [=] (6.4-3)
DAB ρ DAB rapidez del transporte difusivo de masa
El número de Schmidt sólo contiene propiedades de la mezcla fluida donde ocurre la

transferencia de masa. La viscosidad cinemática ν así como las propiedades que la
componen, la viscosidad µ y la densidad ρ son las de la mezcla. Es importante enfatizar esto,
ya que la difusividad, la viscosidad y la densidad pueden variar apreciablemente con la
concentración y la temperatura. No obstante, cuando las mezclas son diluidas, sus
propiedades físicas pueden aproximarse por las del solvente puro.
En el Capítulo 5 vimos que el número de Peclet de transferencia de masa a menudo aparece

como variable independiente en correlaciones. Al igual que con el número de Reynolds
aquí también existe la necesidad de definir claramente la distancia característica L0:
L0 v0 rapidez del transporte de masa por convección

( Pem ) L0 = Re L0 Sc = [=] (6.4-4)
DAB rapidez del transporte de masa por difusión
Cuando existen situaciones en las que hay transferencia simultánea de calor y masa,
frecuentemente aparece el número de Lewis definido como
α kT kT rapidez del transporte difusivo de calor

Le = =
ˆ
=
[ =] rapidez del transporte difusivo de masa (6.4-5)
DAB ρ C p DAB CC p DAB
donde α es la difusividad térmica, kT es la conductividad térmica y Cp es la capacidad

calorífica, ya sea en unidades másicas o molares, todas propiedades de la mezcla.
El número de Sherwood es el parámetro adimensional que contiene al coeficiente de

transferencia de masa más común; se utiliza ampliamente como variable dependiente en las
correlaciones, y se define por:
R. LOBO/2007
kc L0 L0 DAB resistencia al transporte de masa por difusión

Sh L0 = = [ =] resistencia el transporte de masa por convección
DAB 1 kc
(6.4-6)
rapidez del transporte de masa por convección
=
rapidez del transporte de masa por difusión
Nuevamente, debemos observar que L0 representa una distancia característica particular de

cada situación física (cuando este número adimensional se correlaciona con el de Reynolds,
debe usarse en ambos la misma distancia característica). Debe ser claro, por ejemplo que
ShL0=3 significa cosas diferentes para la transferencia de masa desde burbujas de un gas a
un líquido, y para la disolución de un sólido desde una placa plana hacia un líquido en flujo
turbulento.
Debemos tener cuidado de no confundir el número de Sherwood con el número de Biot de

transferencia de masa, Bim, introducido en la Sec. 2.5, ya que ambos tienen una definición
semejante. El número de Biot apareció en relación a una interfase sólido-fluido como una
razón de la resistencia al transporte difusivo dentro del sólido a la resistencia al transporte
convectivo en la fase fluida; en consecuencia la difusividad DAB en el número de Biot se
refiere al soluto en el sólido. En contraste, y como lo indica la Ec.(6.4-6), el número de
Sherwood indica una razón de resistencias al transporte difusivo y convectivo, en la fase
fluida donde ocurre la transferencia y por ello DAB es la difusividad del soluto en la fase
fluida.
Hay otros números adimensionales usados en problemas de transferencia de masa y algunos

de ellos pueden obtenerse por combinación de algunos de los hemos definido aquí. Dos de
ellos son particularmente importantes porque en su definición aparece el coeficiente de
transferencia de masa y también se usan con frecuencia como variables dependientes en
correlaciones y analogías. El número de Stanton se define por
Sh L 0 kc
St = = (6.4-7)
Re L 0 Sc v0
Como tendremos ocasión de ver más adelante, a menudo se utiliza el denominado factor jD
para correlacionar datos experimentales del coeficiente de transferencia de masa. Este
factor se define por
kc 2 3
jD = St Sc 2 3 = Sc (6.4-8)
v0
y su relación con el número de Sherwood es:
R. LOBO/2007
Sh L 0 = jD Re L 0 Sc1 3 (6.4-9)
En la Tabla 6.4-1 se sintetizan los números adimensionales relevantes para correlacionar

coeficientes de transferencia de masa.
R. LOBO/2007
310
Tabla 6.4-1. Números adimensionales de uso común en transferencia de masa.
Número Significado físico Observaciones

adimensional
Reynolds fuerzas inerciales Poner cuidado en la definición de la distancia característica y velocidad. Se usa en
convección forzada.
Re L0 = L0 v0 ν fuerzas viscosas
Schmidt difusividad de momentum
Sc = ν DAB Contiene propiedades físicas de la mezcla. Aparece en analogías y correlaciones.
difusividad de masa
Sherwood rap. de transferencia de masa Contiene al coeficiente de transferencia de masa. Variable dependiente en
correlaciones. También llamado No. de Nusselt de transferencia de masa
rap. de difusión
Sh L0 = kc L0 DAB
Stanton rap. de transferencia de masa Usado en la analogía de Reynolds.
St = kc v0 St= Sh L Re L Sc
rap. de flujo 0 0
Factor jD
--- Analogía de Chilton-Colburn
jD = St Sc2 3
Lewis difusividad de calor
Le = α DAB Uso frecuente en problemas de transferencia de calor y masa simultáneos.
difusividad de masa
Peclet Masa rap. transp. masa convección ( Pem ) L = Re L Sc
0 0
PRINCIPIOS DE TRANSFERENCIA DE MASA/ 2a. EDICIÓN
Pem = v0 L0 DAB rap. transp. masa difusión

Biot Masa rap. de transferencia de masa externa Distinguir de Sh. L0 y DAB son de la fase sólida interna. kc rige la transferencia de
( Bim ) L0 = kc L0 DAB rap. de difusión interna en el sólido masa al fluido externo al sólido. Ver definición en Sec. 2.5.
Damköhler II rapidez de la reacción
Importante en problemas de difusión y reacción química simultáneas.
Da II = k R CAn −1 L20 DAB rapidez de la difusión
Grashof fuerzas boyantes
2 Importante en convección natural
Gr = gL30 ( ∆ρ ρ ) ν fuerzas viscosas
R. LOBO/2007
B. Análisis dimensional
En la Sec. 5.5 realizamos el análisis dimensional de las ecuaciones de cambio y logramos

identificar algunos de los parámetros más relevantes en las transferencias de momentum,
masa y calor. En este apartado extenderemos nuestro estudio del análisis dimensional para
obtener la forma general de las correlaciones.
El teorema π de Buckingham1, es la base del análisis dimensional2 y permite determinar los

grupos adimensionales relevantes para la descripción empírica de determinados procesos.
Este teorema establece que
la relación funcional entre n variables cuyas dimensiones puedan escribirse en

términos de m dimensiones fundamentales, puede expresarse como función de (n−m)
grupos adimensionales compuestos de combinaciones de las variables.
Aunque el teorema señala el número de grupos adimensionales, no informa sobre la

funcionalidad de éstos. El método de los índices de Rayleigh supera esta dificultad, por lo que
será el que usemos para ejemplificar la forma de obtener correlaciones del coeficiente de
transferencia de masa.
Supongamos que tenemos un sistema descrito por las variables x0, x1, ..., xn . Entonces,
podemos suponer que el comportamiento del sistema podrá ser descrito por alguna relación
de la forma
x0 = Φ( x1 , x2 ,..., xn ) (6.4-10)
Rayleigh supuso que la función Φ podía ser expresada por una serie de potencias infinita de
las variables independientes de la forma
i =∞
y y y
x0 = ∑ Ai x1 1i x2 2 i ...xn ni (6.4-11)
i =−∞
donde las Ai son coeficientes numéricos y yki son índices numéricos. Esta expresión
implica de hecho la continuidad de x0 y sus derivadas parciales con respecto a las otras
variables, condición que generalmente se cumple. La Ec.(6.4-11) representa cantidades
físicas que tienen dimensiones. Entonces, todos los términos de la serie deben tener las
mismas dimensiones o, dicho de otra manera, la ecuación debe ser dimensionalmente
homogénea. Esto impone restricciones sobre los valores que pueden tener los índices y1i
, ..., yni en la serie infinita. Para encontrar estas restricciones, sólo necesitamos
considerar un término típico de la serie.
1 Buckingham, E. (1914) Phys. Rev., 2, 345.

2 McAdams, WH. (1964) Transmisión de Calor. 3ª. Edición. McGraw Hill. Madrid
R. LOBO/2007
Veamos ahora con un ejemplo cómo trabaja el método de los índices de Rayleigh.
Consideremos un tubo cilíndrico largo recubierto en su pared interna de un material A
ligeramente soluble, dentro del cual fluye un líquido en flujo turbulento al que se está
transfiriendo el soluto. El sistema es semejante al ilustrado en la Figura 6.3-2. El
problema consiste en conocer cómo está relacionado el coeficiente de transferencia de
masa promedio kc que rige la disolución del recubrimiento del tubo con el resto de las
variables del sistema. Dividamos el procedimiento en varios pasos.
Primer paso: Seleccionar las variables de las cuales depende el coeficiente de transferencia
de masa.
a. Si consideramos un tubo lo suficientemente y suponemos que el flujo está

completamente desarrollado, podemos ignorar los efectos de entrada; entonces, la
longitud del tubo no afectará. Sin embargo, el diámetro del tubo D será importante.
b. El flujo está determinado por la velocidad del fluido vb,z, por su densidad ρ y por su
viscosidad µ.
c. Es razonable suponer que en la pared del tubo la velocidad del fluido es cero, por lo que
en la interfase fluido-pared la transferencia de masa se da únicamente por difusión. En
consecuencia, la difusividad DAB del soluto en la mezcla será importante.
Entonces, el coeficiente de transferencia de masa será una función de la forma
kc = Φ (vb, z , ρ , µ , DAB , D) (6.4-12)
Únicamente necesitamos un término de la serie para determinar los exponentes. A fin de

simplificar la notación usaremos para éstos a, b, etc. Por consiguiente, la forma de la
Ec.(6.4-12) será
kc = Avba, z ρ b µ c DABd D e (6.4-13)
Segundo paso: Seleccionar las dimensiones fundamentales en que se expresarán las

variables.
Las dimensiones de cada una de las variables de la Ec. (6.4-13) son:
L L M
kc [ = ] ; A [ = ] adimensional ; vb, z [ = ] ; ρ [ =]
t t L3
(6.4-14)
M L2
µ [ =] ; DAB [ = ] ; D [ =] L
Lt t
R. LOBO/2007
La inspección de las Ecs.(6.4-14) nos permite concluir que las dimensiones relevantes en el
presente caso son la masa [M], la longitud [L] y el tiempo [t].
Tercer paso: Encontrar los grupos adimensionales.
En primer lugar, usemos el teorema π. El número de variables n, de acuerdo a la Ec.(6.4-13)

es seis; el número de dimensiones fundamentales presentes en ellas es m=3, tal como
acabamos de concluir. Entonces, el número de grupos adimensionales necesarios para
describir el sistema es n−m=3. Para encontrar cuáles son dichos grupos adimensionales,
construyamos una matriz dimensional (Tabla 6.4-2) donde se expresan las dimensiones de
cada variable presente en la Ec.(6.4-13), como se muestra a continuación.
Tabla 6.4-2. Matriz de dimensiones de las variables que afectan al coeficiente

de transferencia de masa.
Variable →
kc vb, z ρ µ DAB D
Dimensiones ↓
[M] 0 0 1 1 0 0
[L] 1 1 −3 −1 2 1
[t] −1 −1 0 −1 1 0
Con la ayuda de esta matriz, es posible escribir la Ec.(6.4-13) en forma dimensional
a b c d
⎛L⎞ ⎛M⎞ ⎛ M ⎞ ⎛L ⎞
2
⎛L⎞
⎜ ⎟ [ =] A ⎜ ⎟ ⎜ 3 ⎟ ⎜ ⎟ ⎜ ⎟ ( L )
e
(6.4-15)
⎝t⎠ ⎝ t ⎠ ⎝ L ⎠ ⎝ Lt ⎠ ⎝ t ⎠
Si igualamos los exponentes de cada una de las dimensiones, obtenemos
Para [ M ] 0 =b+c
Para [ L ] 1 = a − 3b − c + 2d + e (6.4-16)
Para [ t ] − 1 = −a − c − d
El sistema de las Ecs.(6.4-16) consta de 3 ecuaciones y 5 incógnitas, por lo que es posible

resolver para 3 incógnitas en términos de las dos restantes. Por ejemplo, si escogemos
resolver en términos de a y d obtenemos
b = a + d −1
c = a−a−d (6.4-17)
e = a −1
R. LOBO/2007
Si sustituimos estos resultados en la Ec.(6.4-13) y rearreglamos términos podemos obtener
a −d
kc ρ D ⎛ ρ D vb, z ⎞ ⎛ µ ⎞
= A⎜ ⎟ ⎜ ⎟ (6.4-18)
µ ⎝ µ ⎠ ⎝ ρ DAB ⎠
El grupo del lado izquierdo de esta ecuación contiene al coeficiente de transferencia de

masa y es adimensional. Sin embargo, es uno que no se usa comúnmente. Si multiplicamos
ambos lados de la Ec.(6.4-18) por el número de Schmidt Sc=µ/ρDAB, después de rearreglar
llegaremos a:
a 1− d
k D ⎛ Dv ρ ⎞ ⎛ µ ⎞ 1− d
= A ( Re D ) ( Sc )
a
Sh D = c = A ⎜ b, z ⎟ ⎜ ⎟ (6.4-19)
DAB ⎝ µ ⎠ ⎝ ρ DAB ⎠
Podemos reconocer en esta expresión al grupo que contiene el coeficiente de

transferencia de masa como el número de Sherwood, donde la distancia característica
usada es el diámetro del tubo D. Asimismo, reconocemos el número de Reynolds basado
en la misma distancia característica y el número de Schmidt. La Ec.(6.4-19) es la
correlación deseada; el coeficiente A y los exponentes a y d se determinan a partir de
datos experimentales. Para el presente caso, a partir de datos experimentales, Linton y
Sherwood1 determinaron el coeficiente A y los exponentes para finalmente obtener
0.83 13
kD ⎛ D vb, z ρ ⎞ ⎛ µ ⎞
Sh D = c = 0.023 ⎜ ⎟ ⎜ ⎟ = 0.023Re0.83
D Sc
13
DAB ⎝ µ ⎠ ⎝ ρD ⎠
2000 < Re D < 70000 (6.4-20)
1000 < Sc < 2 260
Es fundamental darnos cuenta que todas las correlaciones tienen rangos de validez en
cuanto los valores de las variables independientes, que no son otros que los rangos de
valores de dichas variables en los experimentos en que se basan las correlaciones. El uso de
las correlaciones para estimar la variable dependiente fuera de dichos rangos no es
recomendable. Una correlación está completa sólo cuando junto con ella aparecen los
rangos de validez de los números adimensionales usados como variables independientes.
Es evidente que podríamos haber decidido resolver el sistema de Ecs.(4.3-9) en términos de

otros dos exponentes, y haber obtenido otra correlación en términos de grupos
adimensionales distintos; en este caso, la nueva correlación hubiera sido totalmente
equivalente a la expresada en la Ec.(4.3-12). Sin embargo, es preferible expresar los
resultados en términos de números adimensionales comúnmente utilizados.
1 Linton, WH y Sherwood, TK (1950) Chem. Engr. Prog., 46, 258.
R. LOBO/2007
Como ya hemos señalado, el análisis dimensional permite una planeación más racional de
experimentos con un número más pequeño de variables que las que afectan al sistema de
interés. En el ejemplo recién visto, son 5 las variables (vb,z, ρ, µ, DAB, D) que afectan al
coeficiente de transferencia de masa , pero gracias al análisis dimensional hemos logrado
reducirlas a dos (ReD, Sc). Debe ser evidente que necesitamos muchos menos experimentos
para caracterizar a kc usando dos variables independientes que cinco. Además, hemos
logrado otra ventaja adicional: para variar los números adimensionales, es posible
manipular la variable dentro de ellos que sea más sencilla de cambiar. También, en caso de
que los materiales de que está hecho el sistema de interés sean caros, pueden utilizarse en
su lugar materiales baratos, ya que lo que importa es el valor del número adimensional y no
el valor de las propiedades específicas del sistema en cuestión.
No obstante su sencillez, este método empírico del análisis dimensional tiene limitaciones. La
principal de ellas es la selección de las variables que afectan al sistema. Si no comprendemos
bien el problema, es probable que no seleccionemos adecuadamente todas las variables
significativas, quizá omitamos alguna importante o incluyamos una irrelevante. También
puede darse el caso de incluir variables cuyo efecto esté implícito en otras variables, como
sería la temperatura en nuestro ejemplo (la temperatura afecta propiedades como la densidad
y la viscosidad). Otra limitación es la dificultad en la selección apropiada de las dimensiones;
normalmente, podemos usar las dimensiones convencionales, pero en ocasiones las ideas
implícitas en ellas no son las apropiadas. Por ello, es preferible el método del análisis
dimensional de las ecuaciones de cambio expuesto en la Sec. 5.5 para encontrar los números
adimensionales relevantes, en combinación con el método de Rayleigh para encontrar la
forma funcional de la correlación.
Para mayores detalles sobre la aplicación de los métodos del análisis dimensional el
estudiante puede consultar los textos de Johnstone y Thring1 y Broadkey y Hershey2.
6.5 Correlaciones del coeficiente de transferencia de masa

Hasta aquí hemos mostrado cómo obtener valores del coeficiente de transferencia de masa
que pueden derivarse de las expresiones analíticas del flux interfacial de masa cuando el
flujo es laminar. En la sección anterior, mostramos que mediante experimentos y el análisis
dimensional es posible obtener correlaciones empíricas del coeficiente de transferencia de
masa para situaciones donde el flujo es laminar y, particularmente, cuando el flujo es
turbulento. En la presente sección expondremos algunas de las múltiples correlaciones que
se han publicado en la literatura especializada y que han mostrado su utilidad práctica para
1 Johnstone, RE. y Thring, MW (1957) Pilot Plants, Models, and Scale-up Methods in Chemical Engineering,
McGraw-Hill, New York.
2 Broadkey, RS y Hershey HC (1988) Transport Phenomena. A Unified Approach. McGraw-Hill, New York.
R. LOBO/2007
estimar valores del coeficiente de transferencia de masa con niveles aceptables de

incertidumbre1.
Por otra parte, debemos tener en cuenta que las correlaciones se han desarrollado con datos
experimentales provenientes de muy diversas investigaciones y muchas de ellas se
publicaron hace tiempo, con los datos disponibles hasta ese momento. Por ello es
recomendable que el estudiante acuda a diversos trabajos de revisión crítica de la
investigación en el campo que se publican ya sea en artículos o libros2. Cabe señalar que en
las correlaciones que se presentan a continuación el coeficiente de transferencia de masa es
uno promedio; cuando sea local, así lo haremos notar.
A. Superficies planas
El caso considerado es el de la transferencia de masa desde o hacia una superficie plana

conforme un fluido fluye paralelo a la superficie. La capa límite que se forma sobre la
superficie puede ser únicamente laminar, o puede presentar zonas de flujo laminar, de
transición y de flujo turbulento, o ser únicamente turbulenta. La mayor parte de los datos
experimentales se han obtenido en gases midiendo la rapidez de evaporación de un líquido
o la sublimación de un sólido. Las correlaciones son:
kc z
Flujo laminar: Sh z = = 0.332Re1z 2 Sc1 3 (6.5-1)
DAB
kc Lz
Flujo laminar: Sh Lz = = 0.664Re1Lz2 Sc1 3 (6.5-2)
DAB
Re z y Re Lz < 2 × 105 ; 0.6<Sc<2500
donde Lz es la longitud de la superficie. Los números de Reynolds se definen por :
ρ vb, z z ρ vb, z Lz
Re z = ; Re Lz = (6.5-3)
µ µ
1 La situación física implícita en cada correlación presenta una fenomenología compleja que éstas solo pueden
representar de manera aproximada. Para cálculos de ingeniería, una guía sobre el nivel aceptable de error es
de ± 20 %, aunque algunas correlaciones muestran desviaciones de los valores experimentales hasta del 40 %.
Otras, en cambio, tienen un error de ± 10 %.
2 Las revisiones críticas son importantes porque son elaboradas por expertos en el campo, discuten la
fenomenología específica de la transferencia de masa, contrastan los resultados los experimentales de diversos
autores y las circunstancias experimentales bajo las cuales las correlaciones se han obtenido, y plantean las
limitaciones de éstas. Ver, por ejemplo: Sherwood, TK, Pigford, RL y Wilke, CR (1975) Mass Transfer.
McGraw Hill, New York. Geankoplis, CJ (1972) Mass Transport Phenomena. Holt, Rinehart and Winston.
New York.
R. LOBO/2007
y vb,z es la velocidad promedio de la corriente libre. Observemos que la Ec.(6.5-1)

proporciona el valor de un coeficiente de transferencia de masa local a la longitud z medida
desde el inicio de la superficie, en tanto que la Ec. (6.5-2) da el coeficiente promedio sobre
la longitud Lz de la superficie plana donde existe el flujo laminar.
Para flujo turbulento sobre una superficie plana también existen correlaciones para calcular
el coeficiente local y el coeficiente promedio. Las correlacione para el coeficiente local y
promedio son:
kc z
Flujo turbulento: Sh z = = 0.0296 Re0.8
z Sc
13
(6.5-4)
DAB
kc Lz
Flujo turbulento: Sh Lz = = 0.0365Re0.8
Lz Sc
13
(6.5-5)
DAB
R e Lz > 2 × 105 ; 0.6<Sc<2500
La especificación de ReLz >2×105 para considerar que el flujo es turbulento es un tanto

arbitraria, ya que la zona de transición está delimitada ente este valor de Re y 5×105. La
Ec.(6.5-4) debe usarse particularmente cuando sobre una placa ocurren ambos flujos
laminar y turbulento; la distancia z en esta correlación debe ser contada a partir del
establecimiento del flujo turbulento. Cuando la longitud total de la superficie es mayor que
la distancia crítica zc a la cual ocurre la transición a flujo turbulento, el coeficiente
promedio sobre toda la longitud de la placa debe ser algún tipo de promedio pesado de los
valores del coeficiente laminar y el turbulento, como se ilustrará en el Ejemplo 6.5-1. Si la
longitud total de la superficie plana es mucho mayor que la distancia crítica a la cual el
flujo se torna turbulento (Lz>>zc) debe usarse la Ec.(6.5-5).
B. Películas líquidas descendentes y flujo en tubos
Aquí consideraremos dos casos. (a) absorción de un soluto hacia o desorción del mismo
desde la película líquida descendente, donde la transferencia de masa es controlada por la
resistencia en el líquido. El proceso de lleva a cabo en una columna de pared mojada donde
el líquido desciende en flujo laminar en forma de una película estable, para tiempos de
contacto gas-líquido cortos [ver Ec.(6.2-23)]. Las columnas de pared mojada son muy
ubicuas en estudios de transferencia de masa debido a que el área interfacial de
transferencia es fácilmente discernible de su geometría (ver ejemplo 6.5-2)
La correlación a usar es teórica y fue desarrollada en la Sec. 6.2, aunque aquí la

presentamos con una pequeña modificación en la definición del número de Reynolds, que
también impacta de definición del número de Péclet:
R. LOBO/2007
kc Lz
( ) (L ) ( ) (L )
12 12 12 12
Sh Lz = = 0.691 Re Lz Sc z Ly = 0.691 Pem Lz z Ly (6.5-6)
DAB
donde
4Γ 4w 4 ρ vb, y Lz 4vb, y Lz
Re Lz = = = ; Pem,Lz = ; Re Lz < 1200 (6.5-7)
µ π Dµ µ DAB
En estas ecuaciones, w es el flujo másico del solvente, Γ es el flujo másico del solvente por
unidad de perímetro de la columna (M/tL), D es el diámetro interno de la columna, Lz es el
grosor de la película, vb,y es la velocidad promedio de la película líquida descendente(=〈vy〉
= 2/3 vmax), y el número de Reynolds está basado en el grosor de la película líquida Lz, el
cual puede calcularse mediante:
13
⎛ 3ν Γ ⎞
Lz = ⎜ ⎟ (6.5-8)
⎝ ρg ⎠
donde g es la aceleración de la gravedad. El límite superior del número de Reynolds a usar

es aún materia de discusión, ya que diversos autores señalan distintos valores en los cuales
el flujo cesa de ser laminar y se torna turbulento. Aquí usamos la recomendación de
Geankoplis1 señalando que los valores del coeficiente obtenidos de estas ecuaciones deben
ser multiplicados por un factor de 1.5 para obtener un estimado más cercano a los datos
experimentales con que se cuenta. Para números de Reynolds entre 1200 y 3000, se
recomienda multiplicar el valor del coeficiente obtenido por un factor de entre 1.5 y 2.5.
Para tiempos de contacto largos, Sherwood y col. recomiendan una aproximación que se
obtiene de la solución del problema para tiempos largos expresada por la Ec.(4.6-17)
kc L y ρ ⎛ 2DAB ρ Ly ⎞
≈ 0.24 + 5.1⎜ ⎟ (6.5-9)
Γ ⎝ 3ΓLz ⎠
(b) El segundo caso considera la evaporación de un líquido desde una película líquida en
una columna de pared mojada, y la evaporación (o disolución) de un sólido desde la pared
de un tubo, hacia un fluido en el cual se localiza la resistencia a la transferencia de masa y
se encuentra en flujo turbulento. La siguiente correlación la obtuvieron Gilliland y
Sherwood2 y es aplicable para gases:
1 Geankoplis, CJ (1972) Mass Transport Phenomena. Holt, Rinehart and Winston. New York.
2 Gilliland, ER y Sherwood, TK(1934) Ind. Eng. Chem., 26,516.
R. LOBO/2007
kc D pBM
gases: Sh D = = 0.023Re0.83
D Sc
0.44
DAB P (6.5-10)
2, 000 < Re D < 35, 000 ; 0.6 < Sc < 2.5
En esta ecuación (pBM/P) es un factor de corrección para altas concentraciones o altas tasas
de transferencia de masa semejante a (yB)ln en la Ec.(6.2-7). En mezclas diluidas el factor es
igual a la unidad.
Para el caso de líquidos, aplica la correlación de Linton y Sherwood1:
kc D
líquidos: Sh D = = 0.023Re0.83
D Sc
13
DAB (6.5-11)
2, 000 < Re D < 35, 000 ; 1,000 < Sc < 2, 260
Los datos en que se basan las correlaciones expresadas por las Ecs.(6.5-10) y (6.5-11) han
sido también correlacionados mediante el factor jD para gases y líquidos:
kc
(Sc ) = 0.023Re−D0.17
23
jD =
vb, z
(6.5-12)
Re D > 2100; 0.6 < Sc < 3000
En las correlaciones (6.5-10), (6.5-11) y (6.5-12) el número de Reynolds está definido por:
ρ vb, z D
Re D = (6.5-13)
µ
donde la distancia característica D es el diámetro del tubo y la velocidad es la velocidad

axial promedio vb,z a través del área seccional del tubo.
C. Esferas
La situación física de las correlaciones que se expondrán aquí implica el flujo de un líquido
o un gas sobre una esfera rígida estática. Esto se puede dar en una variedad de situaciones:
(a) un sólido puro en contacto con un gas o un líquido en flujo; (b) una gota esférica de
líquido puro en contacto con una gas en flujo; y (c) una gota esférica de líquido dispersa en
otro líquido. La mayor parte de los datos experimentales se han obtenido en la evaporación
de gotas de un líquido puro hacia un gas o por analogía de transferencia de calor desde
esferas sólidas hacia un gas o a un líquido. Las correlaciones para esta situación son
diversas, por lo que se recomienda al estudiante poner atención en los detalles de cada una.
1 Linton WH y Sherwood, TK (1950) Chem. Eng. Progr., 46, 258.
R. LOBO/2007
Aparentemente, la correlación más completa, basada en datos de transferencia de masa en

gases y líquidos, es la de Steinberger y Treybal1:
kc De
( )
0.62
Sh De = = Sh 0 + 0.347 Re De Sc1 2 (6.5-14)
DAB
donde
Sh 0 = 2.0 + 0.569 ( GrSc )

0.25
(6.5-15)
Los rangos dentro de los cuales esta correlación es válida son:
1 < Re De < 16900 ; 0.6<Sc<3200 ; (GrSc)<108
Esta correlación incluye el efecto de la convección natural alrededor de la esfera y tiene una
desviación promedio de los datos experimentales ±13 %. Sin embargo, a menudo no se
cuenta con los datos necesarios para calcular el efecto de la convección natural
cuantificados por el número de Grashof (Gr). En estos casos, para gases puede utilizarse la
bien conocida ecuación de Frössling2 que, con una pequeña modificación de la
originalmente publicada, es
kc De
gases Sh De = = 2.0 + 0.552 Re0.53
De Sc
13
DAB (6.5-16)
4
1 < Re De < 4.8 × 10 ; 0.6<Sc<2.7
Para líquidos, Geankoplis3 recomienda:
kc De
líquidos Sh De = = 0.347 Re0.62
De Sc
13
DAB (6.5-17)
2000 < Re De < 16900; Sc<3200
Existen otras correlaciones probadas, aplicables principalmente a líquidos, que fueron

originalmente correlacionadas al número de Péclet de transferencia de masa y que pueden
ser usadas con confianza dentro de los rangos indicados. Una de ellas es la de Brian and
Hales4
1 Steinberger, RL y Treybal, RE (1960) AIChE J., 6, 227. Sherwood y col. (1975) sostienen que la mejor
correlación es la de Ranz,WE y Marshall Jr, WR (1952). Chem. Eng. Progr., 48, 141 y 173.
2 Frössling, N (1938) Beitr. Geophys., 52, 170.
3 Geankoplis, CJ (1972) Mass Transport Phenomena. Holt, Rinehart and Winston. New York.
4 Brian, PLT y Hales, HB (1969) AIChE J., 15, 419.
R. LOBO/2007
kc De
( )
12
Sh De = = 4 + Pe 2m,3De
DAB (6.5-18)
( )
Pe m,De = Re De Sc < 1× 104
Para el caso en que los números de Péclet sean más grandes, Levich1 recomienda usar:
kc De
Sh De = = 1.01Pe1m,3De
DAB (6.5-19)
( )
Pem,De = Re De Sc > 1× 10 4
No está por demás reiterar que en todas estas correlaciones la distancia característica para el
cálculo de los números de Reynolds, Sherwood y Péclet de transferencia de masa es el
diámetro de la esfera De, que la velocidad a utilizar es la velocidad promedio del fluido vb,z
perpendicular a la esfera y que los coeficientes de transferencia de masa son el promedio
sobre toda el área de la esfera.
En todas las correlaciones descritas el tamaño de la partícula esférica se considera

constante. Sin embargo, en muchas situaciones prácticas tales como la evaporación de gotas
o la disolución de partículas sólidas el tamaño cambia y el transitorio debe tomarse en
cuenta y quizá sea necesario introducir una corrección por altas concentraciones o altas
tasas de transferencia de masa semejante a (yB)ln.
D. Burbujas y gotas
Las gotas y burbujas pequeñas hasta de 0.1 cm de diámetro se comportan como esferas
rígidas. Cuando la resistencia a la transferencia de masa se encuentra en el exterior de las
partículas (en la fase continua), el coeficiente de transferencia de masa puede ser calculado
de las correlaciones del apartado anterior y en particular de la Ec. (6.5-19), aunque para
burbujas, la siguiente correlación empírica2 da mejores resultados
kc Db ⎛ Db ⎞
Sh Db = = 1.13Pe1m,2Db ⎜ ⎟ (6.5-20)
DAB ⎝ 0.45 + 0.2 Db ⎠
donde Db es el diámetro de la burbuja y se expresa en cm.
Para burbujas y gotas de tamaño moderado, los fenómenos se tornan mucho más
complejos, ya que al interior de las mismas existe una circulación del fluido presente, lo
1 Levich, VG (1962) Physicochemical Hydrodynamics. Prentice Hall. Englewood Cliffs.

2 Johnson,AI, Besic, F y Hamilec,AE (1969). Can. J. Chem. Eng. 47, 559.
R. LOBO/2007
cual provoca un aumento en la rapidez del transporte. Las burbujas y gotas de mayor
tamaño aún, oscilan en su forma y con frecuencia se rompen y también coalescen. Además,
la presencia de impurezas que modifican la tensión superficial disminuye la rapidez de la
transferencia. Estos y otros factores hacen que los datos experimentales existentes
presenten grandes variaciones y que la predicción de los coeficientes de transferencia de
masa sea difícil.
E. Cilíndros
La situación física representada por las correlaciones corresponde a cilindro fijo en el

espacio sobre el cual fluye un fluido perpendicular su eje axial, ignorando la transferencia
de masa en los extremos del cilindro. Los datos experimentales provienen de la sublimación
de sólidos en una corriente gaseosa y pocos de sistemas líquidos; sin embargo, se ha
mostrado una considerable concordancia con datos análogos de transferencia de calor. Las
correlaciones son:
kc D
gases: Sh D = = 0.281Re0.50
D Sc
0.44
DAB
(6.5-21)
400 < Re D < 25, 000 ; 0.6 < Sc < 2.6
Para el caso de líquidos, Geankoplis recomienda la correlación de Linton y Sherwood:
kc
(Sc ) = 0.281Re−D0.4
23
líquidos: jD =
vb
(6.5-22)
400 < Re D < 25, 000 ; 3, 000 < Sc
F. Lechos empacados y fluidizados
La situación física considerada enseguida es la de un lecho empacado con sólidos de forma

regular o irregular a través de la cual fluye un gas o un líquido (una sola fase), donde la
resistencia a la transferencia de masa se localiza en el fluido. Las aplicaciones involucradas
en esta situación física incluyen columnas de adsorción y desorción de gases y líquidos en
sólidos, columnas de intercambio iónico y reactores catalíticos gas-sólido. Las
correlaciones dadas aquí no aplican para flujo de dos fases fluidas en un lecho empacado.
Para la correcta aplicación de las correlaciones debemos definir algunos conceptos relativos
a las características físicas de los lechos empacados, y a la velocidad del fluido dentro de
ellos.
Los lechos empacados contienen un volumen vacío y un volumen sólido ocupado por el
empaque. En esta estructura se forman canales con áreas seccionales variables a través de
R. LOBO/2007
los cuales fluye el fluido. El volumen vacío se define en términos de la fracción vacía,
fracción de huecos del lecho o porosidad del lecho:
volumen de huecos
ε = fracción vacía = (6.5-23)
volumen total del lecho
El área de la superficie del empaque en el lecho se expresa en términos su área específica:
área del empaque

a = área específica del lecho = (6.5-24)
volumen del lecho
Puede entonces definirse un equivalente del radio hidráulico para el flujo del fluido en un
lecho empacado como:
volumen de huecos volumen de lecho ε

Rh = = (6.5-25)
área del empaque volumen de lecho a
Por otra parte, es de utilidad definir el área específica del lecho en términos del área
específica Sp de las partículas de empaque:
Ap
a = S p (1 − ε ) = (1 − ε ) (6.5-26)
Vp
donde Ap es el área de la partícula de empaque y Vp es su volumen. Por ejemplo, para una

partícula de empaque esférica de diámetro Dp tenemos que:
π D p2 6
Sp = = (6.5-27)
π D3p 6 Dp
Si sustituimos esta ecuación en la (6.5-26) obtenemos el área específica de un lecho de

partículas esféricas de diámetro Dp es:
6 (1 − ε )
a= (6.5-28)
Dp
Esta relación es válida para lechos fijos y fluidizados. Cuando las partículas que forman el
lecho no son esféricas, podemos definir un diámetro equivalente como:
6 6V p
D p ,eq = = (6.5-29)
Sp Ap
R. LOBO/2007
y usar la Ec.(6.5-28) para calcular el valor correspondiente de a.
Es necesario advertir que el área por unidad de volumen de un lecho empacado que hemos
definido aquí, a, no debe ser confundida con el área interfacial por unidad de volumen de
lecho entre dos fluidos inmiscibles que fluyen a través de un lecho empacado. La primera es
aproximadamente el área por volumen de lecho del empaque sólido y es independiente del
flujo del fluido, en tanto que la segunda es el área de contacto fluido-fluido en un lecho
empacado y es función, entre otras cosas, de la magnitud de los flujos de las fases fluidas.
La velocidad promedio del fluido en lechos empacados y fluidizados puede definirse de dos
maneras. La más común es la velocidad superficial, vb,sup, que se define como el flujo
volumétrico QV (m3/s, por ejemplo) del fluido dividido por el área seccional del tubo At que
contiene al lecho (como si el tubo estuviera vacío):
QV w ρ
vb,sup = = (6.5-30)
At At
donde w es el flujo másico y ρ es la densidad del fluido. Aunque la velocidad superficial

representa un valor menor que la velocidad promedio real, su uso es muy común debido a
la facilidad con que puede calcularse. La velocidad promedio real del fluido en un lecho
empacado o fluidizado se denomina velocidad intersticial, vb,int, y se define considerando
que la fracción del área seccional del lecho que está compuesta de huecos es igual a la
fracción vacía ε. Entonces,
vb,sup Qv At w ρ At
vb,int = = = (6.5-31)
ε ε ε
Una vez establecidos estos conceptos en lechos empacados, podemos pasar a enunciar las
correlaciones.
Para gases y líquidos, Sherwood y col. recomiendan estimar el coeficiente de transferencia

de masa promedio mediante la siguiente correlación
kc
gases y líquidos: jD = (Sc )2 3 = 1.17 Re−D0.415
vb,sup p
(6.5-32)
10 < Re D p < 2,500 ; 0.6 < Sc < 2, 000
En esta correlación el número de Reynolds se define como:
ρ vb,sup D p
Re D p = (6.5-33)
µ
R. LOBO/2007
donde Dp es el diámetro de una esfera que posee la misma área por unidad de volumen que
la partícula de forma el empaque.
Existen otras correlaciones aplicables sólo a gases o sólo a líquidos con rangos de validez
diversa. Para gases fluyendo en un lecho empacado de esferas, puede usarse la correlación
de Gupta y Thodos1
kc
gases: ε jD = ε ( Sc )2 3 = 2.06 Re−D0.575
vb,sup e
(6.5-34)
90 < Re De < 4, 000
Para líquidos fluyendo a baja velocidad en un lecho empacado de esferas puede usarse la
correlación de Wilson y Geankoplis2
kc
líquidos: ε jD = ε ( Sc )2 3 = 1.09Re−D2e 3
vb,sup (6.5-35)
0.002 < Re De < 55 ; 165 < Sc < 70, 600
Un lecho fluidizado es esencialmente una columna donde se encuentran partículas sólidas

pequeñas en flujo, suspendidas en un flujo de gas o de líquido. Para lechos fluidizados de
partículas esféricas en gases o líquidos, Gupta y Thodos3 desarrollaron la siguiente
correlación:
kc 0.863
gases y líquidos: ε jD = ε ( Sc )2 3 = 0.010 +
vb,sup Re0.58
De − 0.483 (6.5-36)
20 < Re De < 3, 000
En las Ecs.(6.5-34), (6.5-35) y (6.5-36), el número de Reynolds se calcula con base en el

diámetro de la esfera De y la velocidad del fluido es la superficial.
G. Equipos con contacto fluido-fluido, fluido-sólido y fluido-fluido-sólido
Existe una gran diversidad de operaciones industriales donde ocurre la transferencia de

masa gas-líquido, gas-sólido, líquido-líquido y gas-líquido-sólido. Como ejemplos
1 Gupta, AS y Thodos, G. (1963) AIChE J., 9, 751., y___ (1964), IEC Fund.. 3, 218.
2 Wilson, EJ y Geankoplis CJ (1966) IEC Fund., 5, 9.
3 Gupta, AS y Thodos, G. (1962) AIChE J., 8, 609.
R. LOBO/2007
podemos citar absorbedores, columnas de destilación, lechos fluidizados, reactores gas-

líquido con y sin agitación mecánica, reactores gas-líquido-sólido suspendido, reactores
gas-líquido sólido fijo, etc. En estos equipos, además de la problemática para estimar el
coeficiente de transferencia de masa, existe la necesidad de estimar adecuadamente el área
interfacial de contacto entre las fases. Para este tipo de operaciones se han desarrollado
correlaciones específicas que pueden consultarse en algunos textos y en artículos de
revisión especializados que se encuentran en la literatura1.
_______________________
No está de más que advirtamos, en general, que la velocidad del fluido en las correlaciones
no siempre aparece en ellas como tal, sino que se presenta en términos de flujos másicos
(w) o flujos volumétricos (Qv) del fluido. Esto también ocurre en las definiciones del
número de Reynolds. La razón práctica para ello es que la velocidad de un fluido
generalmente no se mide, sino que se calcula a partir de variables que se pueden medir con
la instrumentación disponible. Ello no debe confundir al estudiante, quien sólo debe
recordar que:
flujo volumétrico flujo másico densidad

velocidad del fluido = = (6.5-37)
área transversal al flujo área transversal al flujo
Ejemplo 6.5-1. Cálculo del coeficiente de transferencia de masa en la

evaporación de agua líquida desde una superficie plana
En un experimento a 25 °C y 1 atm, una corriente de aire seco (B) fluye sobre una
superficie sólida plana de 3 m de longitud la cual se mantiene siempre saturada de agua
líquida. Al contacto con la corriente de aire, el agua líquida se evapora y se transfiere al
aire. A estas condiciones, la viscosidad del aire es 0.0184 cp y la difusividad del vapor de
agua (A) en el aire es 0.26 cm2/s. La velocidad del aire es de 2.0 m/s
Calcula:
(a) El coeficiente de transferencia de masa local para la evaporación del agua a los 50 cm
de longitud sobre la superficie de agua y el coeficiente de transferencia de masa promedio
sobre esta distancia;
(b) El valor del flux de agua que se evapora a los 50 cm y el flux promedio sobre esta
distancia;
1 Ver por ejemplo: Kastanek, F, Zahradnik, J, Kratochvil, J y Cermak, J (1993) Chemical Reactor for Gas-
Liquid Systems, Ellis Horwood. Londres. Bailey, JE y Ollis DF (1986) Biochemical Engineering
Fundamentals. McGraw Hill. New York. Cheremisinof, NP (Ed) (1986) Handbook of Heat and Mass
Transfer. Vol 2: Mass Transfer and Reactor Design. Gulf Publishing. Houston. Seader JD y Henley EJ (1998)
Separation Process Principles. John Wiley. New York. Treybal, RE (1980) Mass Transfer Operations. 3a.
Edición. McGraw Hill. New York. Ramachandran, PA y Chaudhari, RV (1983) Three Phase Catalytic
Reactors. Gordon and Breach. Londres.
R. LOBO/2007
(c) La distancia sobre la longitud de la superficie a la cual el flujo se torna turbulento;

(d) El coeficiente de transferencia de masa promedio sobre toda la longitud de la placa;
(e) el coeficiente de transferencia de masa promedio sobre toda la longitud de la superficie
de agua, si la velocidad del aire es de 10 m/s.
Solución. El estudio de este ejemplo es importante porque en él ilustraremos las

definiciones de un coeficiente de transferencia de masa local y uno promedio, y veremos
cuál es la correlación adecuada dependiendo de si el flujo es laminar o turbulento.
(a) Comencemos por determinar las propiedades del fluido involucrado en la transferencia
de masa, para determinar el régimen de flujo y seleccionar la correlación adecuada. En la
vecindad de la interfase, la mezcla gaseosa que fluye sobre la superficie de agua líquida es
principalmente aire y un poco de vapor de agua. La presión de vapor del agua a 25 °C es
0.0313 atm, por lo que podemos considerar que la mezcla vapor de agua-aire (A-B) es
diluida y que sus propiedades pueden ser calculadas como si fueran las del aire (B):
ρ mezcla ≅ ρ B =
MB P
=
( 28.8 g mol )(1 atm ) −3
= 1.18 × 10 g cm
3
RT (82.05 cm atm mol K ) ( 298 K )

3
µmezcla ≅ µ B = 0.0184 cp=1.84 × 10 −4 g cm ⋅ s
µ B 1.84 × 10−4 g cm ⋅ s
ν mezcla ≅ ν B = = = 0.156 cm s
2
ρB 1.18 × 10−3 g cm3
Entonces, los números de Reynolds y Schmidt son
zvb, z (50 cm)(200 cm s)

Re z =50 cm = = = 64,102 < 2 × 10
5
⇒ flujo laminar
νB 0.156 cm s
2
νB 2
0.156 cm s
Sc = = 2
= 0.6
DAB 0.26 cm s
La definición del problema nos pregunta por el coeficiente de transferencia de masa a los
50 cm de longitud de la superficie de agua líquida. Ello quiere decir que debemos calcular
el coeficiente de transferencia de masa local, específicamente a esa distancia. La
correlación para este coeficiente local para flujo laminar la proporciona la Ec.(6.5-1):
kc z
Sh z = = 0.332 Re0.5
z Sc
13
DAB
de donde podemos despejar el coeficiente de transferencia de masa local:
R. LOBO/2007
kc = 0.332 ( DAB z ) Re0.5

z Sc
13
(6.5-38)
Si sustituimos los valores de las variables, obtenemos el valor del coeficiente de

transferencia de masa local a los 50 cm de longitud y cuando la velocidad del fluido es 2.0
m/s:
= 0.332 ( 0.26 cm 2 /s 50 cm ) ( 64,102) ( 0.6)

0.5 13
kc 50 cm = 0.37 cm s
Para calcular el coeficiente promedio, nuevamente debemos determinar el régimen de flujo

mediante el número de Reynolds, ahora definido por
vb, z Lz (200 cm s)(50 cm)

Re Lz = = = 64,102 < 2 × 10
5
⇒ flujo laminar
ν 2
0.156 cm s
Entonces, el coeficiente de transferencia de masa promedio en los primeros 50 cm de la

superficie de agua lo obtenemos de la Ec.(6.5-2):
kc Lz
Sh Lz = = 0.664Re1Lz2 Sc1 3
DAB
de donde podemos obtener
kc = 0.664 ( DAB Lz ) Re0.5

Lz Sc
13
(6.5-39)
Sustituyendo valores en esta ecuación obtenemos:
( kc )0 = 0.664 ( 0.26 cm 2 /s 50 cm ) ( 64,102)

50 cm
( 0.6)
0.5 13
= 0.74 cm s
Aquí, por claridad, hemos introducido el rango de valores de la distancia sobre la cual se
calculó este coeficiente promedio.
(b) El flux local a los 50 cm y el flux promedio sobre esta distancia del agua que se evapora
están dados, respectivamente, por:
(
N A i = k c CA i − CA b ; ) (
N A i = kc CA i − CAb ) (6.5-40)
Los valores de los coeficientes de transferencia de masa ya los calculamos arriba. Nos resta
encontrar el valor de las concentraciones. Dado que el aire es seco, la concentración de
agua en el seno de la corriente gaseosa CAb podemos considerarla igual a cero. La
R. LOBO/2007
concentración de vapor en la interfase aire-agua del lado del aire CAi la podemos calcular
mediante la hipótesis de equilibrio en la interfase:
CA i = CAsat =
pAsat
=
( 0.0313 atm ) −6
= 1.28 × 10 mol/cm
3
RT (82.05 cm atm mol K ) ( 298 K )

3
La sustitución de los valores correspondientes en las Ecs.(6.5-40) resulta en:
NAi 50 cm ( )
= ( 0.37 cm/s ) 1.28 × 10 −6 mol/cm3 − 0 = 4.74 × 10 −7 mol/cm 2 s
= 17.1 mol/m 2 h
( NA, z )0 ( )
50 cm
= ( 0.74 cm/s ) 1.28 × 10 −6 mol/cm3 − 0 = 9.48 × 10 −7 mol/cm 2s
= 34.2 mol/m 2 h
Estos fluxes son relativamente bajos porque el flujo es laminar y la concentración

interfacial del vapor de agua es pequeña.
(c) El cambio de flujo laminar a turbulento no es abrupto, sino gradual. Sin embargo, aquí
consideraremos que el cambio de régimen se da cuando el número de Reynolds es 2×105.
Entonces, podemos definir una distancia crítica zc a la cual el flujo cambia de laminar a
turbulento de la siguiente manera:
zc vb, z
Re zc = 2 × 105 =
ν
De esta relación podemos obtener zc:
2 ×105 ν ( 2 × 10 )( 0.156 cm /s )
5 2
zc = = = 156 cm
vb, z ( 200 cm/s )
Lo anterior quiere decir que la longitud total de la superficie de 300 cm está dividida en una
zona de flujo laminar y una zona de flujo turbulento:
0 < zlam ≤ 156 cm ; 156 < zturb ≤ 300 cm
(d) En vista de la existencia de zonas de flujo laminar y de flujo turbulento, para encontrar
el coeficiente de transferencia de masa promedio sobre la longitud total Lz=300 cm de la
superficie, debemos promediar los coeficientes locales sobre ella. El coeficiente local para
la zona laminar está dado por la Ec.(6.5-38) y el coeficiente local para la zona turbulenta lo
podemos obtener de la Ec.(6.5-4):
R. LOBO/2007
kc z
Sh z = = 0.0296 Re0.8
z Sc
13
DAB
de donde
⎛D ⎞
kc = 0.0296 ⎜ AB ⎟ Re0.8
z Sc
13
(6.5-41)
⎝ z ⎠
Entonces, el coeficiente promedio sobre la longitud total de la superficie es.
1 ⎛
( kc )lam dz + ∫z ( kc )turb dz ⎞⎟
zc Lz
kc = ⎜
Lz ⎝ ∫0 c ⎠
(6.5-42)
Si sustituimos las Ecs.(6.5-38) y (6.5-41) en la (6.5-42):
DAB Sc1 3 ⎡ dz ⎤
0.5 0.8
⎛ vb, z ⎞ zc dz ⎛v ⎞ Lz
kc =
Lz
⎢0.332 ⎜
⎢⎣ ⎝ ν ⎠
⎟ ∫0 z 0.5
+ 0.0296 ⎜ b, z ⎟
⎝ ν ⎠
∫
zc
⎥
z 0.2 ⎥
⎦
(6.5-43)
DAB Sc1 3 ⎡
=
Lz ⎣
0.664 Re0.5 0.8 0.8 ⎤
zc + 0.037 Re Lz − Re zc ⎦ ( )
Evidentemente, podemos escribir esta ecuación en términos del número de Sherwood
promedio:
Sh Lz =
kc Lz ⎡
DAB
= 0.664 Re0.5
⎣
0.8 0.8 ⎤
zc + 0.037 Re Lz − Re zc ⎦ Sc
13
( ) (6.5-44)
Para calcular el coeficiente promedio sobre toda la placa sólo nos resta conocer ReLz , ya
que el resto de las cantidades en la Ec.(6.5-43) son conocidas:
Lz vb, z (300 cm)(200 cm/s)

Re Lz = = = 3.84 × 10
5
ν (0.156 cm /s)
2
Finalmente, podemos obtener el valor del coeficiente promedio sobre toda la superficie:
2 13
⎡0.664 (2 × 105 )0.5 + 0.037 ( (3.84 × 105 )0.8 − (2 × 105 )0.8 )⎤
(0.26cm /s)(0.6)
kc = ⎣ ⎦
300 cm
= 0.32 cm/s
(e) Para el caso en que la velocidad del aire cambia a 10 m/s, nuestra intuición nos dice que
la longitud de la zona laminar será menor y que la turbulenta aumentará. Esto lo podemos
R. LOBO/2007
verificar calculando la distancia crítica zc para el cambio del régimen de flujo de laminar a
turbulento, de manera semejante a como lo hicimos anteriormente:
2 ×105 ν ( 2 × 10 )( 0.156 cm /s )
5 2
zc = = = 31.2 cm
vb, z (1000 cm/s )
Esto significa que la longitud total de la superficie de 300 cm está ahora dividida en una
zona de flujo laminar y una zona de flujo turbulento como sigue:
0 < zlam ≤ 31.0 cm ; 31.2 < zturb ≤ 300 cm
Para la velocidad de 10 m/s, el número de Reynolds sobre toda la superficie es
Lz vb, z (300 cm)(1000 cm/s)

Re Lz = = = 1.92 × 10
6
ν 2
(0.156 cm /s)
Con fines ilustrativos, calcularemos el coeficiente promedio de manera directa a partir de la

Ec.(6.5-5)
kL
Sh Lz = c z = 0.0365Re0.8
Lz Sc
13
DAB
y también lo haremos mediante el método de promediar los coeficientes locales que

acabamos de exponer en el inciso anterior. De la última ecuación conseguimos,
⎛D ⎞
kc = 0.0365 ⎜ AB ⎟ Re0.8
Lz Sc
13
(6.5-45)
⎝ Lz ⎠
donde, al sustituir los valores correspondientes obtenemos:
⎛ 0.26 cm/s ⎞
⎟ (1.92 × 10 ) ( 0.6 ) = 2.84 cm/s
6 0.8 13
kc = 0.0365 ⎜ (6.5-46)
⎝ 300 cm ⎠
Para obtener el valor del coeficiente promedio por vía del promediado de los coeficientes
locales usamos la Ec.(6.5-43), de donde podemos calcular
2 13
⎡0.664 (2 × 105 ) 0.5 + 0.037 ( (1.92 × 106 )0.8 − (2 × 105 ) 0.8 )⎤
(0.26cm /s)(0.6)
kc = ⎣ ⎦
300 cm
(6.5-47)
= 2.62 cm/s
R. LOBO/2007
Estos resultados merecen algunos comentarios. Al aumentar la velocidad del fluido a 10

m/s la zona laminar se redujo considerablemente hasta alrededor del 10 % de la longitud
total de la superficie y la zona turbulenta aumentó de manera correspondiente. Este cambio
se ve reflejado también en la magnitud del coeficiente de transferencia de masa promedio.
Éste aumentó prácticamente un orden de magnitud, respecto del que obtuvimos para la
velocidad de 2.0 m/s. El valor del coeficiente promedio calculado de la Ec.(6.5-47) es
ligeramente menor que el calculado de la Ec.(6.5-46) porque esta última ecuación supone
que toda la superficie está expuesta a flujo turbulento, en tanto que aquélla toma en cuenta
la pequeña zona laminar. Sin embargo, para fines prácticos, ambos valores pueden
considerarse iguales.
Ejemplo 6.5-2. Absorción de CO2 en una columna de pared mojada
Una corriente de agua desciende en forma de película por el interior de un tubo de 10 cm de

diámetro y 1.0 m de longitud, en flujo laminar y a una velocidad promedio de 10 cm/s. El
espesor de la película es de 0.025 cm. Una corriente de CO2 puro asciende lentamente por
el interior del tubo y se absorbe en el agua a 25 ºC y 1 atm. En estas condiciones, la
solubilidad de CO2 es 3.36×10−5 mol/cm3 y su difusividad en agua es 1.96×10−5 cm2/s. La
columna opera en estado estacionario.
Calcula la concentración del CO2 disuelto en el agua a la salida de la columna:

(a) usando un coeficiente de transferencia de masa local;
(b) usando un coeficiente de transferencia de masa promedio
Solución. En la Figura 6.5-1 se muestra una representación esquemática del sistema.

Consideraremos que y es la dirección del flujo de líquido, Lz es el espesor de la película y
Dt es el diámetro de la columna. Debido a que el diámetro de la columna es mucho mayor
que el espesor de la película líquida, puede despreciarse la curvatura de ésta. El CO2 en la
fase gaseosa es puro y por ello no existen gradientes de concentración en ella.
(a) Comenzamos realizando un balance de masa sobre el CO2 (A) en un elemento

diferencial de volumen de la película líquida πDtLz∆y. El área perpendicular al flujo de
líquido es πDtLz y el elemento de área interfacial de transferencia de masa contenida en el
elemento de volumen es ∆Ai=πDt∆z. Este balance se expresa de la siguiente manera:
( que entra en y )( que sale en y + ∆y )

flujo de A disuelto − flujo de A disuelto + ⎛ flujo de A absorbido ⎞ = 0
⎜ a través de ∆A ⎟
⎝ i ⎠
R. LOBO/2007
agua
y CA b , vb, y
N Ai
Ly ∆y
∆y
N Ai
Lz
Dt
CO 2
Figura 6.5-1. Representación de la absorción de CO2 en una columna de pared mojada.
En términos matemáticos este balance es:
(
π Dt Lz N A, y y − N A, y y +∆y )+ N A i ∆Ai =0 (6.5-48)
donde NA,y es el flux total de A en el líquido, πDtLz es el área perpendicular al flujo de

líquido y NAi es el flux de absorción de A en el área ∆Ai=πDt∆y. Como el flujo es laminar,
el flux difusivo en dirección y es insignificante, por lo que
N A, y = J A, y + CA b vb , y ≅ CA b vb , y (6.5-49)
donde vb,y es la velocidad promedio del líquido y es constante, y CAb es la concentración

promedio del soluto en el líquido y varía con y. El flux interfacial está dado por:
(
N A i = k c CA i − CA b ) (6.5-50)
La sustitución de las Ecs.(6.5-49) y (6.5-50) en la (6.5-48) produce:
(π Dt Lz vb, y ) ( CAb y −CA b y +∆y ) + k (C

c Ai )
− CA b π Dt ∆y = 0 (6.5-51)
R. LOBO/2007
Si dividimos por el elemento diferencial de volumen, tomamos lím ∆y→0 y rearreglamos,

llegamos a:
dCA b kc
−
dy
+ (
Lz vb, y
)
CA i − CA b = 0 (6.5-52)
En esta ecuación CAi es la concentración en la interfase líquido/gas, es constante y se

corresponde con la de saturación porque el gas es CO2 puro; la concentración promedio CAb
en el líquido irá aumentando con la longitud de la columna medida desde su parte superior.
Debido a que el agua entra a la columna libre de CO2, la condición de frontera es:
y=0; CA b = 0 (6.5-53)
La Ec.(6.5-52) aún no puede integrarse. Debido a que seleccionamos una diferencia de

concentración local y un área de transferencia local para hacer el balance de masa, el
coeficiente de transferencia de masa que debemos usar es también uno local. La correlación
adecuada para el cálculo de un coeficiente local en películas descendentes laminares es la
Ec.(6.2-19)
kc y
( ) ( )
12 12
Sh y = = 0.691 Re y Sc = 0.691 Pem y (6.5-54)
DAB
de donde podemos obtener que
12 12
⎛ D ⎞ ⎛ y vb, y ⎞ ⎛ vb, y DAB ⎞
kc = 0.691⎜ AB ⎟ ⎜ ⎟ = 0.691⎜ ⎟ (6.5-55)
⎝ y ⎠ ⎝ DAB ⎠ ⎝ y ⎠
Si sustituimos la Ec.(6.5-55) en la (6.5-52), separamos variables e integramos la ecuación

resultante, después de rearreglar obtenemos:
CA b ⎡ ⎛D y⎞ ⎤
12
= 1 − exp ⎢ −1.38 ⎜ 2AB ⎟ ⎥ (6.5-56)

CA i ⎢ ⎜ Lz vb, y ⎟ ⎥
⎣ ⎝ ⎠ ⎦
Esta ecuación proporciona la variación de la concentración promedio del CO2 a lo largo de

la columna de pared mojada. Sólo nos resta evaluar esta concentración a la distancia
y=Ly=100 cm, para obtener la concentración promedio del CO2 a la salida de la columna:
R. LOBO/2007
⎧⎪ ⎡ ⎛ (1.96 × 10−5 cm 2 /s )(100 cm) ⎞ ⎤ ⎫⎪

12
CA b y =100 cm = ( 3.36 × 10 mol/cm ) ⎨1 − exp ⎢−1.38 ⎜

−5 3
2 ⎟ ⎥ ⎬ (6.5-57)
⎪⎩ ⎢
⎣ ⎝ (0.025 cm) (10 cm/s) ⎠ ⎥⎦ ⎭⎪
−5
= 1.809 × 10 mol/cm 3
(b) Este cálculo puede hacerse de otra manera, a partir de un balance global de masa sobre
el soluto en toda la columna, entre la entrada y la salida de ella. La cantidad de soluto total
entre entrada y salida de la columna debe ser igual a la cantidad total transferida. La
cantidad total transferida se expresa en términos del flux interfacial promedio en toda la
columna multiplicada por el área total de transferencia de la columna. Este enfoque ya lo
utilizamos en el Ejemplo 6.3-1. Entonces, el balance de masa sobre el soluto en toda la
columna se expresa por:
( ) ( )
QV CA b Ly − QV CA b ( 0 ) = π Dt Lz vb, y CA b Ly = N A i Ai

( ) (6.5-58)
=0
donde QV es el flujo volumétrico del fluido y Ai=πDtLy es el área interfacial total contenida
en la columna. En el Ejemplo 6.3-1 demostramos que el flux interfacial promedio se
expresa en términos del coeficiente de transferencia de masa promedio y la diferencia de
concentración media logarítmica:
N A i = kc ( ∆CA )ln = kc
( )
[CA i − CA b Ly ] − [CA i − CA b ( 0 )]
(6.5-59)
⎡ CA i − CA b L y
ln ⎢
( ) ⎤⎥
⎢⎣ CA i − CA b ( 0 ) ⎥⎦
Para nuestro caso, en vista de que CAb(0)=0, esta expresión se reduce a:
N A i = kc
( )
−C A b L y
(6.5-60)
⎛ CA b ( L y ) ⎞
ln ⎜ 1 − ⎟
⎜ CA i ⎟
⎝ ⎠
Si sustituimos Ai y la Ec.(6.5-60) en la ecuación del balance de masa global (6.5-58),

después de simplificar y rearreglar llegamos a:
⎡ ⎛ kc L y ⎞⎤
( )
CA b Ly = CA i ⎢1 − exp ⎜ −
⎜ vb, y Lz ⎟⎟ ⎥ (6.5-61)
⎢⎣ ⎝ ⎠ ⎥⎦
Para calcular la concentración promedio de CO2 a la salida de la columna, CAb(Ly),

únicamente nos resta calcular el coeficiente de transferencia de masa promedio kc , ya que
R. LOBO/2007
el resto de las cantidades son conocidas. El coeficiente promedio lo podemos obtener de la

correlación aplicable a películas líquidas laminares dada por la Ec.(6.5-6):
kc Lz
( ) (L )
12 12
Sh Lz = = 0.691 Pem Lz z Ly
DAB
de donde podemos obtener que:
12 12
⎛ D ⎞ ⎛ 4vb, y Lz ⎞ ⎛ Lz ⎞
kc = 0.691⎜ AB ⎟ ⎜ ⎟ ⎜⎜ ⎟⎟ (6.5-62)
⎝ Lz ⎠ ⎝ DAB ⎠ ⎝ Ly ⎠
Podemos sustituir los valores correspondientes para conseguir:
12 12
⎛ 1.96 × 10 −5 cm 2 /s ⎞ ⎛ 4(10 cm/s)(0.025 cm) ⎞ ⎛ 0.025 cm ⎞
kc = 0.691 ⎜ ⎟⎜ −5 2 ⎟ ⎜ ⎟
⎝ 0.025 cm ⎠ ⎝ 1.96 × 10 cm /s ⎠ ⎝ 100 cm ⎠
= 1.93 × 10 −3 cm/s
Estamos ahora ya en posición de calcular la concentración del CO2 a la salida de la

columna de pared mojada, sustituyendo los valores correspondientes en la Ec.(6.5-61).
⎡ ⎛ (1.93 × 10−3 cm/s)(100 cm) ⎞⎤

( )
CA b L y −5
= (3.36 × 10 mol/cm ) ⎢1 − exp ⎜ −
3
⎟⎥
⎣ ⎝ (10 cm/s)(0.025 cm) ⎠⎦
= 1.807 × 10−5 mol/cm3
Si comparamos este resultado con el obtenido en la Ec,(6.5-57) podemos ver que son
prácticamente idénticos, lo cual no debe sorprendernos. Los cálculos del inciso (a) fueron
realizados con base en variables locales y los del inciso (b) fueron hechos con base en
variables promedio. Era necesario demostrar concretamente que los cálculos son
equivalentes, con tal que cualquier método de realizarlos se haga correctamente. No
obstante, seguimos adoptando la posición de que el método que usa variables locales tiene
la ventaja de hacer reflexionar al estudiante sobre la verdadera naturaleza local de los
fenómenos de transferencia de masa.
6.6 Analogías entre las transferencias de momentum, calor y masa

Si revisamos en la Sec. 1.6 las ecuaciones que describen el transporte molecular de
momentum, calor y masa (leyes de viscosidad de Newton, de conducción de calor de
Fourier y de difusión de masa de Fick) notaremos que son semejantes en su forma
R. LOBO/2007
matemática. También, las ecuaciones generales de cambio adimensionales para cada tipo
de transporte obtenidas en el Capítulo 5 muestran grandes analogías en su forma. Ello
obedece a que los mecanismos físicos de transporte difusivo de momentum, calor y masa
son semejantes, y los mecanismos de transporte convectivo de momentum, calor y masa
también los son entre sí. Esta semejanza física en los mecanismos de transporte es de gran
utilidad, ya que con alguna frecuencia permite estudiar uno de los fenómenos de transporte
estudiando otro.
Así, hemos visto en el Capítulo 5 el surgimiento de grupos adimensionales análogos. Por

ejemplo, recordamos que el número de Prandtl (Pr) en transferencia de calor es análogo al
número de Schmidt (Sc) en transferencia de masa, ya que tienen significados físicos
análogos. También hemos visto que en transferencia de calor en número adimensional que
incluye al coeficiente de transferencia de calor (h) es el número de Nusselt (Nu) y que su
análogo en transferencia de masa es el número de Sherwood (Sh), que incluye al
coeficiente de transferencia de masa. La obtención de coeficientes de transferencia de masa
a través de las analogías se basa en la hipótesis de que las transferencias de momentum,
calor y masa ocurren por medio de mecanismos y velocidades semejantes. Reynolds, fue
quien, en 1874, dio inicio al uso de las analogías, en un problema de transferencia de
momentum y calor.
Las analogías son de utilidad para condiciones en que el flujo alrededor de la interfase es
turbulento. En esta sección trataremos dos analogías: la de Reynolds, que es importante por
su sencillez y contenido conceptual; y la de Chilton-Colburn, que es de gran utilidad
práctica. Para un análisis más detallado de las analogías, puede consultarse a Sherwood et
al (1975) y a Welty el al (1976).
A. La analogía de Reynolds
En la Sección 4.2 introducimos el concepto de difusividad turbulenta para tomar en cuenta el

transporte de masa por convección forzada turbulenta en la vecindad de una interfase.
Para lo que sigue, supondremos que las transferencias se dan en dirección y, y que la interfase
desde la cual se originan está en y=0. La expresión del flux de masa para el caso de una
solución diluida es (ver Sec. 5.6)
( ()
N A, y = − DAB + DAB
t
) dCd y
A
(6.6-1)
donde ED es la difusividad turbulenta promedio de masa. Se pueden definir también

expresiones análogas para los fluxes de momentum y calor:
τ (ytotal
z
)
(
= − ν +ν ( )
t
) d (dρyv )
z
(6.6-2)
R. LOBO/2007
(
d ρ Cˆ p ,m T )
qy = − α + α ( (t )
) dy
(6.6-3)
donde ν y α son la viscosidad cinemática y la difusividad térmica, respectivamente;

ν ( t ) y α ( t ) son las difusividades turbulentas promedio de momentum y calor. Si definimos
las siguientes variables adimensionales
vz Ti − T C A i − CA
vz∗ = ; T∗ = ; CA∗ = (6.6-4)
vb, z Ti − Tb CA i − CA b
las sustituimos en las Ecs.(6.6-1) a (6.6-2), las evaluamos en la interfase (y=0), e igualamos
a la definición del flux correspondiente en la interfase, obtenemos
( (t )
N A i = − DAB + DAB CA i − CA b )( )
⎛ d CA∗ ⎞
⎜− (
⎟ = k c CA i − CA b
⎝ d y ⎠ y =0
) (6.6-5)
( t
) ⎛ dT ∗ ⎞
qi = − α + α ( ) ρ Cˆ p ,m (Ti − Tb ) ⎜ − ⎟ = h (Ti − Tb )
⎝ d y ⎠ y =0
(6.6-6)
(
τ i = − ν +ν ( t ) ρ vb, z ⎜ ) ⎛ d vz∗ ⎞
⎝ d y ⎠ y =0
1
⎟ = 2 ρ vb, z f
2
(6.6-7)
donde qi y τi son, respectivamente, el flux de calor y de momentum en la interfase; kc es el

coeficiente de transferencia de masa, h es el coeficiente de transferencia de calor, f es el
factor de fricción y Cˆ p ,m es la capacidad calorífica específica media entre las temperaturas
Ti y Tb. La velocidad, la concentración y la temperatura que se usan en estas ecuaciones son
valores promedio en el tiempo y ello se indica con la barra encima de ellas (ver sec. 7.6
para una discusión más amplia sobre el significado de estos valores promedio en el tiempo.
La analogía de Reynolds establece la igualdad de mecanismos y de rapidez en las

transferencias de masa, calor y momentum. Ello implica que los gradientes adimensionales
de concentración, temperatura y velocidad sean iguales. Entonces, si de las Ecs.(6.6-5) a
(6.6-7) despejamos los gradientes y los igualamos, con un poco de álgebra llegaremos a
kc h 1 f vb, z
= = (6.6-8)
( (t )
DAB + DAB ) (
ρ Cˆ p ,m α + α (t )
) ( 2 ν +ν (t ) )
R. LOBO/2007
Ahora bien, en condiciones de flujo turbulento en la vecindad de la interfase, las

contribuciones difusivas son pequeñas, por lo que pueden considerarse insignificantes
comparadas con las difusividades turbulentas. La teoría de la difusividad turbulenta y los
datos experimentales con que se cuenta a la fecha tienden a indicar que las difusividades
turbulentas son semejantes en magnitud, por lo que es razonable suponer que:
()
= α ( ) =ν ( )
t t t
DAB (6.6-9)
entonces, si hacemos uso de lo anteriormente dicho, podemos obtener la forma final de la

analogía de Reynolds
kc h f
= = (6.6-10)
vb, z ρ Cˆ p ,m vb, z 2
Esta analogía resulta de interés por su sencillez y nos permite estimar coeficientes de
transferencia de masa a partir de coeficientes de transferencia de calor o factores de fricción
de casos análogos. Por ejemplo, si contáramos con un valor del coeficiente de transferencia
de calor para el enfriamiento de una esfera en un fluido que fluye sobre ella en régimen
turbulento, sería posible calcular el coeficiente de transferencia de masa para el caso en que
un soluto se estuviera disolviendo desde una esfera sobre la cual pasa un fluido en flujo
turbulento. Mejor aún, si se cuenta con una correlación para el coeficiente de transferencia
de calor o el factor de fricción para una determinada situación, puede obtenerse una
correlación para el coeficiente de transferencia de masa para la situación análoga.
Cuando se compara la analogía de Reynolds con datos experimentales se encuentra que

funciona bien para gases, para los cuales Sc≈Pr≈1. Sin embargo, la analogía falla
totalmente para líquidos, donde los números adimensionales mencionados difieren en
forma considerable. La conclusión de comparar la analogía de Reynolds con datos
experimentales es que los mecanismos y la rapidez de la transferencia turbulenta de
momentum, calor y masa sí son semejantes en gases, pero no en líquidos.
B. La analogía de Chilton-Colburn
La idea de las analogías es interesante y de consecuencias prácticas. Por ello, se ha hecho

un considerable esfuerzo por desarrollar otras que superen las limitaciones de la de
Reynolds y a la fecha existe una abundante cantidad de analogías publicadas en la literatura
especializada. Sin embargo, la extensión más útil y comúnmente usada es la atribuida a
Chilton y Colburn1. Estos autores, a partir del análisis de datos experimentales, encontraron
que si en la analogía de Reynolds multiplicaban el grupo que contiene el coeficiente de
transferencia de masa por el número de Schmidt elevado a la dos tercios; y si multiplicaban
1 Chilton, TH y Colburn, AP (1934) Ind. Eng. Chem., 26, 1183.
R. LOBO/2007
el grupo que contiene al coeficiente de transferencia de calor por el número de Prandtl

elevado también a la dos tercios, la analogía mejoraba bastante y se podían estimar
razonablemente bien los datos de líquidos. Entonces, la analogía de Chilton-Colburn es
kc h f
(Sc )2 3 = jD = jH = ˆ ( Pr )2 3 = (6.6-11)
vb, z ρ C p ,m vb, z 2
Esta es la analogía más utilizada para estimar coeficientes de transferencia de masa,

principalmente a partir de datos de problemas análogos de transferencia de calor. La
igualdad con f/2 es sólo aproximada, particularmente en situaciones en que el flujo ocurre
alrededor de superficies curvas, las desviaciones pueden ser grandes. El término que
contiene al coeficiente de transferencia de masa con frecuencia se denomina factor jD, y el
que contiene el coeficiente de transferencia de calor, factor jH.
Por último, hay que reiterar que las propiedades del fluido deben evaluarse a las
condiciones promedio entre las que prevalecen en la interfase y en el seno del fluido. Esto
es particularmente importante en cuanto el efecto de la temperatura sobre las propiedades y
éstas deben calcularse a una temperatura promedio entre la que prevalece en la interfase y
la del seno del fluido.
Ejemplo 6.6-1. Estimación del coeficiente de transferencia de masa mediante la

analogía de Chilton-Colburn
Consideremos la evaporación de amoníaco líquido hacia una corriente de aire a 27 ºC y 1

atm. En un experimento bajo las mismas condiciones, se determinó que el coeficiente de
transferencia de calor es 0.108 cal/cm2⋅s⋅ºK. La mezcla gaseosa que forman el aire y el
amoníaco es diluida. A las condiciones del experimento, las propiedades del aire son:
µ=0.0184 cp; Cˆ p,m =0.24 cal/g⋅ºK; kT=6.25×10−5 cal/s⋅cm⋅ºK y el coeficiente de difusión del
amoníaco en aire es 0.23 cm2/s.
Estima el coeficiente de transferencia de masa para la evaporación del amoníaco.
Solución. En este problema podemos aplicar directamente la analogía de Chilton-Colburn

dada por la Ec.(6.6-11)
k h
( Sc )2 3 = ˆ ( Pr )2 3
vb, z ρ C p ,m vb, z
de donde podemos obtener
R. LOBO/2007
23 23
h ⎛ Pr ⎞ h ⎛ DAB ⎞
kc = ⎜ ⎟ = ⎜ ⎟ (6.6-12)
ˆ
ρ C p ,m ⎝ Sc ⎠ ˆ
ρ C p ,m ⎝ α ⎠
Solo nos resta calcular las propiedades que aparecen en la Ec.(6.6-12). Notemos que hemos
supuesto condiciones isotérmicas en el problema, por lo cual la temperatura promedio para
evaluar las propiedades es también 300 ºK. En vista de que la mezcla aire-amoníaco se
considera es diluida, podemos considerar que las propiedades de la mezcla se pueden
aproximar por las del aire:
PM aire (1 atm)(28.8 g/mol)

ρ mezcla ≅ ρaire = = 3
−3
= 1.17 × 10 g/cm
3
RT 82.05 atm ⋅ cm /mol ⋅ K)(300 K)
−5
kT (6.25 × 10 cal s ⋅ cm ⋅ K)
α mezcla ≅ α aire = = 2
= 0.22 cm s
ˆ
ρ C p ,m ( 1.17 × 10
−3
g cm
3
) (0.24 cal g ⋅ K)
Sustituyendo valores en la Ec.(6.6-12)
23
2
(0.108 cal cm ⋅ s ⋅ K) ⎛ 0.23 cm 2 s ⎞
kc = −3 3 ⎜ 2 ⎟ = 396.2 cm s
(1.17 × 10 g cm ) (0.24 cal g ⋅ K) ⎝ 0.22 cm s ⎠
Ejemplo 6.6-2. El hidrómetro de bulbo húmedo. Cálculo de la humedad

relativa
Para el control de diversos procesos y para conocer condiciones climáticas, es necesario

conocer la humedad relativa del aire. La humedad relativa se define como el cociente de la
presión parcial (o fracción mol) del vapor de agua en el aire a una determinada temperatura
sobre la presión de saturación ( o fracción mol de saturación) del agua a esa misma
temperatura:
pAb Tb yAb Tb
Hr = = (6.6-13)
pAsat Tb yAsat Tb
El higrómetro (o hidrómetro) de bulbo húmedo es un dispositivo sencillo y barato, cuyas

mediciones sirven de base para determinar la humedad relativa. Este higrómetro consiste
esencialmente de dos termómetros estándar (o termopares), uno de los cuales está
cubierto de una tela porosa que se mantiene mojada por agua contenida en un pequeño
recipiente. El termómetro que está descubierto mide la temperatura del aire que fluye y el
termómetro que está cubierto por la tela mojada mide una temperatura más baja debido al
R. LOBO/2007
enfriamiento que provoca la evaporación del agua causada por el flujo de aire. El
termómetro descubierto se denomina de bulbo seco y el que está cubierto por la tela
mojada se llama de bulbo húmedo. En la Figura 6.6-1 se muestra una instalación típica.
El problema consiste en desarrollar una expresión para el cálculo de la humedad relativa a

partir de las lecturas de los termómetros de bulbo húmedo y bulbo seco.
termómetro de termómetro de
bulbo seco bulbo húmedo
Tb
superficie cilíndrica de
Tí transferencia de calor y masa
aire
depósito de tela mojada

agua a Ti con agua
Figura 6.6-1. Representación de un hidrómetro de bulbo húmedo
Solución. Comencemos observando que el flujo de agua que se evapora alrededor del
termómetro de bulbo húmedo es el flux interfacial multiplicado por el área lateral de la
superficie cilíndrica expuesta al flujo de aire.
( ) (
WAi = π Dt Lt N A i = π Dt Lt kc CA i − CA b ≅ π Dt Lt kc Cm yA i − yAb ) (6.6-14)
donde Dt y Lt son el diámetro y la longitud del termómetro de bulbo húmedo expuesto a la

corriente de aire. Notemos con cuidado que CAi es la concentración de saturación de vapor
de agua en la superficie de la tela del termómetro de bulbo húmedo, cuya temperatura es Ti;
CAb es la concentración del vapor de agua que prevalece en el seno del aire y está a Tb.
Observemos también que hemos aproximado la concentración molar total C a cada
temperatura con una concentración molar total promedio Cm entre las temperaturas Tb y Ti.
Las fracciones mol yA tienen un significado análogo a las concentraciones.
El flujo convectivo de calor desde el aire hasta la superficie húmeda de la tela que cubre al
termómetro de bulbo húmedo es el mismo que se usa para vaporizar el agua de la tela.
Entonces,
R. LOBO/2007
π Dt Lt qi = π Dt Lt h (Tb − Ti ) = π Dt Lt N A i ∆ H A vap
(6.6-15)
( )
= π Dt Lt kc Cm yA i − yA b ∆ H A vap
donde hemos usado el balance de materia (6.6-14) y ∆ H Avap es el calor molar de

vaporización del agua. Rearreglando la Ec.(6.6-15),
h (Tb − Ti )
yA b = yA i − (6.6-16)
Ckc ∆ H A vap
La relación entre los coeficientes de transferencia de calor y de masa podemos obtenerla de

la analogía de Chilton-Colburn, Ec.(6.6-11), de donde podemos obtener:
23
h ⎛ α ⎞
= ρ m Cˆ p ,m ⎜ ⎟ = ρ mCˆ p ,m Le2 3 (6.6-17)
kc ⎝ DAB ⎠
donde Le es el número de Lewis y ρm es la densidad media entre Tb y Ti. Si hacemos uso

del hecho de que ρ mCˆ p ,m = Cm C p ,m , podemos sustituir la Ec.(6.6-17) en la (6.6-16) y
rearreglar para obtener:
⎛ C p ⎞ 2 3
yA b = yA i − ⎜ ⎟ Le (Tb − Ti ) (6.6-18)
⎜ ∆H A vap ⎟
⎝ ⎠
donde C p,m es ahora la capacidad calorífica molar media entre Tb y Ti. Si dividimos la
Ec.(6.6-18) por la fracción mol de saturación del vapor de agua a Tb, obtendremos la
definición de la humedad relativa, la cual se puede escribir para mayor claridad como
⎛ humedad ⎞ yAb yAsat ⎛ C p ⎞

⎜ relativa ⎟ = H r = sat =
Ti
−⎜ ⎟ Le 2 3 (Tb − Ti ) (6.6-19)
⎝ ⎠ yA Tb yAsat ⎜ ∆ H A vap yAsat ⎟
Tb ⎝ Tb ⎠
o bien, expresado en términos de presiones:
⎛ humedad ⎞ pAb pAsat ⎛ PC p ⎞

⎜ relativa ⎟ = H r = sat =
Ti
− ⎜ ⎟ Le 2 3 (Tb − Ti ) (6.6-20)
⎝ ⎠ pA Tb pAsat ⎜ ∆ H A vap pAsat ⎟
Tb ⎝ Tb ⎠
La Ec.(6.6-19) o la (6.6-20) cumple nuestro objetivo. Pueden utilizarse para un cálculo

confiable de la humedad relativa, a partir de las mediciones de los termómetros de bulbo
R. LOBO/2007
húmedo y bulbo seco. Incluso pueden usarse para sustancias que no sean agua, es decir,
para mezclas de un gas no condensable con otro que sí lo es. Con frecuencia se elimina de
estas ecuaciones en número de Lewis, ya que para el sistema aire-vapor de agua este
número es cercano a la unidad.
6.7 Transferencia de masa a través de interfases

En las secciones anteriores de este capítulo nos hemos referido a la transferencia de masa
que ocurre desde una interfase hacia una fase homogénea. Sin embargo, existen muchas
operaciones industriales en las cuales la transferencia de masa se presenta desde una fase
fluida bien mezclada hacia otra, a través de una interfase. Ejemplos de ello son la
absorción, la destilación, la extracción líquido-líquido y los reactores químicos. En pocas
palabras, para que ocurra la transferencia de un compuesto desde una fase hasta la otra,
deben ocurrir dos procesos de transferencia de masa, uno en cada fase involucrada y en
proceso ocurre a través de una interfase.
En la transferencia de masa entre dos fases, la resistencia al proceso se encuentra localizada

en la vecindad de la interfase. Sin embargo, ahora tendremos resistencia a la transferencia
en dos regiones, una para cada fase en consideración. Esta situación se ilustra en la figura
4.5-1 para una interfase gas-líquido. Lewis y Whitman1 desarrollaron estos conceptos con
el nombre de la teoría de las dos películas.
Para fines de ilustración consideraremos una interfase gas-líquido, pero el desarrollo es

también válido para una interfase líquido-líquido. En la Figura 6.7-1 se ilustra un punto
local P en una columna de absorción que contiene a una interfase. Un soluto A se transfiere
desde el seno del gas hacia la interfase; ahí se disolverá en el líquido y se transferirá hacia
el seno de éste. La presión parcial local del soluto en el seno del gas, donde se encuentra
bien mezclado, es pAb; la correspondiente en la interfase del lado del gas es pAi. La
concentración local del soluto A en la interfase del lado del líquido será CAi y en el seno de
éste será CAb. Debido a que la zona interfacial es delgada, podemos considerar que todo el
proceso de transferencia de masa se da en estado estacionario.
El flux local de soluto A que arriba a la interfase del lado del gas está dado por
(
N A i ,G = kG pA b − pA i ) (6.7-1)
donde kG es el coeficiente local de transferencia de masa en la fase gaseosa, con unidades

típicas de (mol/cm2·s·atm). Como veremos adelante, en problemas de transferencia de masa
1 Lewis, WK y Whitman, WG (1924), Ind. Eng. Chem., 16, 1215.
R. LOBO/2007
a través de interfases es conveniente denominar al coeficiente de transferencia de masa en

cada fase coeficiente individual, que en este caso es un coeficiente local individual.
gas libre de soluto

seno del líquido
CAb ( z )
líquido puro gas líquido
pA b
empaque pAi
inerte
CA i
gas con soluto seno del gas CAb

pAb ( z ) N Ai
líquido con soluto
Figura 6.7-1. Transferencia de masa a través de una interfase gas-líquido
El flux local de soluto A que se transfiere desde la interfase del lado del líquido hacia el
seno del líquido está dado por
(
N A i , L = k L CA i − CA b ) (6.7-2)
donde kL es el coeficiente local individual de transferencia de masa de la fase líquida, con

unidades típicas de (cm/s), y es idéntico a kc. La transferencia de masa entre las fases se da
en estado estacionario, por lo que
N A i ,G = N A i , L = N A i
(6.7-3)
( ) (
kG pA b − pA i = k L CA i − CA b )
De esta ecuación podemos obtener:
kL p − pAi
− = Ab (6.7-4)
kG CA b − CA i
En la Figura 6.7-2 aparece la curva de equilibrio para el soluto A en las fases gas y líquido
(recordar que la pendiente de esta curva representa al coeficiente de partición m). Las
R. LOBO/2007
condiciones en un punto del seno del gas y del líquido están representadas por el punto P1
(ver Figura 6.7-1). Como lo indica la Ec.(6.7-4), una línea recta que parta del punto P y
tenga pendiente −kL/kG tocará la línea de equilibrio en el punto I, que representa las
condiciones locales en la interfase. Las condiciones en la interfase caen sobre la curva de
equilibrio debido a la hipótesis de equilibrio en la interfase que usamos para relacionar las
concentraciones en las fronteras de dos fases. En este diagrama también quedan
establecidas las fuerzas motrices locales para cada una de las transferencias en las fases
individuales: (pAb−pAi )para la fase gas, y (CAi−CAb) para la fase líquida.
pA cuerva de
P equilibrio
pA b pA = mCA∗
fuerza motriz en el gas
kL
pendiente = −
fuerza motriz global
kG
pAi
I
pA∗
CA b CAi CA∗
CA
fza. motríz líquido
fuerza motriz global
Figura 6.7-2. Fuerzas motrices locales para la transferencia de masa local entre fases.
No obstante, las condiciones en la interfase, sea del lado del gas o del líquido, son difíciles
de evaluar o medir experimentalmente y, por lo general, son desconocidas. Sería muy
deseable poder cuantificar la transferencia del soluto del seno del gas al seno del líquido
una fuerza motriz que fuera observable experimentalmente (medible). Para alcanzar este
objetivo, es necesario eliminar pAi y CAi de las Ecs.(6.7-1) y (6.7-2). La forma de hacerlo
es usar la hipótesis de equilibrio en la interfase, suponer que en la interfase se alcanza
instantáneamente el equilibrio, de modo que las concentraciones interfaciales de ambas
fases están relacionadas por
pA i = mCA i (6.7-5)
1 El punto P está fuera de la línea de equilibrio, indicando que en esa localidad el sistema no está en equilibrio
y que, por lo tanto, habrá un proceso de transferencia de masa.
R. LOBO/2007
donde m es el coeficiente de partición de A entre el gas y el líquido (la constante de Henry

en gases diluidos en líquidos). Entonces, podemos sustituir la Ec.(6.7-5) en la (6.7-3) y de
la expresión resultante podemos despejar CAi:
kG pA b + k L CA b
CA i = (6.7-6)
mkG + k L
Podemos ahora sustituir esta ecuación en la (6.7-2) y después de simplificar obtendremos:
1
NAi =
1 m ( pAb − mCAb ) (6.7-7)
+
kG k L
Esta ecuación merece varios comentarios. En primer lugar, notamos que ahora el flux local
está expresado en términos de las concentraciones en los senos del gas y del líquido locales,
las cuales son observables experimentalmente. En segundo lugar, la Ec.(6.7-7) expresa el
flux local en términos de la suma de las resistencias a la transferencia de masa en cada una de
las fases; en otras palabras, hemos logrado expresar el flux en términos de la resistencia local
total a la transferencia de masa dada por (1/kG + m/kL).
Es importante percatarnos de que las resistencias 1/kG y 1/kL individuales en cada fase no
son directamente aditivas, sino que están relacionadas con el coeficiente de partición m, que
es una propiedad termodinámica del soluto en las fases en contacto.
Sería deseable ahora poder expresar el flux que cruza la interfase en términos de un solo
coeficiente de transferencia de masa que incluya las resistencias en el gas y en el líquido,
así como una sola diferencia de concentración en una de las fases. A este tipo de coeficiente
se le denomina coeficiente de transferencia de masa global, para distinguirlo de los
coeficientes individuales en cada una de las fases. Si escogemos expresar la Ec.(6.7-7) en
términos de una fuerza motriz de presiones parciales podemos escribir
(
N A i = KG pA b − pA∗ ) (6.7-8)
donde
1 1 1 m
KG = ó = + (6.7-9)
1 m KG kG k L
+
kG k L
y
pA∗ = mCA b (6.7-10)
R. LOBO/2007
En estas ecuaciones, KG es el coeficiente de transferencia de masa local global basado en

una fuerza motriz de presiones parciales; y p∗A es la presión parcial hipotética del
compuesto A que estaría en equilibrio con un líquido de composición CAb. Observemos
con cuidado que, fijadas la presión y la temperatura, el equilibrio de A entre las fases
gaseosa y líquida es único y, por lo tanto, pA∗ es tan buena medida de CAb como CAb
misma. Notemos también que ( pAb − pA∗ ) indica cuán alejado está el sistema localmente del
equilibrio. Esto se hace evidente en la Figura 6.7-2.
Al inicio de la derivación anterior podríamos haber decidido alternativamente expresar los

resultados en términos de una fuerza motriz local del lado del líquido y llegar a la expresión
(
N A i = K L CA∗ − CA b ) (6.7-11)
donde
1 1 1 1
KL = ó = + (6.7-12)
1 1 KL mkG k L
+
mkG k L
y
pA b
CA∗ = (6.7-13)
m
En estas ecuaciones, KL es el coeficiente de transferencia de masa local global basado en una

fuerza motriz en el líquido; y C ∗A es la concentración hipotética de soluto en el líquido que
estaría en equilibrio con el gas conteniendo una presión parcial del soluto igual a pAb.
Nuevamente, podemos decir que una vez fijadas la temperatura y la presión, C ∗A es tan buena
medida de pAb como pAb misma. Notemos además que ( CA∗ − CA b ) representa el grado de
alejamiento del sistema con respecto al equilibrio. Estos hechos pueden constatarse con un
cuidadoso examen de la Figura 6.7-2.
Las formas de expresar el flux local dadas por las Ecs.(6.7-8) y (6.7-11) son totalmente
equivalentes. Conviene destacar que KG y KL no son independientes, ya que
K L = mKG (6.7-14)
Los coeficientes locales globales KG y KL con frecuencia se denominan, respectivamente,

coeficiente de transferencia de masa global del lado del gas y coeficiente de transferencia de
R. LOBO/2007
masa global del lado del líquido, aunque es muy importante no olvidar que ambos
coeficientes incluyen las dos resistencias a la transferencia de masa.
En sistemas donde el soluto A presente en el gas es muy soluble en el líquido (NH3 en agua,
por ejemplo), m es pequeña y de la Ec.(6.7-9) se obtiene,
KG ≈ kG ; m pequeña; compuesto A muy soluble (6.7-15)
y entonces se dice que el proceso global de transferencia de masa se encuentra controlada

por la fase gaseosa. En contraste, en si el soluto A presente en el gas es poco soluble en el
líquido (C02 en agua, por ejemplo), m es grande y de la Ec.(6.7-12) podemos obtener,
K L ≈ k L ; m grande; compuesto A poco soluble (6.7-16)
y entonces el proceso total de transferencia de masa se encuentra controlado por la fase

líquida.
La razón de la resistencia a la transferencia de masa local en una de las fases respecto a la

resistencia total local puede calcularse como:
resistencia en la fase gaseosa 1 kG

= (6.7-17)
resistencia total 1 KG
resistencia en la fase líquida 1 kL

= (6.7-18)
resistencia total 1 KL
Es muy importante observar que los coeficientes individuales de transferencia de masa kG y

kL en solución diluida son prácticamente independientes de concentración. Sin embargo,
los coeficientes globales KG y KL sí dependen de ella. Esta dependencia proviene de que los
coeficientes globales son función del coeficiente de partición m y éste, a su vez, depende
de la concentración. Sólo cuando la línea de equilibrio es recta (m=constante) los
coeficientes globales KG y KL son independientes de la concentración.
Ejemplo 6.7-1. Estimación de coeficientes de transferencia de masa

individuales y globales a partir de otros conocidos o más fáciles de obtener
A menudo el ingeniero se ve delante del problema de estimar coeficientes de transferencia

de masa para una sustancia X en un equipo específico de contacto entre dos fases (digamos
gas-líquido) y no se conoce el comportamiento de este sistema, ni qué porcentaje de la
resistencia a la transferencia de masa se localiza en el gas ni en el líquido. El problema
consiste en estimar los coeficientes de transferencia de masa para la sustancia X a partir de
otras cuyo comportamiento sí se conoce.
R. LOBO/2007
Solución. Para estimar el coeficiente en el gas, podemos realizar la experimentación no con

la sustancia X, sino con una sustancia que sabemos que presenta la mayor resistencia en el
gas, digamos NH3, en el mismo aparato, con el mismo empaque y a las mismas condiciones
en que necesitamos los coeficientes para la sustancia X. Los datos de absorción del NH3 los
podemos correlacionar con una expresión del tipo
( kG ) NH3 Z
Sh G,NH
Z
3 = = AG Re aZ ,G Sc1NH
3
3G
(6.7-19)
DNH3G
En esta expresión, Z es la distancia característica adecuada, que puede ser el tamaño

nominal del empaque. La constante AG es una constante de ajuste y el número de Reynolds
del gas será únicamente dependiente de las condiciones hidrodinámicas seleccionadas. El
exponente (1/3) en el número de Schmidt es típico (teoría de la capa límite) en correlaciones
donde la resistencia a la transferencia de masa se localiza en el gas. Es de esperarse que una
correlación del coeficiente de transferencia de masa para la sustancia X en el gas siga el
mismo comportamiento, sólo que el número de Schmidt y la difusividad serán de la
sustancia X en el gas, esto es:
( kG ) X Z
Sh G,X
Z = = AG Re aZ ,G Sc1X3G (6.7-20)
DXG
Entonces, si hacemos el cociente entre las Ecs.(6.7-20) y (6.7-19) y rearreglamos,

obtendremos la solución a nuestro problema de estimar kG para la sustancia X:
13 23
( kG )X D ⎛ Sc XG ⎞ ⎛ D ⎞
= XG ⎜ ⎟ ≈ ⎜ XG ⎟ (6.7-21)
( kG ) NH3 DNH3G ⎜ Sc NH G
⎝ 3
⎟
⎠
⎜ DNH G
⎝ 3
⎟
⎠
Podemos proceder de manera semejante para el líquido, sólo que ahora seleccionaremos
para los experimentos un gas poco soluble, digamos CO2, para que la resistencia a la
transferencia de masa se localice en el líquido. La correlación para el CO2 es del tipo
( k L )CO2 Z
Sh Z ,CO2 = = AL RebZ ,L Sc1CO2 2L (6.7-22)
DCO2L
En esta expresión, Z sigue siendo la misma distancia característica usada para el gas. La
constante AL es una constante de ajuste y el número de Reynolds del líquido será
únicamente dependiente de las condiciones hidrodinámicas seleccionadas. El exponente
(1/2) en el número de Schmidt es típico en correlaciones (teoría de penetración) donde la
resistencia a la transferencia de masa se localiza en el líquido. Es de esperarse que una
correlación para la sustancia X en el líquido siga el mismo comportamiento, sólo que el
R. LOBO/2007
número de Schmidt y la difusividad serán de la sustancia X en el líquido, lo cual lleva a

proponer que:
( k L )X Z
Sh L,X
Z = = AL Re aZ .L Sc1X2L (6.7-23)
DXL
Entonces, si hacemos el cociente entre las Ecs.(6.7-23) y (6.7-22) y rearreglamos,

obtendremos la solución a nuestro problema de estimar kL para la sustancia X:
12 12
( k L )X D ⎛ Sc ⎞ ⎛ DXL ⎞
= XL ⎜ XL
⎟ ≈⎜ ⎟ (6.7-24)
( k L )CO2 DCO2L ⎜⎝ ScCO2LG ⎟
⎠
⎜ DCO L
⎝ 2
⎟
⎠
En las Ecs.(6.7-21) y (6.7-24) hemos supuesto que los cocientes de las viscosidades
cinemáticas presentes en los números de Schmidt son aproximadamente iguales a la unidad
ν XG ν XL
≈1 ; ≈1 (6.7-25)
ν NH3 G ν CO2 L
lo cual es una buena aproximación para soluciones diluidas.
Los resultados alcanzados deben usarse con cautela, pues sólo representan estimaciones.
Con frecuencia, en la literatura se encuentra, en ambas estimaciones, que las razones de
difusividades están elevadas al exponente (1/2). La razón para ello es que la estimación se
hace a partir de la teoría de penetración del coeficiente de transferencia de masa, lo cual
vimos en la Sec. 6.2 y volveremos a ver en el Capítulo 7. La explicación para el exponente
(2/3) en el resultado la Ec.(6.7-21) proviene de la teoría de la capa límite, que también
veremos en el Capítulo 7. En términos prácticos, la diferencia es menor y bien podemos
usar el exponente (1/2) en los resultados expresados por las Ecs.(6.7-21) y (6.7-24).
Por último, para poder calcular el coeficiente global para la sustancia X debemos contar
con el coeficiente de partición mX, el cual debe obtenerse de experimentos al equilibrio.
Sólo entonces, podemos obtener el coeficiente global de las definiciones dadas por la Ec.
(6.7-9) o la (6.7-12). Es importante señalar, además, que los coeficientes empleados en este
ejemplo son locales, si bien el procedimiento y los resultados son los mismos para obtener
los coeficientes promedio.
Ejemplo 6.7-2. Absorción de H2S en agua
En una torre empacada que opera a 30 ºC y 2 atm, se lleva a cabo un estudio de absorción
de H2S de una corriente de aire en agua. Un valor estimado del coeficiente de transferencia
R. LOBO/2007
de masa individual del lado del gas es de 9.31×10−2 cm/s. En un punto determinado de la
torre, la fracción mol del H2S en la interfase del lado del líquido es 2×10−5, y la presión
parcial del H2S en el seno del gas es de 0.05 atm. El coeficiente de partición (la constante
de Henry) del sistema es 600 atm/(fracción mol).
Calcula:
(a) el flux de H2S que se absorbe en agua;
(b) el coeficiente de transferencia de masa individual kG;
(c) el coeficiente individual de transferencia de masa del lado del líquido, kL;
(d) el coeficiente de transferencia de masa global del lado del líquido KL; y
(e) el porcentaje de la resistencia a la transferencia de masa que se localiza en la fase
líquida.
Solución. (a) Designemos al H2S como la especie A. Su tasa de absorción está dada por
(
N AGi = kcG CAGb − CAGi ) (6.7-26)
donde hemos usado el superíndice G para indicar la fase gaseosa. Notemos que por
definición de la fuerza motriz kcG ≠ kG ; mientras kcG usa una diferencia de concentración
molar como fuerza motriz, kG usa una diferencia de presión parcial.
La concentración de A en el seno de la fase gaseosa es
pA b (0.05 atm)
CAGb = = 3
= 2.01 × 10
−6
mol cm
3
RT (82.05 atm ⋅ cm mol ⋅ K)(303K)
Necesitamos ahora conocer la concentración de A en la interfase del lado del gas.

Suponemos que en la interfase prevalece el equilibrio y que la relación entre las
concentraciones del gas y del líquido está dada por el coeficiente de partición m. Por las
unidades de éste, podemos ver que la relación es
pA i = m′ xA i = 600 xA i = (600 atm frac.mol)( 2 × 10 −5 frac.mol)

= 1.2 × 10 −2 atm
Entonces, la concentración en la interfase del lado del gas es
pA i (0.012 atm) 3
CAGi = = 3
= 4.83 × 10
−7
mol/cm
RT (82.05 atm ⋅ cm mol ⋅ K) (303K)
Ahora simplemente sustituimos valores en la Ec.(6.7-26) para obtener el flux de H2S que se
absorbe,
R. LOBO/2007
N AGi = (9.31 × 10−2 cm s)[(2.010 − 0.483) × 10−6 mol/cm]3

= 1.42 × 10−7 mol/cm 2s
(b) De la Ec.(1) podemos ver que
( )
N AGi = kcG CAGb − CAGi = kG pA b − pA i ( ) (6.7-27)
de donde es evidente que
∆CAG ∆ pA RT kcG
kG = kcG = kcG = (6.7-28)
∆ pA ∆ pA RT
Sustituyendo valores en la Ec.(6.7-28) podemos obtener el valor de kG buscado
(9.31 × 10 −2 cm s) −6
kG = 3
= 3.74 × 10
2
mol/cm ⋅ atm ⋅ s
(82.05atm ⋅ cm mol ⋅ K) (303K)
(c) Por definición, sabemos que el flux interfacial del lado del líquido es:
(
N AL i = k L CA i − CA b = N AGi ) (6.7-29)
En este caso, para calcular kL debemos hacer dos suposiciones. Primero, que CAb≈0, en
vista de que la solución es diluida y no conocemos ningún valor de la concentración de A
en el seno del líquido. La segunda suposición es que la concentración molar total de la
solución, CL, es aproximadamente igual a la del agua, lo cual es razonable porque la
solución es diluida. Entonces,
ρ H 2O (1 g/cm3 ) −5
CA i = CL xA i ≅ CH 2O xA i = xA i = (2 × 10 )
M H 2O (18 g/mol)
= 1.11 × 10 −6 mol/cm3
Una vez calculadas las concentraciones en el líquido, de la Ec.(6.7-29) obtenemos
N AL i (1.42 × 10
−7
mol/cm s)
2
kL = = = 0.13 cm/s
( CA i − 0 ) (1.11 × 10
−6
mol/cm )
3
(d) El coeficiente de transferencia de masa global del lado del líquido está dado por la
Ec.(6.7-12)
R. LOBO/2007
−1
⎛ 1 1 ⎞
KL = ⎜ + ⎟ (6.7-30)
⎝ mkG k L ⎠
Observemos con cuidado que el coeficiente de partición m en la Ec.(6.7-30) no tiene las

mismas unidades que el que hemos utilizado en este problema. En otras palabras, el dato
que nosotros tenemos en este ejemplo corresponde a una definición dada por pAi=m'xAi, en
tanto que el coeficiente de partición en la Ec.(6.7-30) proviene de una definición pAi=mCAi.
La relación entre ambas formas de expresar el coeficiente de partición está dada por
m′ m′ (600 atm)
m= = = 4 3
= 1.08 × 10 atm ⋅ cm mol
CL ρ L M L (1 g/cm 18 g/mol)
3
Ahora simplemente sustituimos los valores correspondientes en la Ec.(6.7-30) para obtener
−1
⎡ 1 1 ⎤
KL = ⎢ 4 3 −6 2
+ ⎥
⎣ (1.08 × 10 atm ⋅ cm mol)(3.74 × 10 mol cm ⋅ atm ⋅ s) (0.13 cm s) ⎦
= 0.031 cm s
(e) El porcentaje de la resistencia total a la transferencia de masa que se localiza en la fase

líquida es:
resistencia en la fase líquida 1 kL 1 (0.13 cm/s)

× 100 = × 100 = × 100 = 23.9
resistencia total 1 KL 1 (0.031 cm/s)
Recapitulación
La caracterización de la transferencia de masa en las interfases es un problema complejo
debido al desconocimiento de la situación del flujo alrededor de éstas. Ello hace necesario
introducir el modelo del coeficiente de transferencia de masa, que si bien es aproximado,
resulta muy útil. Este modelo es una alternativa al de la ley de Fick, ya que agrega en el
coeficiente de transferencia de masa diversos efectos que ésta no puede tomar en cuenta.
En la Sec. 6.1, planteamos una definición empírica del coeficiente de transferencia de masa.
Señalamos que existe una diversidad de definiciones de él y apuntamos que es necesario
definir con claridad la diferencia de concentración y el área de transferencia.
Posteriormente, en la Sec. 6.2 vimos cómo una serie de problemas que analizamos en los
capítulos 2 y 4 vía la ley de Fick tienen un equivalente si son analizados mediante el
coeficiente de transferencia de masa. En la Sec. 6.3 hicimos una aproximación teórica
R. LOBO/2007
cualitativa a la definición empírica describiendo cómo se agrupan en el coeficiente la

difusión molecular y la turbulenta, así como la convección causada por alta concentración
del soluto y/o altas tasas de transferencia. También demostramos con algunos ejemplos
cómo calcular este coeficiente directamente a partir de datos experimentales.
En las Sec. 6.4 vimos cómo valernos del análisis dimensional para correlacionar datos
experimentales en términos de números adimensionales y en la Sec 6.5 presentamos una
variedad de correlaciones comúnmente usadas para estimar el coeficiente de transferencia
de masa. En la Sec. 6.6 abordamos las analogías entre las transferencias de momentum,
calor y masa para calcular el coeficiente de transferencia de masa. Por último, en la Sec. 6.7
extendimos el uso del modelo del coeficiente de transferencia de masa a situaciones en las
cuales la transferencia ocurre entre dos fases; desarrollamos la llamada teoría de las dos
películas e ilustramos el uso de coeficientes de transferencia de masa globales para
cuantificar el flux de masa que se transfiere entre fases.
Globalmente, el contenido del capítulo es práctico. Nos da bases para analizar y realizar
cálculos en aplicaciones como disolución, evaporación, secado, destilación, absorción y
extracción, entre otras. Es útil también para describir la transferencia de masa en reactores
químicos, aunque no para entender las reacciones químicas que suceden dentro de estos
equipos.
Revisión de conceptos
R6.1. ¿Por qué es importante el estudio de la transferencia de masa en las interfases? Cita
tres ejemplos.
R6.2. Explica las razones por las que es necesario emplear modelos empíricos aproximados
para cuantificar la transferencia de masa en las interfases.
R6.3. Explica el modelo del coeficiente de transferencia de masa.
R6.6. Analiza, con base en ejemplos, los problemas que existen para lograr una definición
precisa del coeficiente de transferencia de masa.
R6.5. Explica cuáles diferencias de concentración son preferibles para definir el coeficiente
de transferencia de masa. ¿Por qué?
R6.6. Describe en qué tipo de sistemas no es posible obtener el coeficiente de transferencia

de masa, y por lo tanto se hace necesario obtener el producto ka. ¿Qué nombre recibe este
tipo de coeficiente?
R6.7. Explica el proceso de transferencia de masa ilustrado en la figura 6.2-1.
R. LOBO/2007
R6.8. ¿Qué es la difusión turbulenta? ¿Cuál es su origen? ¿Cómo se cuantifica? ¿De qué
depende la difusividad turbulenta?
R6.9. ¿Cuál es la expresión del flux total (difusión + convección) en la interfase en

términos del coeficiente de transferencia de masa? ¿A qué se reduce dicha expresión
cuando las soluciones son diluidas y la tasas de transferencia de masa bajas?
R6.10. Describe en qué se basan los métodos de obtención del coeficiente de transferencia
de masa.
R6.11. Define los principales números adimensionales de interés para transferencia de

masa. Explica su significado físico.
R6.12. ¿Qué es una correlación? ¿Cómo se obtiene?
R6.13. ¿Qué significa la afirmación de que una correlación implica una situación física
determinada?
R6.16. ¿Cuáles son las ideas básicas que sustentan el uso de las denominadas analogías
para obtener coeficientes de transferencia de masa?
R6.15. Explica las hipótesis en que se basa la analogía de Reynolds. ¿A qué tipo de
sistemas es aplicable y en cuáles falla totalmente?
R6.16. ¿Cuál es la analogía de Chilton-Colburn? ¿Cuáles son sus alcances y limitaciones?
R6.17. Distingue entre lo que es un coeficiente de transferencia de masa local, uno

promedio, uno individual y uno global.
R6.18. ¿Es siempre cierta la hipótesis de que en la interfase entre dos fases existe
equilibrio?
R6.19. En un sistema de dos fases, ¿son las resistencias a la transferencia de masa en cada
fase directamente aditivas? Explica.
R6.20. ¿Cómo depende el coeficiente de transferencia de masa global de la concentración?

¿Cuándo es independiente de ella?
R6.21 ¿Por qué en operaciones de contacto gas-líquido y líquido-líquido en lechos

empacados el área interfacial no es igual al área física del empaque?
R.22. Estudia de nuevo el Ejemplo 6.2-1 y explica el resultado expresado por la Ec.(6.2-42)
R. LOBO/2007
Problemas
6.1. Determinación del coeficiente de transferencia de masa. En un recipiente de 20 L,
inicialmente lleno de aire seco, se ponen 1.6 L de agua y ésta comienza a evaporarse. El
recipiente se mantiene a 25 ºC y el área superficial del agua dentro de él es de 150 cm2.
Después de 3 minutos el aire está 0.05% saturado de vapor de agua. ¿Cuál es el coeficiente
de transferencia de masa? ¿En cuánto tiempo se alcanzará una concentración del vapor de
agua del 90 % de su valor de saturación?
6.2. Determinación del coeficiente de transferencia de masa. Se tiene un tubo empacado

de esferas de ácido benzoico de 0.2 cm de diámetro; las esferas tienen 23 cm2 de superficie
por cada cm3 de lecho. Una corriente de agua fluye dentro del tubo a una velocidad
superficial de 5 cm/seg. Al final de 100 cm de lecho el agua está 62% saturada con ácido
benzóico. Calcula el coeficiente de transferencia de masa.
6.3. Determinación del coeficiente de transferencia de masa. En un recipiente se agitan

agua y bromo hasta formar una emulsión. La concentración del bromo en el agua alcanza
en 3 minutos un valor de alrededor del 50% del valor de saturación. Calcula el coeficiente
de transferencia de masa del bromo al agua.
6.6. Uso de correlaciones. Considera un lecho empacado con esferas de naftaleno de 1.0
pulg. de diámetro. Una corriente de aire seco pasa a través del lecho a una velocidad de 4
ft/s. Encuentra el coeficiente de transferencia de masa para la sublimación del naftaleno, si
la temperatura es 60 ºF y su coeficiente de difusión en aire es 0.21 ft2/hr.
6.5. Uso de correlaciones. Se derrama un combustible sobre una superficie plana y forma
una capa de 0.25 pulg. de espesor sobre una longitud de 7 ft. El combustible está a 60 ºF; a
esta temperatura su presión de vapor es 2 psia; su viscosidad cinemática es 1.05 ft2/hr, su
densidad es 54 lb/ft3 y su coeficiente de difusión en aire es 0.55 ft2/hr. Si una brisa a 75 ºF
sopla paralela a la superficie del combustible derramado con una velocidad de 5 millas por
hora, determina la tasa de evaporación del combustible. Repite el cálculo para velocidades
del viento de 10, 15 y 20 millas/hr.
6.6 Uso de correlaciones. Una corriente de aire pasa a una velocidad de 50 ft/s a través de
un tubo cilíndrico de 1.0 pulg. de diámetro interno hecho de naftaleno. La temperatura es
50 ºF y la presión 1 atm. A esta temperatura, la presión de vapor del naftaleno es de 0.0209
mm Hg y su difusividad en aire es de 0.2 ft2/hr. Determina a qué longitud del tubo la
concentración del naftaleno será de 3.7 × 10−8 lbmol/ft3.
6.7. Uso de analogías. Repite el problema 6.6 utilizando:
(a) la analogía de Reynolds.

(b) la analogía de Chilton-Colburn.
R. LOBO/2007
Analiza tus resultados a la luz de los postulados de cada analogía.
6.8. Determinación de la humedad relativa de una corriente de aire. Una corriente de

aire pasa por un ducto que contiene un termómetro de bulbo húmedo y uno de bulbo seco.
Las lecturas en estos termómetros son de 17 y 37 ºC, respectivamente. Determina la
humedad relativa de la corriente de aire.
6.9. Determinación de la tasa de evaporación de agua. Una corriente de aire fluye a una
velocidad de 1.20 m/s, tiene una temperatura de 25 ºC y una humedad relativa de 30 %.
Esta corriente fluye paralelamente a una superficie de agua que está a 15 ºC por un trayecto
de 20 cm. Calcula la tasa de evaporación del agua.
6.10. Absorción de amoníaco en agua. Se tiene un tubo vertical de 5.0 cm de diámetro

interno y 1.2 m de largo, dentro del cual desciende una película de agua de 0.2 mm de
espesor a una velocidad de 0.17 m/s. Dentro del espacio libre en el interior del tubo
asciende una corriente de aire con amoníaco. La fracción mol del amoníaco en la interfase
del lado del líquido es 0.06. Calcula la tasa de absorción del amoníaco en el agua.
6.11. Evaporación de una gota. Una gota esférica de agua de 0.04 cm de diámetro cae a
una velocidad de 2.0 m/s en aire seco a 1 atm. La temperatura del aire es de 60 ºC y la de la
superficie de la gota es 30 ºC. En estas condiciones no es posible ignorar la parte
convectiva del flux total (la mezcla aire-vapor de agua en la vecindad de la interfase no es
diluida). Suponiendo pseudo estado estacionario, calcula la tasa de evaporación de la gota y
el tiempo en que ésta se evaporará por completo.
6.12. Coeficientes de transferencia de masa globales. Deriva expresiones semejantes a

las desarrolladas en la Sección 6.5 para los coeficientes de transferencia de masa globales
utilizando como fuerzas motrices para el gas y el líquido las respectivas fracciones mol.
6.13. Coeficientes de transferencia de masa globales y fuerzas motrices. En una torre de

absorción se remueve SO2 de una mezcla aire-SO2 mediante su absorción en agua. La torre
opera a 1 atm de presión y 20 ºC. En un punto dentro de la torre se analiza el gas y el
líquido y se encuentra que la presión parcial del SO2 es de 30 mm Hg y la concentración en
la solución es de 0.0344 moles SO2/ft3. Los coeficientes individuales de transferencia de
masa son kG=0.295 lbmol/hr·ft2·atm y kL=1.3 ft/hr. En experimentos por separado se
determinaron datos de equilibrio gas-líquido de SO2 sobre soluciones acuosas de SO2,
obteniéndose los siguientes datos
pSO2 mm Hg 0.5 3.2 8.5 26 59
CSO2 lbmol/ft3 0.0191 0.0911 0.174 0.388 0.681
R. LOBO/2007
(a) Grafica los datos de equilibrio en papel milimétrico. Localiza en este diagrama el
punto de operación de la columna correspondiente a los análisis del gas y del líquido.
(b) Determina las concentraciones de la interfase gas-líquido correspondientes al punto de

operación. Calcula el valor de la constante de Henry en este punto.
(c) Localiza en el diagrama (pAb−pAi), (CAi−CAb), (pAb−pA*), (CA*−CAb).
(d) Calcula los coeficientes de transferencia de masa globales KG y KL.
(e) Calcula el porcentaje de la resistencia a la transferencia de masa del lado del líquido.
6.13 Coeficientes de transferencia de masa globales Deriva las Ecs. (6.7-11), (6.7-12) y
(6.7-13)
6.14 Coeficientes de transferencia de masa globales Demuestra la Ec. (6.7-14)
6.15 Coeficientes de transferencia de masa globales. Deriva las expresiones para los
coeficiente de transferencia de masa globales Ky y Kx, a partir de los coeficientes
individuales ky y kx y las respectivas fuezas motrices en fracción mol.
6.16 Coeficientes locales y coeficientes promedio. A partir de las ecuaciones que definen
los coeficientes de transferencia de masa locales dados por las Ecs.(6.2-18) y (6.2-28),
obtener el coeficiente promedio aplicando la Ec.(6.2-32). Compara tus resultados con las
Ecs.(6.2-21) y (6.2-30), respectivamente.
6.17 Coeficientes locales y coeficientes promedio. Se tiene un pequeño lecho cilíndrico de

60 cm de longitud, empacado de esferas de naftaleno (A) de 0.2 cm de diámetro, como el
ilustrado en la Figura 6.3-4. El área específica del lecho de las esferas empacadas es
aproximadamente de 25 cm2 por cada cm3 de volumen del lecho y se representa por a. A este
lecho se alimenta aire (B) puro y seco a 1 atm y 75 ºC, a una velocidad superficial de 2 m/s.
La presión de vapor del naftaleno a la temperatura del aire es 5 mm de Hg. El coeficiente de
transferencia de masa promedio en el lecho es 6.81×10−2 cm/s. Puede suponerse que las
esferas de naftaleno no cambian su tamaño.
Calcula el flujo molar de naftaleno a la salida del lecho. Sugerencia: Primero, define la
diferencia de concentración a emplear con el coeficiente de transferencia de masa promedio
dado como dato. Después, calcula la concentración promedio a la salida del lecho.
Finalmente, calcula el flujo molar de naftaleno a la salida del lecho.
R. LOBO/2007
R. LOBO/2007

Unidad Ii - Lobo

Cargado por

Copyright:

Formatos disponibles

Unidad Ii - Lobo

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Unidad Ii - Lobo

Cargado por

Copyright:

Formatos disponibles

PRINCIPIOS DE TRANSFERENCIA DE MASA/ 2a.

Así, por ejemplo, la operación de secado (interfase sólido-gas) transfiere un compuesto

Aunque apenas nos iniciamos en el estudio de la transferencia de masa convectiva, es

En este capítulo y el siguiente estaremos analizando problemas en los que la transferencia

6.1 Definición del coeficiente de transferencia de masa en una fase

WAi = k A∆CA (6.1-1)

Pasemos ahora a precisar la definición de k que adoptaremos en este texto. Refirámonos a

"congelado" un proceso típico de transferencia de masa interfacial. En este sistema existe un

Al dividir esta ecuación por el área interfacial Ai, obtenemos:

solución bien CAb

Observemos que ahora el coeficiente de transferencia de masa es kc, donde el subíndice c es

Definición Unidades del coeficiente de

La definición dada por la Ec.(6.1-4) es común, ya sea en un líquido o de un gas; es

La diferencia de concentración que impulsa la transferencia de masa (∆CA u otra) debe

tienen dificultades para definir correctamente esta diferencia de concentración,

La diferencia de concentración a utilizar es aquella que prevalece entre la superficie

Figura 6.1-2. Disolución de un sólido en un líquido. La diferencia de concentración [CAi−CAb(t)] varía en el

Consideremos primero un sistema donde la diferencia de concentración varía en el tiempo,

gas libre de soluto

gas con soluto seno del gas

líquido con soluto

Figura 6.1-3. Diferencias de concentración en la absorción de un soluto. Las diferencias de concentración en

Examinemos ahora un caso en el que la diferencia de concentración varía en el espacio, pero

Un parámetro esencial en los problemas de transferencia de masa convectiva es el área de

dAi área interfacial

donde el producto (kca) recibe el nombre de coeficiente de transferencia de masa

Podemos resumir nuestra discusión sobre la definición del coeficiente de transferencia de

Ejemplo 6.1-1. Cálculo de diferentes coeficientes de transferencia de masa

Un líquido A se evapora a 25 ºC y 1 atm de presión dentro de un tubo de 0.3 cm de

Solución. El sistema se encuentra representado en la Figura 4.4-1. Las definiciones de los

N A = kc ( CA1 − CA2 ) ; N A = k y ( yA1 − yA2 ) ; N A = kG ( pA1 − pA2 ) (6.1-12)

El flux total NA se puede obtener de dividir la rapidez de evaporación por el área de

pA1 = pAsat = 0.15 atm ; pA2 = 0.01 atm

RT 82.05 cm ⋅ atm mol ⋅ K × 298 K

yA1 = pA1 P = 0.15 atm 1 atm = 0.15

Ejemplo 6.1-2. Conversión entre diferentes coeficientes de transferencia de

De esta relación podemos obtener:

donde hemos usado el hecho de que para un gas CA = pA RT . Sustituyendo valores:

En cuanto a los coeficientes de transferencia de masa del líquido, podemos tratarlos de

Aquí hemos hecho uso de CA = CL xA , donde CL es la concentración molar total de la

Entonces, sustituyendo valores:

k x = ( 0.13 cm s ) ( 0.056 mol cm 3 ) = 7.28 × 10 −3 mol cm 2 ⋅ s

6.2 Algunas expresiones teóricas del coeficiente de transferencia de masa

Difusión en solución diluida y contradifusión equimolar en soluciones concentradas.

N A,z = J A,z (6.2-1)

Si igualamos la Ec.(6.2-2) con la Ec.(6.1-1), podemos hacer las siguientes equivalencias

( CAi − CAb ) = ( CA0 − CA L ) (6.2-3)

Como podemos apreciar en la Ec.(6.2-4), el coeficiente de transferencia de masa es

Difusión y convección en una película estancada en soluciones concentradas. La

C DAB ⎛ 1 − yAL ⎞ C DAB ⎛ yBL ⎞

La Ec. (6.1-1) no es directamente comparable con la Ec.(6.2-5) ya que en ésta la diferencia

C DAB (1 − yAL ) − (1 − yA0 ) ⎛ 1 − yAL ⎞

En la Ec.(6.2-6) podemos reconocer una diferencia de concentración media logarítmica,

(1 − yAL ) − (1 − yA0 ) = yBL − yB0

( CAi − CAb ) = ( CA0 − CA L ) (6.2-9)

La Ec.(6.2-10) es la relación buscada. Es importante notar que en este caso el coeficiente de

Difusión en un medio seminfinito en solución diluida (teoría de penetración). En este

Igualando esta ecuación con la Ec.(6.1-1), es evidente que

( CAi − CAb ) = ( CA0 − CA ∞ ) (6.2-12)