Lectura Complementaria 2 PDF

FUNDAMENTOS DE LA
TRANSFERENCIA DE MASA
LOS capítulos anteriores relacionados con los fenómenos de transferencia

se han referido a fases de una componente con tendencia natural a alcanzar
las condiciones de equilibrio.
Cuando un sistema contiene dos o más componentes cuyas concentra-
ciones varían de un punto a otro, presenta una tendencia natural a transferir
la masa, haciendo mínimas las diferencias de concentración dentro del siste-
ma.
La transferencia de un constituyente de una región de alta concentraci6n
a una de baja concentración se llama transferencia d e masa.
Muchas de nuestrasexperienciascotidianasestánrelacionadascon la
transferencia de masa. Un terrón de azúcar en una taza de café negro se di-
suelve finalmente y después se difunde de manera uniforme en el café. El
a p a se evapora de l o s estanques, incrementando la humedad de la corriente
de aire que pasa por el lugar. Del perfume emana una agradable fragancia, que
se esparce por la atmósfera circundante.
L a transferencia de masa juega un papel muy importante en muchos pro-
cesos industriales: la r e m o s i h d e materiales contaminantes de las corrientes
de descarga de los gases del agua contaminada, l a difusión de neutrones den-
tro delos reactores nucleares, la difusión de sustancias quelos poros del carbbn
activado absorben, la rapidez de las reacciones químicas catalizadas y biológi-
cas así como el acondicionamiento del aire, son ejemplos típicos.
Si se piensa en el terrón de azúcar puesto dcntro de l a taza de café ne-
gro, la experiencia nos enseña queel intervalo de tiempo que se requiere para
destribuir el azúcar depende desi el líquido está en reposo0 se le agita m&-
nicamente por medio de una cucharita.
El mecanismo de transferencia de masa, tal como se ha observado en el
d e transferencia de calor, depende de l a dinámica del sistema en el que se 1 1 ~ -
va a cabo.
533
534 Fundamentos dela transferencia demasa
La masa puede transferirse por medio del movimiento molecular fortuito

en los f-luidos en reposo o puede transferirse de una superficie a un fluido en
movimiento, ayudado por las características dinámicas del flujo.
Estos dos modos diferentes de transferencia de masa: molecular y con-
vectiva, son análogos a la conduccibn calorífica y a la transferencia convectiva
de calor.
Cada uno de estos modos de transferencia de masa se va a describir y
analizar.Igual que en el caso de la transferencia de calor, inmediatamente
debe uno percatarse de que ambos mecanismos actúan, a menudo simultánea-
mente.Sinembargo, al ocurrirambos,unopuede ser cuantitativamente
dominante, de manera que es necesario usar soluciones aproximadas que in-
cluyan solamente el modo dominante.
24.1 T R A N S F E R E N C I A D E M A S A M O L E C U L A R
Ya en el aiio de 18 15 Parrot observi, cuantitativamentequecuando

una mezcla de gases contiene dos o más especies moleculares cuyas concentra-
ciones relativas varían de un punto a otro, resulta un proceso, aparentemente
natural, que tiende a disminuir cualesquiera desigualdades de composicibn.
Estatransferenciamacroscópicademasa,independientedecualquier
convección que se lleve a cabo dentro del sistema, se define con el nombre de
difusión melocular.
En el caso específico delas mezclas gaseosas se puede deducir una expli-
caci6n l6gica de este fenómeno de transferencia a partir de la teoría cinética
de los gases.
A temperaturas superiores al cero absoluto, las moléculas individuales se
encuentran en un estado de movimiento continuo, aunque fortuito. Dentro
de las mezclas de gases diluidos, cada una de las moléculas de soluto se com-
portaenformaindependientede las otrasmoléculasdesoluto,yaque
rara vez se topa con ellas. Estin ocurriendo contlnuamente colisiones entre
el solvente y el soluto. Como resultado de estas colisiones, las moléculas del
soluto describen trayectorias en zigzag, a veces hacia una regibn de mayor
concentración, a veces hacia una de concentración m i s baja.
Examinemos una seccibn hipotktica que pase en forma normal al gra-
diente de concentración dentro de una mezcla gaseosa isobárica e istérmica
que contenga moléculas de solvente y soluto.
I,os dos elementos delgados e iguales de volumen que se encuentran so-
bre y por debajo de la sección, contienen el mismo número de moléculas, tal
Como lo estipula la ley de Avogadro.
Aunque n o es posible establecer la dirección específica en la queviajará
una mol6cula particular en un intervalo dado de tiempo, puede decirse que
un número definido de moléculas que se encuentren enel elemento inferior de
volumen cruzará la sección hipotética desde abajo y el mismo nimero de mo-
léculas abandonar6 el elemento superior y atravesará la sección desde arriba.
Transferencia demasa molecular 535
Con la existencia del gradiente de concentración, habrá más moléculas

de soluto en uno delos elementos de volumen que en el otro; así pues, resul-
tará una transferencia total neta de unaregión de concentracihn nlayor a una
de cocentración menor.
El flujo de cadaunade las especiesmolescularesocurreen la direc-
ción del gradiente negativo de la concentracibn.
Como se señaló en los Capítulos 7 y 15, la transferencia molecular de
momento y la transferencia de energía por conducción se deben también al
movimiento al azar de las moléculas. De manera que es de esperarse que los
tres fenbmenos de transferencia dependan de muchas propiedades y caracte-
rísticas iguales, tales como la de trayectoria media libre y el hecho de que los
análisis tehricos de los tres fenómenos tienen muchoen común.
LA ECUACION DE RAPIDEZ DE FICK
Las leyes de transferencia demasa ponen de manifiesto la relacibn entre

el flujo de la sustancia que se esti difundiendoy el gradiente de concentración
responsable de esta transferencia de masa. Desafortunadamente, la descrip-
ción cuantitativa de la difusión molecular es considerablemerlte mis comple-
ja que las dcscripciones análogas correspondientes a la transferencia molecular
de momento y energía que tienen lugar en una fase de una componente.
Como la transferencia de masa o difusión, como se le llama también,
ocurre solamente en mezclas, su evaluación debe inclulr un examen del efec-
to de todaslas componentes. Por ejemplo, a menudo se desea conocer la rapi-
dez de difusión de una componente específica relacionada conla velocidad de
la mezcla en la cual se esta moviendo. Conlo cada una de las componentes
puede poseer unamovilidad diferente, la velocidad de la mezcla debe evaluar-
se promediando las velocidades de todas las componentes que se encuentren
presentes.
Parapoder establecer una base común para estudios futuros, estudiemos,
primero, las definiciones y relaciones que se utilizan a menudo para explicar
el papel de los componentes de la mezcla.
Concentraciones. En una mezcla múltiple de componentes la concentra-
ción de una especie molecular se puede expresar de muchas maneras. En la fi-
gura 24.1 aparece un volumenelemental dL', quecontieneuna mezcla de
componentes, incluyendo a la especie A . Comocadauna de las moléculas
de cada especie tiene una masa, puede definirse una concentracibn de masa
para cada especie, así como para la mezcla.
La concentraci6n de masa p . 4 , correspondiente a la especie A , se define
como la masa de A por unidad de volumen de la mezcla. La conccntraci6n
total de masa O densidad, p , es la masa total de la mezcla contenida en la uni-
dad d e volumen, esto es:
c Pt
n
P= (24-1)
I = I
536 Fundamentos de la transferencia demasa
Figura 24.1 Volumenelemental que contieneunamezcladecomponentesmúltiples.
donde n es el número de especies presentes en la mezcla.La fracción de masa,

W A , es la concentración de la especie A , dividida entre la densidad total de
masa:
(24-2)
Por definción, la suma de las fracciones de masa, debe ser uno.

n
w,=1 (24-3)
t=l
La concentración molar de la especie A , C A , se define como el número

de moles de A,presentes por unidad de volumen dela mezcla. Por definición,
un mol de cualquier especie contiene una masa equivalente a su peso mo-
lecular.. Los términos de l a concentración de masa y de l a concentracih molar
están relacionados por medio de la siguiente expresión:
(24-4)
dende MA es el peso molecular de la especie A.

Cuando se está tratando con una fase gaseosa, a menudo,las concentra-
ciones se expresan en términos de las presiones parciales. Bajo las condiciones
en las cuales es válida la ley de los gases ideales, PA V=nA RT, la concentra-
ción molar es:
nA PA
C A =-=- (24-5)
V RT
donde PA es la presión parcial de la especie A en la mezcla, n A es el número
de moles de la especie A, I/ es el volumen de gas , T es la temperatura absolu-
ta y R es la constante delgas. La concentración molar total,c , es el número to-
tal de moles de la mezcla, contenidos en la unidad de volumen, esto es:
n
c= c, (24-6)
i=l
Transferencia de masa molecular 537
o, en una mezclagaseosa que obedezca la ley de los gases ideales:

cat, tal/VzP/RT, donde P es la presión total. La fracción molar correspon-
diente a las mezclas de líquidos o sólidos, x A , y la correspondiente a mezclas
gaseosas, Y A , son las concentraciones molares de la especie A divididas entre
la densidad molar total.
CA
XA = - (líquidos y sólidos)
C
(24-7)
CA
yA = - (gases)
C
Cuando una mezcla obedece la ley de

los gases ideales, la fracción molar,
YA, se puede escribir en f u n c i h d elas presiones:
(24-8)
La ecuación (24-8) es una representación algebraica de la ley de Dalton

que corresponde a mezclas de gases. La suma de las fracciones molares, debe
ser igual a 1, por definción:
n
x,=1
¿= 1
(24-9)
n
yi=1
1=1
En la Tabla 24.1 aparece un sumario de los diferentes términos de con-

centración y de las interrelaciones correspondientes a un sistema binario que
contenga Ias especies A yB.
A menudo la composición del aire da en función delas dos especies principales que
forman parte dela mezcla gaseosa:
oxigeno, O,, yo, = 0.2 1

nitrogeno, N2, yN2= 0.79
Determínese la fracción de masa, tanto de oxígeno como de nitrógeno y el peso

molecular medio del aire, cuando se le mantiene a 25"C(298 K ) y a 1 a t m (1.O13X1O5 Pa).
El peso molecular del oxígeno es de 0.032 kg/mol y el del nitrógeno es de 0.028 kg/mol.
Como base para nuestros cálculos, tómese 1 mol de mezcla gaseosa;
oxígeno presente= (1 mo1)(0.21) = 0.21 mol

(0.032 kg)
= (0.21 mol) = 0.00672 kg
mol
538 Fundamentos dela transferencia demasa
TABLA 24.1 CONCENTRACIONES EN UNA MEZCLA BINARIADE A Y B
Concentración de masa
p = Densidad total demasa de la mezcla

PA de masa de la especie A
= Densidad
PB =Densidad de masa de la especie B
U, = Fracción de masa de la especie A = p a / p
wg = Fracción de masa de la especie B = P B / P

P=PA+PB
1=6J,+w,
Concentraciones molares
Rlezcla líquida o sólida Mezcla gaseosa
c=Densidad molar dela mezcla= n / V c=n/V=P/RT

c, =Densidad molar de la especie A = n,/ v C, = n,/ v= p , / R T
cB = Densidad molar de la especieB= nB/ v Cn = nB/ V = p B / R T
X , =Fraccibnmolar de la especie A YA = cA/c = n , / n = p , / P
= cA/c= n,/n
XB =FracciGn
molar de la especie B yn = cB/c = n B / n= p n / P
= cB/c = nn/n
c = c, + cn
1 =x, +x,
Interrelaciones
(24-10)
nitrógeno presente = (1 mo1)(0.79)= 0.79 mol
= (0.7Y mol) kg) = 0.0221 kg

mol
+
masa total presente = 0.00672 0.022 1 = 0.0288 kg
0.00672 kg
6 J 0 2 = = 0.23
0.0288 kg
0.0221 kg
6Jh: = = 0.77
()..O288kg
Como un mol de la mezcla gaseosa tiene una masa de 0.0288 kg, el peso molecular
medio del aire, debe ser de 0.0288. Cuando se toman en cuenta los otros constituyentes
que forman parte del aire, el peso molecular medio del aire se toma, a menudo, en núme-
ros redondos, como 0.029 kglmol.
Este problema podría, también resolverse utilizando la ley de los gases ideales PV=
nRT. En condicionesideales, O°C o 273 K y 1 atm o 1.O13X1O5 Pa de presión la constan-
te del gas será:
P V (1.013 X 105)(22.4m')
R=-= P a . m'
= 8.314 (24-12)
mol . K
~
nT (1 mo1)(273 K)
El volumen de la mezcla gaseosa a 298K3,es
V="=nRT
P
¡
(1 mol) 8.314- Pa ' m3)(298 K)
mo1.K
1.013x lo5Pa
= 0.0245 m'
La concentraciones son:
-
0.21 mol mol O,
= 8.59 -
'O2 - 0.0245 m' m3
- 0.79 mol = 32.3 mol N,

-
CN2 - 0.0245 m3 m'
c=C=8.59+32.3=40.9mol/m3
La densidad total, p , es:

0.0288 kg
= 1.180 kg/m3
= 0.0245 m'
y el peso molecular medio de la mezcla es:
M = -p= 1.180 kg/m' = 0.0288 kg/mol

c 40.9 mo1/m3
velocidades. Enunsistemadecomponentesmúltiples, las diferentes

especies se moverán de manera normal a diferentes velocidades; por lo tanto,
para evaluar la velocidad de la mezcla de gases, se necesitan promediar las ve-
locidades de cada una de las especies presentes.
La velocidad promedio o media de lamasa correspondiente a una mez-
cla de componetes múltiples se define en funcibn de las densidades y veloci-
dades de la masa, de todas las componentes, en la forma:
(24-13)
i= I
donde vi denota la velocidad absoluta de la especie i con relación a ejes esta-

cionarios de coordenadas.
540 Fundamentos de la transferencia de masa
Esta es l a velocidad que mediría un tubo de Pitot y la que se encontró

anteriormente en las ecuaciones de transferencia de momento. La velocidad
molar media o promedio de una mezcla de componentes múltiples, se define
en función de las concentraciones molares de todas las componentes, por me-
dio de la expresi6n:
(24- 14)
La velocidad de una especie particular con relación a la masa promedio

o velocidad molar media se llama velocidad de difusión. Se pueden definir
dos velocidades diferentes de difusión:
ui-u, que es la velocidad de difusihn de l a especie i con relación a la velo-

cidad media de l a masa.
vi-V. que es l a velocidad de difusión de l a especie i con relación a l a velo-

cidad molar media.
De acuerdo con l a ley de Fick, una especie puede tener una velocidadrelativa
a la masa o velocidad molar media, solamente si existen los gradientes en la
concentración.
Flujos. El flujo de masa ( o molar) de una especie dada es una cantidad

vectorial que denota l a cantidad de l a especie particular, ya sea en unidades
de masa o molares, que pasa en un incremento dado de tiempo atravPs de un
área unitaria normal al vector.
El flujo se puede definir refiriéndose a las coordenadas que permanecen
fijas en el espacio, a las coordenadas que se están moviendo con la velocidad
promedio de l a masa o con la velocidad molar promedio.
La relacibn básica correspondiente a l a difusión molecular define el flu-
jo molar relativo a l a velocidad molar media, la. l"ick* fue quien primero pos-
tuló una relación empírica para este flujo molar y , por lo tanto, se le llama
primera ley de Fick. Esta define l a componente A de la difusiGn en un siste-
ma isotdrmico e isohririco. Si la difusibn se lleva a cabo únicamente en la direc-
cibn de z, la ecuacihn de I:ick de la rapidez es:
(23- 15)
donde , / A , G es el flujo molar en la direccih de z relativa a l a velocidad molar

promedio, d c A /dr es el gradiente de la concentración en la direccihn de L y
*A. b'ick, Ann. Physik, 94, 59 (1855)
Transferencia demasa molecular 541
DA E , el factor de proporcionalidad, es la difusividad d e la masa o coeficiente

de difusión correspondiente a una componenteA quese difunde a través de la
componente B.
Groot* propuso una relación más general de flujo que no está restringi-
da a sistemas isotérmicoso isobáricos.
Groot escribió:
concentración difusión
) ( gradiente de )
o
Como la concentración total, c es constante bajo condiciones isotérmicas

e isobáricas, la ecuación (24-15) es una forma especial de la relación (24-16),
que es más general. Una expresión equivalente, que corresponde ajA .z,quees
el flujo d e masa en la dirección de z , relativo a la velocidad promedio de la
masa, es:
(24-17)
donde d w , / d z es el gradiente de concentración en funcihn de la fracción de

masa.
Cuando la densidad es constante, esta relación se simplifica, quedando
así:
En un sistema binario con una velocidad media constante en la direc-

ción de z , el flujo molar en la dirección de z , relativo a la velocidad molar
media también se puede expresar de la manera siguiente:
JA,z = C A (VA.=- Vz) (24-18)
Si se igualan las expresiones (24-16) y (24-18 ) , se obtiene:

542 Fundamentos de latransferencia demasa
Se puede evaluar V, para este sistema binario, por medio dela ecuación
(24-14), en la forma:
AI sustituir esta expresión enla relación que se tenía, se obtiene:
Como las velocidades componentes, V A y U B , ~ son , velocidadesrelati-

vas al eje fijo z , las cantidades C A , V A ,z y c~ U B , ~son flujos de las compo-
nentesA y B con relación a l eje fijo de coordenadas, z . Así se simboliza este
tipo de flujo, relativo a un conjunto deejes estacionarios, por medio de:
AI sustituir estos símbolos en la ecuación (24-19), se obtiene una rela-

ción que corresponde al flujo dela componente A , relativa al eje z .
(24-20)
Esta relación se puede generalizar y escribir en forma vectorial de la ma-

nera siguiente:
N A =-cDAB V Y A + Y A (NA+ NB) (24-21)
Es importante notar que el flujo molar NA es la resultante de las dos

cantidades vectoriales:
DAB VYA molar, JA , que

El flujo resulta delgradiente
de
la
concentración. Este término se llama contribución .
del gradiente dela concentración.
Y
Y A ( N A +N,)=,, V El flujomolarqueresultacuando la componente
molar A circula con el flujo global. Este término
delflujo se llama contribucibndelmovimiento
global.
Cualquiera de estas cantidades puede ser una parte importante del flu-
j o molar total,N, . Cuando la ecuación (24-21) se usa para describir una difu-
sión molar, la naturaleza vectorial de los flujos individuales N, y N,, se debe
analizar y después debe evaluarse l a dirección de cada una de las dos cantida-
des vectoriales.
Si la especie A se estuviera difundiendo en una mezcla de componentes
múltiples, la expresión equivalente a la ecuación (24-21), sería:
donde D AM es el coeficiente de difusión de A de la mezcla.

El flujo de masa, n, , relativo a un sistema fijo de coordenadas espacia-
les, se define, para un sistema binario, en función de la densidad de masa y
d i la fracción de masa, por medio de:
n A = - p D A B V wA + w A ( n+
An R ) (24-22)
donde
n~ =PAV A
En condiciones isotérmicas e isobáricas, esta relaciónse reduce a:
Las cuatro ecuaciones que definen los flujos JA ,j,, N, y nA son enun-
ciados equivalentes de la ecuación de Fick de la rapidez.
El coeficiente de difusión DAB ,es idéntico en todas.las ecuaciones, cual-
quiera de las cuales es adecuada para describir la difusión molecular; sin em-
bargo, ciertos flujos son más fáciles de utilizar en casos específicos. Los flujos
de masa, n, y j, se usan cuando también se requiere que las ecuaciones de
Navier Stokes describan el proceso. Ya que las reacciones químicas se descri-
ben en función de los moles de los reactivos que participan,los flujos molares,
JA y NA se usan para describir operaciones de transferencia de masa en las
que hay reacciones químicas. Los flujos relativos a las cordenadas fijas en el
espacio, nA y N, , se usan a menudo para describir operaciones ingenieriles
que se llevan a cabo dentro del equipo de procesamiento. Los flujos J n y j,
se utilizan para describir la transferencia de masa en celdas de difusión em-
pleadas para medir el coeficiente de difusión.
La Tabla 24.2 resume las formas equivalentes de la ecuación de Fick de
__ la rapidez
'TABLA 24.2 FORMAS EQUIVALENTES DE LA ECUACION DE FLUJO

DE MASA CORRESPONDIENTE AL SISTEMA A Y B.
Flujo
Gradiente
Ecuacihn
de
rapidez
de
Fick
Restricciones
JA 'YA
VCA
De acuerdo con la segunda ley de la termodinámica, los sistemas que nose

encuentran en equilibrio tienden a éste después de algún tiempo. Existe una
fuerza motriz generalizada, que, en términos de química termodinámica es
-dpc/dz donde ,uc es el potencial químico. La velocidad de difusión molar de
la componente A se define en función del potencial químico por medio de:
(24-23)
donde uA es la "movilidad" de la componente A o la velocidad resultante de

la molécula mientras se encuentra bajo la influencia de una fuerza motriz
unitaria.
La ecuación (24-23) se conoce como la relación de Nerst-Einstein. El
flujo molar deA se transforma en:
(24-24)
La ecuación (24-24) se puede usar para definir todos los fenómenos de

transferencia de masa molecular. Como ejemplo, estúdiense las condiciones es-
pecificadas para la ecuación (24-15); el potencial químico de una componen-
te en una solución ideal homogénea a temperatura y presión constantes, se
define por medio de:
@, = E*.'+R T In cA (24-25)
dende p o es una constante: el potencial químico del estado estándar. Cuando

se sustituye esta relación en la ecuación (24-24), se obtiene la ecuación de
Fick de la rapidez correspondiente a una fase homogénea:
(24-15)
Hay otrascondicionesfísicasademás de las diferenciasdeconcentra-

ción,queproducenungradiente de potencialquímico: las diferenciasde
temperatura y de presión así como las diferencias en las fuerzas creadas por
campos externos, tales como el gravitacional, el magnético y el eléctrico. Por
ejemplo, se puede obtener la transferencia de masa aplicando un gradiente de
temperatura a un sistema de componentes múltiples;este fenómeno de trans-
ferencia, el efecto Soret de difusión térmica, aunque usualmente pequeño en
relación con otros efectos de difusión, se usa con éxito en la separación de
isótopos. Las componentes de una mezcla de líquidos se pueden separar me-
diente una fuerza centrífuga pordifusión de presión..
Existen muchos ejemplos bien conocidos en los cuales los flujos de ma-
sa se inducen en una mezcla sujeta a un campo externo de fuerzas: la separa-
ciónporsedimentaciónbajo la influencia de la gravedad,laprecipitación
electrolítica debida a un campo electrostático de fuerzas y la separación mag-
nética de mezclas minerales a través de la acción de un campo magnético de
fuerzas. Aunque estos fenómenos de transferencia de masa son importantes,
son procesos muy específicos.
La transferenciademasamolecularqueresulta de las diferencias de
concentración y a la cual describe la ley de Fick, es el resultado del movi-
miento molecular fortuito en pequeñas trayectorias medias libres, indepen-
dientes de las paredes del recipiente. La difusión de los neutrones rápidos y
de las moléculas en poros pequefiísimoso a una densidad muy baja degas no
pueden describirse mediante esta relación.
Los neutrones producidos en un proceso de fisión nuclear poseen, ini-
cialmente, energías cinéticas altas yse llaman neutrones rúpidos a causa de SUS
altas velocidades, hasta de 15 millones de metros por segundo. Con altas ve-
locidades, los neutrones pasan a través de las cubiertas electrónicas de otros
átomos o moléculas con muy poco problema. Para deflectarse, los neutrones
rápidosdebenhacercolisiónconunnúcleo,elcualrepresentaunblanco
muy pequeño comparado con el volumen de la mayoría de los átomos y mo-
léculas. La trayectoria media libre de los neutrones rápidos es aproximada-
mente un millón de veces mayor que las trayectorias medias libres de 10s ga-
ses a presiones ordinarias. Después de que los neutrones rápidos reducen su
velocidad por medio de colisiones elásticas de dispersión entre 10s neutrones
Y el núcleo del moderador del reactor. Estos neutrones más lentos, quese Ila-
man neutrones térmicos, emigran de las posiciones de alta concentración a
posiciones de más baja concentración y la ley de Fick de la difusión describe
este movimiento.
Si la cantidad de gas es baja, o si los poros a través de los que viaja el gas
son muy pequeños, las moléculas chocarán con las paredes con mayor fre-
cuencia que entre sí. Esto se conoce con el nombre de flujo de Knudsen o di-
fusión d e Knudsen. Despuks dechocarcontralapared, las moléculas son
absorbidas, momentáneamente y después expelidas en direcciones fortuitas.
El flujo de gas se reduce debido a las colisiones con la pared. Se pueden en-
contrar relaciones para describir la difusión de Knudsen en los gases median-
te el uso de la teoría cinética. La constante de proporcionalidad que relaciona
el flujo de Knudsen con el gradiente de concentración es independiente de la
presión. Se demostrará que la constante de proporcionalidad correspondiente
a la difusión de gases, tal como la describe la ley de Fick, es inversamente
proporcional a la presión.
La transferencia de masa descrita por la ley de Fick, que resulta de las
diferencias de concentración es el proceso principal con el que se topan los
ingenieros así como el fenómeno de transferencia de masa que se estudiará
en capítulos siguientes. Es importante, sin embargo, percatarse de que estos
otros fenómenos relacionados de transporte molecular existen y de que pue-
den ser importantesen las futurasoperacionesdetransferenciademasa.
24.2 E L C O E F I C I E N T E D E D l F U S l O N
La proporcionalidad de la ley de Fick, D AB , se conoce con el nombre

de coeficiente de difusión. Sus dimensiones fundamentales, que pueden obte-
nerse a partir de la ecuación(24-15),
sonidénticasa las dimensiones fundamentales de las otras propiedades de

transferencia: la viscosidad cinemática, v y la difusividad térmica, 01, o su ra-
zón equivalente, K / p c p . La difusividad de la masa se ha dado en cm2/seg., las
unidades SI son m2/seg, o sea un factor IOp4 veces menor. En el sistema in-
glés, ft2 /h, son las unidades utilizadas. Para hacer conversiones de uno a otro
de estos sistemas se utilizan las siguientes relaciones:
(24-26)
El coeficiente de difusión depende de la presión, de la temperatura y de

la composición del sistema.En el Apéndice de Tablas, Tablas J 1,J 2 y 53,pue-
den consultase los valores experimentales correspondientes a las difusivida-
El coeficiente de difusión 547
des de los gases líquidos y sólidos, respectivamente. Como es de esperarse, de

acuerdo con la movilidad de las moléculas, los coeficientes de difusión son
generalmente mayores en relación con los gases (entre los valores de 5X lop6
y 1X m' /seg.), queenrelacióncon los líquidos (entre losvalores: 10"'
y lo-' m' /seg.), que son mayores a los valores obtenidos en relación con los
sólidos (entre y .lo"' m' /seg.).
En ausencia de datos experimentales, se han obtenido expresionessemi-
teóricas que aportan aproximaciones cuya validez es tan buena como la de
los valores experimentales debido a las dificultades que existen para la medi-
ción de estas últimas.
DIFUSIVIDAD DE LA MASA GASEOSA
Las expresiones téoricas correspondientes al coeficiente de difusión en

las mezclas gaseosas de baja densidad, en función de las propiedades molecu-
laresdelsistema, fueronobtenidasporJeans",Chapman?ySutherlandS
usando la teoría cinética de los gases.Si se utiliza el razonamiento de estos
científicos, para explicar los fenómenos de transferencia, se puede examinar
el movimiento de las moléculas de gas, tal como se hizo en las secciones 7.3 y
15.2, y después obtener una expresión que relacione el coeficiente de difu-
sión con las propiedades de los sistemas gaseosos.
Véase ahora el volumen de control de la figura 24.2. Si se especifica
que el gas es estático o que fluye en forma laminar en la dirección de x, se
puede considerar que la transferencia de masa de la especie A en la dirección
de y se lleva a cabo solamente aescala molecular. Si se aplica la ecuación(4-1),
'I
p = p (y)
A/
F i g u r a 24.2 Movimiento molecular en la superficie de un volumen de control.
*Sir James Jeans, Dynamical Theory of Gases, Cambridge Univ Press, Londres, 1921.
f S. Chapman y T. G. Cowling, Mathematical Theory of Non-Uniform Gases, Cambridge Univ. Press.,
Londres. 1959.
$ W. Sutherland, Phil. Mag. 36, 507 (1893), 38, 1 (1894).
a un flujo constante demasa a través de la cara superior del elemento,
Esta ecuacibn establece sencillamente, que el flujo ascendente de masa

debe ser igual al flujo descendente.
Como primera aproximación estudiemos un sistema que contenga molé-
culas de igual tamaño y masa, que posea velocidades mediasiguales. Sólo una
mezcla gaseosa formada por isótopos del mismo elemento se parecería apro-
ximadamente a este sistema. Si se vuelve a examinar la ecuación obtenida so-
bre una vase microscópica, se puede concluir que el número de moIéculas que
atraviesa la cara superior desde abajo deber ser igual al número de moléculas
que atraviesa esta cara desde arriba. Ya que existe una concentración dela es-
pecie A , tal como aparece en la figura 24.2, se transportarán más moléculas a
travgs de la superficie de control desde arriba, que desde abajo. Esto da como
resultado un flujo neto de moléculas deA en la dirección de y .
Tal como se hizo anteriormente en los Capítulos 7 y 15, se usarán las si-
guientes ecuaciones, obtenidas de la teoría cinética de los gases de baja densi-
dad:
donde? es la velocidad molecular fortuita media,¿?/4 es la velocidad de una

molécula individual mientras pasa a través del área A x Az, A es el paso me-
dio libre, K es la constante de Boltzman, m es la masa de una molécula, d es
el diámetro de las moléculas esféricas, N es la concentra.ción molecular y Z es
la frecuencia con la que llegará una molécula al área Ax Az. La ecuacibn de
continuidad, descrita en funcibn delas moléculas en movimiento, es:
c c
m,- Ax Azl,--C m,- Ax Azly+= O
n 4 n 4
o, al sumar las N moléculas en un volumen unitario:

Si van a tomarse en cuenta solamente las mcdéculas de A que atraviesan

esta superficie, la ecuación debe ser válida para un flujo neto de masa en la
dirección de y
(24-27)
De nuevo, como en los Capítulos 7 y 15, se supondrá que el perfil de la

concentración es esencialmente lineal en una distancia de varias trayectorias
medias libres. Entonces:
donde
y-=y-6
Y
y+=y+6
La sustitución de estas relaciones, correspondiente apAly- andpA(,+ en la

ecuación (24-27), da como resultado:
(24-28)
Ahora, nuestra ecuación se transforma en:
I - aPA
jA,y = -3c.4- (24-29)
ay
Si se compara la ecuación (24-29) conla (24-17),
(24-17)
se pone de manifiesto que el Coeficiente de difusióncorrespondientea la

mezcla de moléculas semejantes, o sea, A y el isótopo A * , es:
1 - (2430)
DAA* =3 C
A
Este coeficiente de difusión se denomina, a menudo, coeficiente de au-
todifusibn y se le utiliza para explicar la difusión de las moléculas trazadoras.
La sustitución de los resultados cinéticos correspondientes a y A, en la c
ecuación (24-30), da:
(24-31)
En un gas ideal se puede reemplazar a N , usando la relación
NKT=cRT= P
y obteniéndose:
La ecuación (24-32) demuestra que el coeficiente de difusión se puede

expresar totalmente en función delas propiedades del gas. De manera distinta
a lo que ocurre en otros coeficientes de transferencia de dos moléculas,l a vis-
cosidad y la conductividad térmica, el coeficiente de difusión depende de la
presión así como de un orden más alto de temperatura absoluta. El significa-
do de esto no debe pasarse por alto aunquese haya usado en el desarrollo un
modelo muy simplificado.
Las versiones modernas de la teoría cinética han intentado tomar en
cuenta las fuerzas deatracción y repulsiónexistentesentre las moléculas.
Nirschfelder, Bird y Spotz", utilizando el potencial de Lennard Jones para
evaluar la influencia de las fuerzas intermoleculares, encontraron una ecua-
ción adecuada al coeficiente de difusión correspondiente a parejas gaseosas
de moléculas no polares? y no reactivas.
(24-33)
/"donde DA, es la difusividad de la masa de A , que se difunde a través de B ,

,/ en /seg., T es latemperaturaabsolutaen K;M,, M, son los pesos mo-
leculares de A y B , respectivamente; P es la presión absoluta en atmósferas;
*J. O. Hirschfelder, R. B. Bird y E. L Spotz, Chemm. Revs. 44, 205-231 (1949).
Para un estudio introductorio de estructuras polares y no polares, colsúltense F. Daniels y R. A. Al-
berty, Physical Chemistry, Wiley, Nueva York, 1955.
IEl coeficiente de difusión 551
u A B es el “diámetro de colisión”, que es un parálmetro de Lennard-Jones en

Angstroms; R, es la “integral de colisión” correspondiente a la difusión mo-
lecular, que es una función adimensional de la temperatura y del campo po-
tencial intermolecular correspondiente a una molécula de A y una de B. En
la Tabla K 1 del apéndice, aparece una lista de valores de a, en función de
; k es la constante de Boltzmann,igual a 1.38 X 10” ergios /I( y E
es la energía de la interacción molecular que corresponde al sistema binarlo
4B
AB,un parámetro de Lennard-Jones en ergios; ver ecuación (24-40). La ecua-
ción (24-33) también señala la dependencia de la! presión y la temperatura,
tal como se demostró en la ecuación (24-32).Como la ecuación (24-33) defi-
ne la propiedad de transferencia de una mezcla gaseosa binaria, no es sor-
prendente encontrar que el coeficiente de difusión también depende de las
componentes que forman el sistema. Cuando se examinó el proceso de trans-
ferenciaenuna fase deuna sola componente, 110 se encontró ninguade-
pendenciade lacomposiciónen la ecuación(2.4-32) ni en las ecuaciones
similares correspondientesa la viscosidad y a laconductividadtérmica.
Los parámetrosdeLennard-Jones, u y , se obtienenusualmentede
datos acerca de la viscosidad. Desafortunadamente esta información puede
obtenerse solamente para unos cuantos gases puros. La Tabla K2 Apéndice
tabula estos valores. En la ausencia de datos experimentales,los valores de las
componentes puras se pueden calcular a partir ‘de las siguientesrelaciones
empíricas :
u = 1.18 vp3 (24-34)

u = 0.841 Vc1’3 (24-35)
(24-36)
EA/K = 0.77 T, (24-37)
EA/K = 1.15Tb (24-38)
donde l/’b es el volumenmolecualren el punto normaldeebullición,en

(cm)3 /g mol (esto se calcula utilizando la Tabla 214.3); V, es el volumen mo-
/ g T, es la temperaturacríticaen
lecular críticoen ( ~ m ) ~mol; K; es la
temperatura normal de ebullición en K y P, es a l presión crítica, en atm-ós-
feras.
Para un sistema binario compuesto de pareja:s moleculares no polares, se
pueden combinar los parámetros de Lennar-Jones de la componente pura,
emperíricamente, por medio de las siguientes relaciones:
(24-39)
Y
EAR = JEAER (24-40)
Estas relaciones deberán modificarse para parejas moleculares: polar-po-

lar y polar-no polar. Hirschfelder, Curtiss y Bierd* analizaron las modifica-
ciones propuestas.
La ecuación de Hirschfelder, Bierd y Spotz, ecuación (24-33), parece
ser la mejor ecuación de que se dispone actualmente para el cálculo de las di-
fusividades demasa quecorrespondena los sistemasgaseosos no polaresbi- ' 8
narios. Reid y Sherwood**, comparando los valores experimentales con los

valores calculados para 80 parejas de gases, encontraron que los valores predi-
chos por las ecuaciones eran correctos, con un error máximo de un 6% con
relación a los valores medidos.
Como los coeficientes de difusión son difíciles de medir con exactitud,
esta desviación se puede deber, parcialmente, a los errores experimentales.
Otras ecuaciones empíricas sugeridas por Arnoldt y Gillilandt y por Slattery
y Bird 5 , no concuerdan tanto.
De acuerdo conlo anterior, se recomienda utilizar la ecuación deHirsch-
felder cuando no pueda hacerse uso de los valores experimentales.
La ecuación de Hirschfelder, que es la (24-33), se utiliza a menudo para
extrapolar los datos experimentales. Para valores moderados de la presión,
hasta de 25 atmósferas, el coeficiente de difusión varía inversamente a la pre-
sión. Las presiones más altas requieren, aparentemente, correcciones corres-
pondientes alos gases, densos; desafortunadamente no existe ninguna relación
satisfactoria para las presiones altas. La ecuacion (24-33) también establece
que el coeficientededifusiónvaríacon la temperatura,deacuerdocon
l a relación: ? I 2 laD.Si se simplifica la ecuación (24-33), se puede predecir
el coeficiente de difusión a cualquier temperatura y a cualquier presión me-
nor de 25 atmósferas, a partir de un valor experimental conocido, por medio
de :
(24-41)
En la Tabla J 1 del Apéndice, los valores experimentales del producto

DA P aparecen en una lista que corresponde a varias parejas de gases a una
temperatura particular. Cuando se usa la eucación (24-41) se pueden exten-
der estos valores a otras temperaturas.
*J. O. Hirschfelder, C. F. Curtiss y R. B. Bird, Molecular Theory of Cases and Liquids, Wiley &Sons,
Inc., Nueva York 1954.
**R.C. Reid y T. K. Shenvood, The Properties o f Gases and Liquids, McGraw-Hill Book Company,
Nueva York 1954. Cap. 8.
tJ. H. Arnold, J. Ann Chem. Soc. 52, 3937 (1930)
$E. R. Gilliland, in Eng. Ghem. 26, 681 (1934)
§J. C. Slattery y R. B. Bird, A. I. Ch. E. J. 4, 137 (1958)
Calcúlese el coeficiente de difusión del bióxido de carbono en aire a 20°C y presión

atmosférica. Compárese este valor con el valor experimental que aparece en la Tabla J 1
del Apéndice. Los valores de (T y E/K se obtienen de la Tabla. K2 del Apéndice.
u, inA EA/K, in K
Bióxido de carbono 3.996 190
aire 3.167 97
Los diversos parámetros que se necesitan para resolver la ecuación (24.33) se pue-
den calcular de la siguiente manera:
~~+u~-3.996+
- 3 . 6 1 7 - ~ . ~ ~ ~ ~
UAB=--
2 2
”
EA,/. =J(EA/K)(EB/K) J(190)(97) = 136

T = 20+273 = 293 K
P = 1 atm
KT
”
- 2.16
€AB
S
Z
,(TablaK.l) = 1.O47
M,, = 44
Al sustituir estos valores en la ecuación, se obtiene:
+ i/A4rB)”2
o.oO~8~8T3”(1/MA
DAB
PuAB2fin,
- - O,,47 cm’/s
- (0.001858)(293)”2(1/44+ 1/29)”2
”
(1)(3.806)2(1.047)
De la Tabla J1 para CO en aire a 2 7 3 K , a una atmósfera, $;etiene:
2
DAB= O. 136 cm‘/s

La ecuación (24-41) se utilizará para corregir las diferencias de temperatura,
DAB,,,
Los valores de 2
!, se pueden calcular así:
136
at T2= 273 T-= 0.498
e A e / ~=
273
.&ITz = 1.074
at TI= 293 R,lT, = 1,047 (cálculo anterior)

El valor corregido correspondiente a l coeficiente de difusión a 20cC es:
(0.136)=0.155cmZ/s ( 1 . 5 5 ~ 1 0 - ~ m ~ / s )
Se ve claramente que la dependiencia de la “colisiGn integral” de la tem-

peratura, es muy pequeña. Por esto, la mayoría de los valores de las difusivi-
dades con relación a la temperatura sólo incluyen la razón ( T I / T),3 ’ 2 .
La transferencia de masa en las mezclas gaseosas de varias componentes
puede describirse por medio de ecuaciones teóricas que incluyen los coefi-
cientes de difusión Correspondientes a las diversas parejas binarias que for-
man la mezcla. Hirschfelder, Curtiss y Bird*, elaboraron una expresión en su
forma más general. Wilkei simplificó la teoría y demostró que una buena
aproximación a la formacorrecta, es la dada por la relación:
donde D,- m e r c l a es la difusividad del a masa que corresponde ala componen-

te 1 en la mezcla gaseosa; L), - n es la difusividad de la masa para el par binario,
la componente 1 que se está difundiendo en la componente n ; y y,‘ es l a frac-
ción molar de la componente n en la mezcla gaseosa calculada con base libre
en la componente 1, esto es:
Y2
Y2’ =
Y2+Y3+. . .+Yfl
EJEMPLO 3
Determínese la difusividad del monóxido de carbono en una mezcla de oxígeno en

la cual las fracciones molares de cada una de las componentes, son:
yoz = 0.20
yN2= 0.70
yco = 0. 10
La mezcla gaseosa está a 298 K y a 2 atmósferas de presión total.

En la Tabla J 1 del apéndice se encuentra:
I \’ ,-!’
D,,,, = 0.185 x mZ/s a 273 K, 1 atm‘
,.h
D,,..,, = 0.192 X m2/s a 288 K, 1 atm
* J. A. Hirschfelder, C. F. Cirtiss y R. B. Bird, ~VfolecularTheory of Gases and Liquids, Wiley, Nueva

York, p. 7 18.
? C . R. Wilke, Chem. Engr. Prog., 46, 95-104 (1950).
El coeficientededifusión 555
Los dos coeficientes binarios de difusión se pueden corregir. para diferencias de presión y
temperatura utlilizando la ecuación (24-41),
D A B condición 1 /
I
DA, condición 2
<\ 52\
Para 298 K y 2 atm, se tiene:
Las composiciones del oxígeno del nitrógeno con base en el CO-libre .son:
o.2o -0.22
yo2' = ___ -
1-0.10
0.70
yNI' = -- - 0.78
1-0.10
Y sustituyendo estos valores en la ecuación (24-42) se obtiene
D c o - o ~ ,-
-N ~ 1
0.22 0.78
0.105 X 1 0 - ~
' 0.101X 1 0 - ~
= 0.102 x loe4mz/s (0.395 ft2/hr)
DIFUSIVIDAD DE LA MASA LIQUIDA
En contraste con los gases, para los cualesexisteuna teoría cinética

avanzada para explicar el movimiento molecular,las teorías de que se dispone
para explicar la estructura de los líquidos y sus características de transferen-
cia aún son inadecuadas para permitir un tratamiento riguroso.La inspección
de los valores experimentales publicados en el Apéndice 52,relativos a los
coeficientes de difusión de los líquidos, revelan que existen diversos órdenes
de magnitud menores que los coeficientes de difusión de los gases y que depen-
den de l a concentración debida a los cambios de viscosidad con l a concentra-
ción así como a modificaciones en el grado en que la solución sea ideal.
Ciertas moléculas se difunden como moléculas, en tanto que otras, que
se llaman electrólitos, se ionizan en las soluciones y se difunden como iones.
Po: ejemplo, el cloruro de sodio, N a C1, se difunde en agua como los iones de
N a y C1-. Aunque cada uno de los iones posee distinta movilidad, la neutra-
lidad eléctr~ica de la solución indica que los iones deben difundirse en alguna
proporción; por lo tanto, es posible hablar de un coeficiente de difusión de

10s electrólitos tales como el DjaC1. Sin embargo, si se encuentran presentes
varios iones, se debe tener en cuenta la rapidez de difusión de los cationes y
aniones individuales y los coeficientes de difusión molecular no tienen signifi-
cado. No es necesario decir que se necesitarán correlaciones separadas que se
utilizan en la predicción de l a relación entre las difusividades de la masa lí-
quida, para los electrcilitos, así corno para 10s no electr6litos.
Se han postulado dos teorías, l a del “hoyo” de Eyring y la hidrodinámi-
ca, como posibles explicaciones al fenómeno de difusión de solutos no elec-
trolíticos en soluciones de baja concentración. En el concepto de Eyring el
líquido ideal se trata como un modelo de capas casi cristalinas interespacia-
das con hoyos. Entonces se describe el fenómeno de transferencia por medio
de un proceso que incluyeel salto de las moléculas de soluto en los hoyos del
modelo de capas. Estos saltos se relacionan empíricamente con la teoría de
Eyring de rapidez de reacción*. La teoría hidrodinámica establece que el coefi-
ciente de difusión de los líquidos se relaciona con la movilidad de las molécu-
las de soluto; esto es, con l a velocidad neta de la molécula mientras ésta se
encuentra bajo l a influencia de una fuerza motriz unitaria. Las leyes de l a hi-
drodinámica establecen relaciones entre l a fuerza y l a velocidad. Una de las
ecuaciones desarrolladas a partir de la teoría hidrodinámica es l a ecuación de
Stokes-Einstein:
donde D A B es l a difusividad de A en una solución diluida de B,k es la conso-

nante de Boltzmann, T es la temperatura absoluta,7 es el radio de l a partícula
de soluto y pB es l a viscosidad del solvente. Esta ecuación ha sido bastante
útil en l a descripción de la difusibn de las patículas coloidales de moléculas
redondas y grandes, en un solvente que se comporta como un continuo enre-
lación con l a especie en difusión.
Los resultados de ambas teoríasse pueden reordenar en la forma general:
(24-44)
en la cual f(u) es una función del volumen molecular del soluto en difusión.
Se han obtenido correlaciones empíricas en un intento de predecir el coefi-
ciente de d i f u s i h de los líquidos en función de las propiedades del soluto y
del solvente. Wilke y Changt propusieron l a siguiente correlación, que es la
mejor de que se dispone para los no electrólitos en una solución diluida
* S. Glasstone, K. J. Laider y H. Eyring, Theory o f Rate Processes, McGraw-Hill Book Company, Nue-
va York, 1941, Capítulo IX.
?C. R. Wilke y P. Chang, A. I. Ch E. J., 1, 264 (1955).
donde D A B es la difusividad de la masa de A que se difunde en el solvente lí-

quido B en cm2/seg, pB es la viscosidad de la solución, en centipoises; T es la
temperatura absoluta, en K, M B es el peso molecular del solvente; v6 esel
volumen molar del soluto en el punto normal de eb'ullición, en cm3/g mol y
aBes el parámetro de "asociación" correspondiente al solvente B.
Los volÚmenes mokculares a los puntos normalesdeebullición, v6 ,
aparecen tabulados en la Tabla 24.3. Para otros compuestos, los volúmenes
atómicos de cada uno de los elementos presentes se suman mediantelas fórmu-
lasmoleculares.EnlaTabla24.4aparece un lista'dode las contribuciones
correspodientes a cada uno de los átomos constituyentes. Cuando se inclu-
yen ciertas estructuras de anillo, deben hacerse correcciones que correspondan
a la configuración específica de anillo. Se recomiendan las siguientes correc-
ciones:
para el anillo triple, como el óxido de etileno restar 6
para el anillo cuádruple como el ciclobutano restar 8.5
para el anillo quintuple,como el furanorestar 11.5
15 para
restarla piridina
para el anillo de benceno
restar 15
para el
naftalina
restar
de
anillo 30
el de
paraantraceno
anillo
47.5
restar
TABLA 24.3 VOLUMENES MOLECUALRES A LA TEhlPERATURA DEL
PUNTO NORMAL DE EBULLICION DE ALGUNOS
COMPUESTOS COMUNES.
Volumen
ar molecular
Compuesto
en
cm3 /g mol
Compuesto
en cm3 /g mol
Hidrógeno H, 14.3 Oxido
nítrico NO 23.6
Oxígeno, O, 25.6 Oxido nitroso!NzO 36.4
Nitrógeno N, 31.2 N H,, 25.8
Aire 29.9 d;bo $ 5 18.9
Monóxido decarbono Co 30.7
Bióxido de carbono CO, 34.0
hidrógeno
Bromuro Br,
& 32.9
53.2
Sulfuro decarbonilocos 51.5 Cloruro Cb 48.4
Bioxido
sulfúrico SO, 44.8 Yodo I, 71.5
LOS valores recomendados del parimetro de asociación, aBpara unos Po-

COS solventes se proporcionan en seguida:
E.JEhlPL0 4
Calcidese el coeficiente de difusióndel líquido, del etanol, 5 €1, OH en agua a
loot. El volumen molecular del etanol se puede evaluar utilizando 10s valores de la Tabla
24.4, en la forma siguiente:
TABLA 24.4 VOLUMENES ATOMICOS CORRESPONDIENTES A

VOLUMENES MOLECULARESCOMPLEJOS DE
SUSTANCIAS SIMPLES
atómico
Volumen
atómico
Volumen
Elemento
cm3
en /g Elemento
mol en cm3 /g mol
Bromuro 27.0 Oxígeno,
excepto en la 7.4
Carbono 14.8 forma anotadaabajo
Cloro 21.6 Oxígeno, en los ésteres 9.1
Hidrógeno 3.7 metílicos
Yodo 37.0 Oxígeno en los éteres 9.9
Nitrógeno, doble 15.6 metílicos
enlace Oxígeno enéteres
los 11.0
Nitrógeno en los 10.5 altos y otros ésteres
aminos primarios Oxígeno, en los 12.0
Nitrógeno en los 12.0 ácidos
aminos secundarios Azufre 25.6
* G. Le Bas, The Molecular Volumes of Liquid Chemical Compounds, Long mans, Green y Co., Ltd.,
Londres, 1915.
Solvente @B
-~
-~~~~ ~ ~ ~ ~ ~
Agua 2.6
Metano1 1.9
1.5
Etanol
Benceno, éter, heptano y
otros solvevtes no asociados 1.o
A 10°C, la viscosidad de una solución que contiene 0.05 mol de alcohol porcada li-
tro de agua es de 1.45 centipoises; los parámetros restantes quevan a usarse, son:
F283K
correspondiente al agua=2.6
para agua=l8
si se constituyen estos valores en la ecuación (24-45),se obtiene:
7.4~
10-*(2.6x 18)”’
DC~H~OH-H~O=(59.2)’.6
= 8.5 x crn2/s (8.5 X 10”’ rn’/s)
Este valor concuerda de manera excelente con el valor experimental de 8.3x10”*

m’ /seg, que aparece en el Apéndice J.
Laecuación(24-45) se recomiendasolamenteensoluciones diluidas

de solutonodisociativos.Para estas soluciones,predice valores cuyomar-
gen de error es de 13 % con relación a los datos experimentales*. La ecua-
ción ( 2 4 - 2 5 ) n o es aplicable a la difusión de las moléculas grandes; esto es:
si V , > O .2 7 ( 4 M, ) ' 7 , entonces se recomienda utilizar una forma redu-
cida de la ecuacion de Stokes-Einstein:
DAB =
1.o5 x 10" T
Pvb'/3
(24-46) I
en el caso de las moléculas grandes de soluto.
Las propiedades de las soluciones conductoras de la electricidad se han
estudiado ya en forma intensiva durante más de 75 años. Aún así, las relacio-
nes conocidasentrelaconductividadeléctricayelcoeficiente de difusión
del líquido son válidas únicamente para soluciones de sales diluldas en agua.
El coeficiente de difusión de una sal equivalente en solución dilulda, está da-
da por la ecuaciónde Nernst:
(24-47)
donde DA, es el coeficiente de difusión basado en la concentración molecu-

lar de A ,en cm2 /seg;R es la constante universal de los gases, 8.31 6 joules/(K)
(g mol); T es la temperatura absoluta en K ; X + O , h-.O son las conductancias
iónicas límite (concentración cero) en (amp/cm2) (volt/cm) (g equivalente/
cm3 ) y es la constante de Faraday, de 96,500 coulombs/g equivalente. Esta
ecuación se ha ampliado a iones polivalentes reemplazando la constante nu-
mérica2por (l/n' !-1/n-)donde n y n- son las valencias delcatióny del
anión, respectivamente
DIFUSIVIDAD DE LA MASA SOLIDA
En cualquier estudiodel movimiento molecular en el estado sólimdo, la ex-

plicación de la transferencia de masa se divide autornáticamente en dos cam-
pos mayores de interés, la difusión de los gases o líquidos en los poros del
sólidoylainterdifusión de los constituyentes sólidos por medio del movi-
miento atómico. La primera clase de difusión tiene un papel preponderante
en la catálisis y es importante para el ingeniero qu.ímico. Los metalurgistas
son los principales investigadores de la difusión de los átomos en los sólidos.
La difusión en los poros se puede llevar a cabo por medio de tres o más
mecanismos: difusión de Fick, difusión de Knudsen y difusión superficial. Si
los poros son grandes y el gas relativamente denso, la transferencia de masa
se llevará a cabo por medio de una difusión de Fick. Dentro del catalizador,
las trayactorias de difusión son como canales de forma irregular; por lo tanto,
el flujo es menor de lo que sería en poros uniformes de la misma longitud y

radio promedio. El flujo de masa se describe en función del coeficiente "efec-
tivo" de difusión, por medio de la ecuación:
J A = - c D A .V~
Y A~ (24-48)
La magnitud del coeficiente depende de las variables que influyan en la

fase de difusión, tales como la temperatura y la presión, así como de las pro-
piedades del catalizador, tales como espacio fracciona1 vacío,6 , un factor an-
gular de longitud, L ' y un factor de forma S ' ,
(24-49)
donde 7 es la tortuosidad, o sea, un factor que describe la relación entre la

longitud real de la trayectoria y la longitud nominal del medio poroso. Satter-
field" encontró valores experimentales correspondientes a: 0 ,T y D A , ef.
La d i f u s i h d e Knudsen ocurre cuando el tamaño de los poros es del or-
den del de la trayectoria media libre dela molécula en difusibn. Se ha encon-
trado una relación que describe la difusión de Knudsen y su coeficiente de
difusión, usando la teoría cinética de los gases.
(24-5 O)
Dk,eff =
7S'p (24-5 1)
donde x. es la longitud de la trayectoria de difusibn, p es la densidad de la

partícula catalizadora y 111 el peso molecular del gas en difusión
La difusión superficial tiene lugar cuando las moléculas que se han ab-
sorbido son transportadas a lo largo de la superficie como resultado de un
gradiente bidimensional de concentración superificial. f Usualmente la difu-
sión superficial tiene un papel menor en la difusión en un sólido catalizador a
menos que exista unagran cantidad de absorción.
La introduccibn de un átomo extraño puede alterar las propiedades im-
portantes de un metal.
El edurecimiento del acero es un proceso industrial basado en la difu-
sión del carbono y otros elementos en el hierro. Se han producido semicon-
ductores por medio de la difusión de átomos "dopados" dentro de cristales.
Los cuatro mecanismos que se postulan para la difusión en los sólidos son: la
difusión de vacancia, la intersticial, la de intersticialidad y la difusión por in-
tercambio directo de moléculas adyacentes.
* C. N. Satterfield, Mass Transfer in Heterogeneus Catalysis, M. I. T. Press, Cambridge, Mass, 1970.

-k v. C . Levich, Physicochemical Hydrodynamics, Prentice Hall, Inc., Englewood Cliffs, N. J. 1962.
En todos los cristales que se encuentran en equilibrio térmico a tempe-

raturas superiores al cero absoluto, existen sitios desocupados o vacantes en
las redes. Un átomo puede saltar de su posición a una vacante vecina, tal co-
mo puede verse en la figura 24.3 a.
0 . 0 . . e....
e o 0 0 0 .*.e e 0
..e 0 e a e a .
0 . 0 . . ...e.
(a)Difusiónlugar
de vacante (b) Difusión intersticial
O.... ...e.
0 .A.
0 . 0 . .
0 . 0 . .
0
(C) Difusión de intersticialidad

A.
....e
(d) Intercambio directo
e
Figura 24.3 Difusión en un cristal.
El átomo en difusión continúa moviéndose a través del cristal por medio

de ufia serie de intercambios con sitios vacantes qu.e aparecen adyacentes a é1
de vez en cuando. Esto usualmente requiere de una distorsión de la red. Este
mecanismo se ha descrito matemáticamente suponiendo un proceso de rapidez
unimolecular y aplicandoel concepto de Eyring de “estado activado”, tal como
se estudió en la teoría de los “hoyos” correspondiente a la difusión de los lí-
quidos. La ecuación resultante es complicada y no se estudiará aquí. Rela-
ciona la difusividad con la relación geométrica entre las posiciones en la red,
la longitud de la trayectoria de salto, la frecuencia de vibración del átomo
en la dirección del salto y la energía de activación asociada al salto. La difu-
si& de vacunciu es el mecanismodominante de difusión en los metales y
aleaciomes cúbicos centrados en las caras. El mecanismo se ha usado para ex-
plicar la difusión en materiales cúbicos centrados en el cuerpo así como en
compuestos iónicos y óxidos.
Un átomo se mueve en difusibn intersticial saltando de un sitio inters-
ticial a otro vecino, como puede verse en la figura 24.3 b. Esto usualmente
incluye una dilatación o distorsión de la red; así pues, el método matemático
que incluye la teoría de Eyring de rapidez unimolecular también se usa pa-
ra definir la difusividad intersticial. Este coeficiente de difusión se define en
función de las mismas variables que en el caso de la. difusión de vacancia.
Cuando el átomo intersticial en difusión posee un tamaño comparable al
de átomos de la red, tiene lugar un proceso diferente del de difusión intersti-
cial. Si se empuja el átomo de una red vecina hacia el sitio intersticial adya-
cente, el átomo intersticial se puede mover hasta un sitio normal en l a red.

Este mecanismo, tal como aparece en la figura 24.3c, se llama difusi6n de
intersticialidad.
El cuarto mecanismo que se propone para explicar la autodifusión en
los metales y aleaciones incluye el intercambio directo de más átomos por
medio de una rotación en forma de anillo. Los tres átomos giran como en un
tiovivo y después de dos intercambios, el átomo A se corre hasta una posi-
cicjn nueva, A ’ . Después efectúa una rotación con los otros dos átomos, Ile-
gando finalmente a A ” . El mecanismo de anillo no se sabe que ocurra en
ningim metal o aleación, pero se ha sugerido como un mecanismo que puede
explicaralgunasanomalíasaparentesen los coeficientesdedifusión de los
metales con redes centradas en el cuerpo.
Existen referencias excelentes para un estudio más detallado de la trans-
ferencia de masa de los átomos a través de cristales.*
24.3 T R A N S F E R E N C I A C O N V E C T I V A D E M A S A
L a transferencia de masa debida a la convecci6n consiste enla transferen-

cia entre un fluido en movimento y una superficie o entre dos fluidos en mo-
vimiento, relativamente no miscibles.Este modo de transferencia depende,
tanto en las propiedades de transferencia como de las características dinámi-
cas del fluido que está fluyendo.
Se debe lhacer una distinción entrelos dos tipos deflujo, tal como se hizo
en el caso de la transferencia convectiva de calor. Cuando una bomba u otro
aparato semejante ocasiona. el movimiento del fluido,el proceso se llama con-
vección forzada. Si el movimiento del fluido se debe a una diferencia de den-
sidades, que puede haber surgido como resultado de una concentración o de
una diferencia de temperatura, el procesose llama convección libre o natural.
La ecuación de rapidez correspondiente a l a transferencia convectiva de
masa, generalizada en forma análoga a la “ley” de Newton del enfriamiento,
ecuación (15-1l ) , es:
NA kc ACA (24-52)
donde NA es l a transferencia de masa molar de la especie A , medida con re-

lación a coordenadas especiales fijas; AcA es la diferencia entre l a concentra-
ción de l a superficie límite y la concentración media de l a corriente de fluido
de la especie A e n difusión y K c es el coefieciente de transferencia convecti-
va de masa.
* A t o m Movement, American Society for Metals, 1951, R. M. Barrer, Diffusion in and Through Solids,
Cambridge University Press, Londres, 1941. P. G. Shewmon, Diffusion of Solids, McGraw-Hill Book
Company Nueva York, 1963. L. A. Girifalco Atomic Migration in Crystak, Blaisdell, Nueva York,
1964. B. L. Sharma, D f f w i o n in Semiconductors, Trad. Tech Publications, D-3392, Clausthal-ZeUer-
feld, Alemania, 1970.
Conclusión 563
Como se hizo en el caso de la transferencia de masa molecular, la trans-

ferencia convectiva de masa, la transferencia se lleva a cabo en la dirección de
una concentración decreciente. La ecuación (24-52) define el coeficiente Kc
en función del flujo de masa y la diferencia de concentración desde el princi-
pio hasta elfinalde la trayectoria de difusión. Por lo tanto, el coeficiente
incluye las características de las regiones de flujos laminar y turbulento del
fluido, en cualesquiera proporciones en las que se encuentren. En los capítu-
los 28 y 29 se estudiarán los métodos de determinación de este coeficiente.
En general es una función de: la geometría del sistema, las propiedades del
fluido del flujo, y la diferencia de concentración AcA .
Puede recordarse, a partir de la experiencia que se ha tenido con fluidos
que circulan alrededor de una superficie, que siempre hay una capa, a veces
extremadamente delgada, cercana a la superfice, donde el flujo es laminar y
las partículas próximas a la frontera sólida se enculentran en reposo. Como es-
t o siempre se cumple, el mecanismo de transferencia de masa entre una super-
ficie y un fluido debe incluir una transferenciade masa molecular a través de
las capas estancaday laminar de fluido. L a resistencia que controla la transfe-
renciaconvectiva de masaes,amenudo,elresultado de esta "película"
de fluido y el coeficiente Kc se llama, de acuerdo con esto, coeficiente de la
película.
Es importante para el estudiante percatarse de la estrechasimilitud entre
el coeficiente de transferencia convectiva de calor y de masa. Esto, de inme-
diato, sugiere el hecho de que las técnicas desarrolladas para la evalución del
coeficiente de transferencia convectiva de calor s e puede aplicar a la transfe-
rencia convectiva de masa. En el Capítulo 28 se hará un estudio completo de
los coeficientes de transferencia convectiva demasa así como de su evalua-
ción.
~~
24.4 C O N C L U S I O N
En este capítulo se han introducido los dos modos de transferencia de

masa, molecular y convectivo. Como la difusión de masa incluye una mezcla
múltiple de componentes, se presentaron las relaciones fundamentales corres-
pondientes da las concentraciones y velocidades (de la especie individual así
como de la mezcla. Se han estudiado la propiedadde transferencia molecular,
D A , , el coeficiente de difusión o difusividad de la masa en sistemas gaseosos,
líquidos y sólidos y se han presentado las ecuaciones que los correlacionan.
Las ecuaciones de rapidez correspondientes a la transferencia demasa
de la especie A en una mezcla binaria, son:
Transferencia molecular de masa
JA = - c D A B V Y A flujo
molar
relativo
la
aveloci-
dad molar media
564 Fundamentos de iatransferencia de masa
j, = -pDA~VoA flujo
masa
de
relativo
velo-
la a
cidad de la masa promedio
N, = -cDABVyA + yA(NA+ N B ) flujomolarrelativoacoorde-
nadas espaciales fijas
nA = - P D A B V W A + W A (nA + n B ) flujode masarelativoacoor-
denadas espaciales fijas
transferencia convectiva d e masa:
_____
496 P R O B L E M A S
24.1 Se va a llevar gas natural líquido, GNL, de la Península de Kenai, en

Alaska, a un lugar de la Bahía Yakina, en Oregon, por medio de un
transportemarítimo. La composiciónmolardelGNLcomercial es:
metano Ch, 94.9%

etano C, H, 4 .O%
propano C, H, 0.6%
bióxidodecarbono CO, 0.5%
Determinese:
(a) la fracción de peso del metano,
(b) el peso molecular medio de la mezcla de GNL
(c) la densidad de la mezcla gaseosa cuando se le calienta a 193 K y
1.O13X1O5 Pa.
(d) la presión parcial del metano cuando la presión total del sistema es
de 1.O13X1O5 Pa.
(e) la fracción de masa del propano en partes por millón.
24.2 Se analiza un gas natural, resultando: CO, ,4%;CH4, 77% , C, H,, 5.0%
y N,, 14% . Determine las siguientes propiedades de la mezcla gaseosa:
(a) fracción molar de C, I & ,
(b) fracción de peso de N,
(c) peso molecular medio de la mezcla de gas natural
(d) presión parcial de cada una de las componentes si la presión total
es de 1,500 mm Hg.
24.3 En un recipiente de 30 m3 hay aire a 400 K y 1.O13X1O5 Pa. Determi-
ne las siguientes propiedades de la mezclagaseosa:
(a) fracción molar de O,,
(b) fracción volumétrica de O,,
(c) peso de la mezcla
25
ECUACIONES DIFERENCIALES DE
LA TRANSFERENCIA DIE MASA
En el capítulo 9 se obtuvieron las ecuaciones diferenciales generales de

transferenciademomentoutilizando el conceptodevolumendiferencial
de control. En el capítulo 16 se generaron en forma ;análogalas ecuaciones di-
ferenciales generales de la transferencia de calor. Se usará de nuevo este mé-
todo para obtenerlas ecuaciones diferenciales de la transferencia de masa. Se
establecerá, también, la ecuación de continuidad de una especie dada haciendo
un balance de la masa en un volumen diferencial de control.
Se obtendrán ecuaciones diferenciales adiciona.les al incluir en la ecua-
ción de continuidad, las relaciones correspondientes al flujo de masa encon-
tradas en el capítulo anterior.
25.1 L A E C U A C I O N D I F E R E N C I A L D E
T R A N S F E R E N C I A DE M A S A
Analicemos el volumen decontrol Ax Ay A,z, a travésdelcualestá

fluyendo una mezcla, incluyendo la componente A , tal como se aprecia en
la figura 25.1. La expresión de volumen de control clue corresponde a la con-
servación de la masa, es:
lo cual puede decirse con palabras, de la siguiente manera:
Rapideznetadeflujoderapideznetadeacumulación
1masa que sale delvolumende

control de
masa dentro delvolu-
control
men de
57 1
572 Ecuaciones diferenciales de la transferencia de masa
Figura 25.1 U n volunlende control diferencial
Si se estudia la conservacih de una especie determinada, esta relación debe

incluir un término que corresponda a la produccihn o desaparición de A por
medio de una reacción química que ocurra dentro del volumen. La relación
general para el equilibrio de la masa de la especie A , correspondiente al vo-
lumen de control bajo estudio, puede enunciarseasí:
Rapidez neta de rapidez neta

de rapidez
pro-
de
de A , desde el
/flujo masa
volumende de control 1 -+
de A dentro del
[acumulación
volumen
de
control
volumen
de I[de A dentro del
ducción
control
- química
(25-1)
1
Se evaluarán los términos correspondientes a la especieA y se explicarán
sus significados mis adelante.
I,a rapideznetadeflujo de masadel volumendecontrol se puede
calcular tomando en cuenta la masa transferida a través de la superficie de
control. Por ejemplo, la masa de A transferida a través del irea As Az,en x,
1
ser5 pA uAx A y Az x , o , en funci6n del vector de flujo. nA = p A uA , sería
A y Az /x.L a rapidez neta de flujo de masa de l a especie A , seri:
en la direccibn de x : nA,xAy A Z I ~ + -&n ,~, x Ay Azlx
en la direccih de y : n,,? A x A Z / , + &-~nA,?Ax Azjy
Y
en l a dirección de z : nA,i Ax Ayjz+Az- nA,zAx Ay[,
~a rapidez de acumulacibn de A en el volumen de control es:
JPA
AX Ay AZ
at
~
La ecuación diferencial dt?transferencia de masa 573
Si A se produce dentro del volumen de control por medio de una reacción

química a una rapidez r A , donde r A tiene las unidatdes: (masa producida de
A)/(volumen)(tiempo), la rapidez de producción de .4 es:
rA AX Ay AZ
Este término, correspondiente a la producción, es análogo al término que

corresponde a la generación de energía, que apareceen la ecuación diferencial
de la transferencia de energía, tal comose estudió en el capítulo 16.
Si se sustituye cada uno de los términos dela ecuación (25-l), se obtiene:
(25.2)
Si se evalúa esta expresión en el límite, cuando Ax,, Ay y AZ tienden a cero,

se obtiene:
(25-4)
La ecuación (25-4) es la ecuación de continuidad de la componente A . Como

n A , y ,y nA,z son las componentes rectangulalres del vector de flujo de
masa nA , la ecuación (25-4) se puede escribir de la siguiente forma:
V.n,+--rA=O
aPA (25-5)
at
Se puede obtener una ecuación semejante de continuidad para una se-
gunda especie, B , en la misma forma. Las ecuaciones diferenciales son:
a ns,x+ -
- a ng,y+ -
a nB,Z+ -
aPB - r13 = 0
(25-6)
dX ata y az
V-nB+--rB=O
aPB (25-7 )
at
574 Ecuaciones diferenciales de la transferencia demasa
donde rB es la rapidez con la quese produciráB dentrodel volumen d t control

por medio de una reacción química. Sumandb las ecuaciones (25-5) y (25-7),
se obtiene:
v * (nA + nB) + a ( P Aa+t P B ) - ( r A + r B ) = o (25-8)
Para una mezcla binaria de A y B, se tiene
debido a la ley de la conservación de la masa y sustituyendo estas relaciones

en la ecuación (25-8), se obtiene:
aP
v.pv+-=o (25-9)
at
que es la ecuación de continuidad correspondiente a la mezcla. La ecuación

(25-9) es idéntica a la ecuación de continuidad (9-2), correspondiente a un
fluido homogéneo.
La ecuación de continuidad de la mezcla y de unaespecie dada, se pueden
escribir en función de la derivada sustancial. La ecuación de continuidad de
la mezcla se puede reordenar y escribir de la siguiente manera, tal como se
demostró en el capítulo 9:
2 + p v .v = o (9-5)
Dt
A través de pasos matemáticos semejantes se puede obtener la ecuación de

continuidad dela especie A en función dela derivada sustancial.La ecuaciónes:
r, =o (25-10)
El mismo método se puede emplear en función de las unidades molares.

Si RA representa la rapidez de producción molar de B por unidad de volumen,
las ecuaciones molares equivalentes, son:
Componente A :
(25-11)
Formas especiales de la ecuación diferencial 575
Componente B:
(25-12)
Mezcla :
(25-13 )
Sin embargo sólo cuando la estequiometría de la reacción

es:
lo cual establece que se produce una molécula de B por cada mol de A que
desaparece y se puede decir que R, = -RB. En general, la ecuación de con-
tinuidad de la mezcla, en unidades molareses:
25.2 F O R M A S E S P E C I A L E S D E L A E C U A C I O N
DIFERENCIAL DE TRANSFERENCIA
DE MASA
Las formas especiales de la ecuación de contilnuidad válidas en las situa-

ciones comunes, son las que aparecen a continuación. Para poder usar las
ecuaciones en la evaluación de los perfiles de concentración, se reemplazan
los flujos, nA y N, por las ecuaciones apropiarlas obtenidas en el capítulo 24.
Estas expresiones son:
NA = -CDABVYA + Y A (NA+ r w (24-2 1)
o su equivalente
Y
(24-22)
o su equivalente
n A- P D A B V W A+ P A v
Al sustituir la ecuación (24-22) en la (25-5) se obtiene:
(25-15)
y al sustituir la ecuación (24-21) en la (25-11) se obtiene:
Cualquiera de las dos ecuaciones: (25-15) o (25- 16) se pueden usar en la

descripción de perfiles de concentración dentro de un sistema de difusión.
Ambas ecuaciones son completamente generales, sin embargo, son relativamen-
te complejas. Estas ecuaciones se pueden simplificar haciendo suposiciones
restrictivas. Las formas importantes de la ecuación de continuidad, con las
suposiciones adecuadas, incluyen:
i ) Si la densidad, p y el coeficiente de difusión, D A B ,se puede suponer

constante, la ecuación (25-15) se trhnsforma en:
Si se divide cada uno de los términos entre el peso molecular de A y se reor-

dena la expresión, se obtiene:
v VCA + 3CA
~
at
+ RA
= DABV2cA (25-17)
I ii)Si no hay término de producción, RA = O y si

coeficientesde
los
J
densidad y difusión se suponen constantes, la ecuación (25-17) se reduce a:
JCA
-
at
+v VCA =D + , ~ V ' C ~ (25-18)
Como (&,/at) + v VcA es la derivada sustancid de cA , cuando se reescribe el

lado izquierdo de la ecuación (25-18), se obtiene:
DCA
"
- D,,v'c* (25-19)
D
t
Formasespecialesde la ecuación diferencial 577
que es una ecuación análoga a la (16-14) de transferencia de calor,
(16-14)
DT
--=aV2T
Dt
donde 01 es la difusividad térmica. L a semejanza entre estas dos ecuaciones

sirve de base para las analogías señaladas entre las transferencias de calor y
de masa.
iii) En una situación en la cual no exista movimiento del fluido, v = O,

ni término de producción, R, = O , ni variación alguna en la difusividad o en
la densidad, la ecuación (25-18) se reduce a:
(25-20)
L a ecuación (25-20) se denomina segunda "ley" de Fick de la difusión. La

suposición de que no hay movimiento del fluido restringe su aplicabilidad a
la difusión en sólidos o líquidos estacionarios y a sistemas binarios de gases
o líquidos, donde N, es igual en magnitud pero actúa en dirección opuesta
a NB ; esto es, en el caso de la contradifusión equimolar. La ecuación (25-20)
es análoga a la segunda "ley" de Fourier de la conducción de calor:
aT
-= aV2T (16-18)
at
iv) Las ecuaciones (25-17), (25-18)y (25-20) se pueden simplificar aún

más cuando el proceso que se va a definires un proceso en estado permanente;
esto es: =O. La ecuacióncorrespondiente
una
a densidad constantey
una difusión también constante, es:
v VCA= DABV2cA+ RA (25-21)

Cuando la densidad es constante, lo mismo que la difusividad y no hay pro-
ducción química, RA = O , se obtiene:
v .VCA = DABV 2 C A (25-22)
Si además, v = O , la ecuación se reduce a
V2CA= o (25-23)
L a ecuación (25-23) es la ecuación de Laplace en función de la concentración
molar.
Todas las ecuaciones, desde la (25-15) hasta la (25-23) están escritas en

forma vectorial, por lo cual, todas son vilidas en cualquier sistema ortogonal
de coordenadas. La transformación de la ecuación al sistema deseado de coor-
denadas se realiza escribiendo el operador Laplaciano, V2, en la forma adecua-
da. La segunda ''ley'' de Fick de l a difusión, en coordenadas rectangulares es:
(25-24)
y en coordenadas cilíndricas
(25-25)
y en coordenadas esféricas es:
L a ecuación diferencial general de transferencia de masa de la componente A ,

o la ecuación de continuidad de A , en coordenadas rectangulares, es:
(25-27)
en coordenadas cilíndricas es:
(25-28)
y en coordenadas esféricas
25.3 C O N D I C I O N E S D E F R O N T E R A
ENCONTRADAS USUALMENTE
Se puede describir un proceso de transferencia de masa resolviendo una

de las ecuacionesdiferenciales de transferencia de masa,ya seausando las
condiciones iniciales de frontera, o ambas, que resulten convenientes para la
determinación de las constantes d e integración. Las condiciones iniciales y de
frontera utilizadas en el caso de la transferencia de masa son muy similares a
Condiciones de frontera encontradas usualmente 579
las usadas en la sección 16.3 para la transferencia de energía. Tal vez sea de
utilidad al estudiante repasar, en esa sección, lo referente a condiciones ini-
ciales y de frontera.
La condición inicial en los procesos de transferencia de masa es la con-
centración de la especie en difusión al principio del intervalo de tiempo bajo
estudio, expresada ya sea en unidades molares o en unidades de concentración
de masa. La concentración puede ser, sencillamente, una constante:
a t = O, cA = cA,,en unidades molares
a t = 0, PA = PA,, en unidades de masa
o puede ser una expresión más complicada, si la ,distribución de concentra-
ciones es una función de las variables al iniciar la medición del tiempo.
Las condicionesdefrontera que generalmente se encuentran, son las
siguientes:
1. Se puede especificar la concentración de u:na superficie. Esta concen-
tración puede estar en función de la concentración molar, cA = cA1de frac-
ciones molares, yA = yAl , en gases o x, = x A 1 ,en líquidos y sólidos; de
concentración de masa, p A = p A l o de fracción de masa, wA = w A 1 . Cuando
el sistema es un gas, la concentración se relaciona con la presión parcial por
medio de la ley de Dalton;asimismo, la concentración puede estaren función
de la presión parcial, pA =p = y A l P. En el caso específico de la difusión de
un líquido a una fase gaseosa, la condición de frontera en la superficie del lí-
quido está definida para una solución líquida ideal, por medio de la ley de
Raoult, de la siguiente manerap, = x A PA , donde x, es la fracción molar en
el líquido y PA es la presión del vapor de la especie A evaluada a la tempera-
tura del líquido.
2. El flujo de masa en la superficie puede estar especificado; por ejemplo:
.
j, = j, o NA = NA1 Los casos de interés en inge:niería incluyen aquel en el
que el flujo es cero debido a una superficie impermeableaquellos
y exi los que
el flujo de masa está especificado en la forma siguiente:
3 . Larapidezdeunareacciónquímicapuedeestar especificada, por

ejemplo, si A desaparece en la frontera debido a una reacción química de
primer orden, se puede escribir N A I = k 1 C ~ 1 ,donde kl es la rapidez constante
en una reacción de primer orden. Cuando laespecie en difusión desaparece en
la frontera debidoa una reacción instantánea, la concentración de esa especie
se supone, a menudo,igual a cero.
4. Cuando un fluido fluye en la fase para la cual está escrita la ecuación
diferencial de masa, la especie puede perderse de la fase de interés por medio de
la transferencia convectiva de masa. El flujo de masa en la frontera, se define,
entonces, así:
donde CA, es la concentración de la corriente de fluido adyacente a l a super-

ficie y kc esel coeficiente detransferenciade masa deconvección que se
definió en la sección 24.3.
Se aplicarán estas condiciones iniciales y de frontera en la solución de
problemas de difusión molecular en los capítulos 26 y 27. En seguida se pre-
sentará un ejemplo de reducción en la ecuación diferencial general de trans-
ferencia de masa de manera que solamente quedenlos términos relevantes.
EJEMPLO 1
En una barra cilíndrica de combustible nuclear que contiene material fisionable, la

rapidez de producción de neutrones es proporcional a la producción de los mismos. Use
una de las ecuaciones diferenciales generales de transferencia de masa para escribir la ecua-
ción diferencial que describa el proceso de transferencia de masa. Haga una listade las
condiciones obvias de frontera.
En el caso de la componente A , la ecuación (25-1 1) estable lo siguiente:
ac,
V.N,+--R,=O
at
Como la rapidez de producción es proporcional a la concentración neutrónica R , = kc,.

En los casos de difusión en sólidos, en los que la contribución del movimiento global es
cero,
N, = -DABVc,
Al sustituir estas relaciones en la ecuación ( 2 5 - 1 1 ) , se obtiene:
at
Si el cilindro es relativamente largo comparado con su radio J 2 ~ A / = O y2 si la concentra-

J~
ción no varía con el ángulo, 8, ¿?cA/d8z = O ; la ecuación se reduce a:
La única condición de frontera que resulta obvia es:
la cual requiere que la concentración de la especie en difusión sea finita en el centro del
cilindro.
Problemas 581
L a ecuación diferencial general de transferencia de masa se ha obtenido

en este capítulo para la componente en difusión de una mezcla. Se presentaron
las formas de la ecuación general, válidas en situaciones específicas. También
se listaron las condiciones de frontera encontradas con mayor frecuenciaen los
procesos de transferencia de masa.
PROBLEMAS
25.1 Obtengalaecuación (25-11) correspondientealacomponente A en

función de lasunidades molares. empezartdo con la expresión de vo-
lumen de control correspondiente a la conservación de la masa.
25.2 Obtenga la ecuación de continuidad de la. componente A , utilizando
un volumen de control cilíndrico e infinitamentegrande y coordenadas
cilíndricas.
25.3 Transforme la ecuación (25-27) de coordenadas rectangulares:
a una ecuación equivalenteen coordenadas cilíndricas.

"$5.4 Una celda de difusión de Arnold es un apara.to sencillo usado para medir
coeficientes de difusión de gases. Un estanque de líquido se mantiene
en el fondo de un tubo de diámetro pequeño. Un gas, insoluble en el
líquido, fluye a través de la boda del tubo extrayendo los vapores de A
que se difunden a travésdel espacio de gas que se encuentra sobre el
estanque de líquido. En condiciones isotérmicas e isobáricas, la evapo-
ración de A es un proceso en estado permanente.
Reduzca la ecuación diferencial general de transferencia de masa para
escribir la ecuación diferencial específica que describa este proceso de
transferencia de masa.
iCuál sería la forma de la ley de Fick, correspondiente a la especie A
que se substituiría en la ecuación diferenc:ial? Proporcione dos condi-
ciones de frontera que podrían ser utilizadas para resolver la ecuación
diferencial resultante.
-Flow of gas B
Líquido A puro
28
TRANSFERENCIA CONVECTIVA
DE MASA
La transferencia de masa por convección es la que se realiza entre una

superficie límite y un fluido en movimiento o entre dos fluidos en movimiento,
relativamente no miscibles. La ecuación de rapidez correspondiente a la trans-
ferencia de masa ya ha sido expresada anteriormente en la forma:
N* = kc AcA 2) (24-5
en la que el flujo de masa, NA , se lleva a cabo en la dirección en la que dismi-
nuye la concentración. Esta sencilla ecuación es la relación que define a kc,
que es el coeficiente de transferencia conuectiua de masa. Es análogo ala
ecuación que define el coeficiente de transferencia convectiva de calor:
q/A=h AT (15-11)
Si se hace memoria, se recordará que el estudio para encontrar el coe-
ficiente de transferencia de calor que se hizo en el capítulo 15, no fue una
empresa fácil.Amboscoeficientes de transporte están relacionados con las
propiedades del fluido, con sus características dinkmicas y con la geometría
del sistema específico bajo estudio.
A causa de la estrecha semejanza entre las ecuaciones de transferencia
convectiva de calor y de masa, es de esperarse qwe el coeficiente de transfe-
rencia de calor del capítulo 19 se pueda aplicar al coeficiente de transferen-
cia de masa y esto es, precisamente, lo que se hará en los análisis siguientes.
Se hará un uso considerable de las secuelas y conceptos de los capítulos 9 a14.
28.1 C O N S I D E R A C I O N E S F U N D A M E N T A L E S A C E R C A
DE LA TRANSFERENCIA CONVECTIVA DE MASA
Se recordará el estudio anterior referente a un fluido escurriendo drede-

dor de una superficie, que existe una capa, a veces delgada, cerca de la super-
657
658 Transferenciaconvectiva de masa
ficie en la que el flujo es laminar. Por lo tanto, la transferencia molecular de

masa siempre estará presente y siempre iendrá un papel importante en cual-
quier proceso de convección. Si el flujo de fluido es laminar, entonces toda
la transferencia entre la superficie y el fluido en movimientose llevará a cabo
por medios moleculares. Si por otra parte, el flujo de fluido es turbulento,
habrá un movimiento físico de paquetes de materia alo largo de las líneas de
flujo transportadas por los remolinospresentesen el flujo turbulento. Así
como en el caso de la transferencia de calor, a unagran rapidez de transferen-
cia de calor se asocian condiciones turbulentas. La distinción entre los flujos
laminar y turbulento será de vital importancia en cualquier situación convec-
tiva.
La capa límite hidrodinámica que se analizó en el capítulo 12, tiene un
papel muy importante en la transferencia convectiva de masa. También se
definirá y analizará una capa límite de concentración que resultará vital para
el análisis del proceso de transferencia convectiva de masa. Esta capaes similar,
pero no necesariamente del mismo espesor que la capa térmica límite que se
estudió en el capítulo 19.
Cuando la transferencia de masa incluye a un soluto que se disuelve con
rapidez constante desde una superficie sólida y después se difunde a un fluido
en movimiento, el coeficiente de transferenciaconvectiva de masa se define así:
En esta ecuación, el flujo NA representa los moles de soluto A que abandonan

la intercara por unidad de tiempo y por unidad de área interfacial. La com-
posición del soluto en el fluido, en la intercara, c ~ ,es~la, composición del
fluido en equilibrio con el sólido a la temperatura y presión del sistema. La
cantidad cA , representa la composiciónen algún punto dentro dela fase fluida.
Cuando se define la capa límite de la concentración, se puede escoger a cA
como la concentración de la componente A en el borde de la capa límite y
sele puede llamar cA-. Si el flujo tuviera lugar en un conducto cerrado, la
composición, c A , podría serla concentración global o concentración de
nzezcla homogénea. La composición de mezcla homogénea es la concentra-
ción que se mediría si se recogiera y se mezclara perfectamente el flujo en
un plano, esto es,se tendría una composición promedio del flujo global.
Existen cuatro métodos de evaluación de los coeficientes de transferencia
convectiva de masa que se estudiarán en este capítulo, que son:
1. análisis dimensional acoplado con experimentos

2. análisis exacto .de capa límite
3. análisis aproximado de capa límite
4. analogía entre la transferencia de momento, energía
y masa.
En las secciones siguientes se estudiará cada uno de estos métodos.

Parámetros importantes en la transferencia convectiva de masa 659
28.2 P A R A M E T R O S I M P O R T A N T E S E N L A
TRANSFERENCIA CONVECTIVA
DE MASA
A menudo se utilizan parámetros adimensionales para relacionarlos da-

tos relativos a la transferencia convectiva. En la transferencia de momento,
nos hemos topadoya con los números de Reynolds, y de Euler. En la relación
de datos de transferenciaconvectiva de calor, son importantes los números de
Prandtl y Nusselt. Algunos de los mismos parámetrosjunto con ciertas razones
adimensionales recién definidas, resultarán de utilidad en la correlación de
los datos correspondientes a la transferencia de masa. En esta sección se es-
tudiará la interpretación físicad e tres de estas razo:nes.
Las difusividades moleculares de los tres fenilmenos de transferencia se
han definido así:
difusividad de momento = u =p/p
k
difusividad térmica =cy = -
PC"
difusividad de la masa = D A R
Tal como se hizo notar anteriormente, dada una de las difusividades tiene las
dimensiones L 2 / t , por lo que la razón de cualesquiera dos de ellas debe ser
adimensional. La razón de la difusividad molecular de momento a la difusi-
vidad molecular de la masa se denomina número de Schmidt:
difusividad de momento
-Sc="=-
- v El. (28-2)
dqusividad de masa DAB PDAB
El número de Schmidt tiene una importancia en la transferencia convectiva

de masa, análoga a la del número de Prandtl, enla transferencia convectiva de
calor. La razón de la difusividad molecular del calor a la de la masa, se deno-
mina número de Lewis:
difusividad térmica = = k (28-3)

difusividad de masa PCpDAB
El número de Lewis aparece cuando un proceso con,sta de la transferencia con-

vectiva simultánea de masa y energía. Los número,s de Schmidt y Lewis s o n
combinaciones de las propiedades de los fluidos, por lo tanto, cada uno de
ellos se puede considerar como una propiedad del sistema endifusión.
660 Transferencia convectiva de masa
Pensemos ahora en la transferencia demasa del soluto A desde un sólido

hasta un fluido que pasa alrededor de la superficie del sólido. En la figura
28.1 aparece el perfil de la concentración.
Figura 28.1 Perfiles de la velocidad y la Concentración correspondientes a un

fluido que circula alrededor de una superficie sólida.
Para este caso, la transferencia demasa entre la superficie y el fluido se puede

escribir en la forma:
Como la transferencia de masa en la superficie se realiza por difusión mo-

lecular, la transferencia de masa también se puede escribir por medio de:
cuando la concentración de frontera,CA,,,es constante, esta ecuación se sim-

plifica, reduciéndose a:
(28-5)
Las ecuaciónes (28-4) y (28-5) se pueden igualar, ya que definen el mismo

flujo de la componente A que abandona la superficie y entra al fluido. Esto
produce la relación:
(28-6)
Si se multiplican ambos lados de la ecuación (28-6) por una longitud signi-

ficativa, L , se obtiene la siguiente expresión adimensional:
(28-7)
Análisis de la transferencia convectiva de masa 661
El lado derecho de la ecuación (28-7) es la relación entre el gradiente de la

concentración en la superficie y un gradiente globad o de referencia de la con-
centración. Por lo tanto, se puede considerar como la relación entre la resis-
tencia a la transferencia molecular de masa y la resistencia del fluido a la
transferencia convectiva de masa. Como la secuela de la ecuación (28-7) es
semejante a la de la ecuazión (19-5) del número cle Nusselt encontrado en la
transierencia convectiva de calor, la razónk,L/Dp,B se denomina número de
Nusselt de transfrencia de masa, NuAB.También se le conoce como número
de Sherwood, Sh.
Estos tres pargmetros Sc, Nu,, o Sh y Le, aparecerán en los análisis de
transferencia convectiva de masa enlas secciones siguientes.
28.3 A N A L I S I S D I M E N S I O N A L D E L A T R A N S i F E R E N C I A
CONVECTIVA DE MASA
El a d i s i s dimensional predice los diversos parámetros adimensionales

que resultan útiles en la correlación de los datos experimentales. Hay dos
procesos importantes de transferencia de masa que se estudiarán: la trans-
ferencia de masa hacia una corriente que fluye en condiciones de convección
forzada y de transferencia demasa hacia unafase (quese mueve en condiciones
de convección natural.
TRANSFERENCIA HACIA UNA CORRIENTE QUE FLUYE

EN CONDICIONES DE CONVECCION FORZADA
Estudiemos,ahora, latransferenciademasadesde las paredes de un

conducto circular hacia un fluido que circula a través del conducto. La trans-
ferencia es el resultado de la fuerza impulsora de la concentración, c ~ -cA.
, , ~
Las variables importantes, sus símbolos y sus represent?ciones dimensionales,
son los que aparecen a continuación:
Vmiab le Símbolo Dimensiones
diámetro del tubo D

uido del densidad P
fluido del viscosidad fi
fluido del velocidad 0
fluido del difusividad DAB
coeficiente de transferencia de masa kc
Las variables anteriores incluyen términos descriptivos de la geometría del

sistema, las propiedades del flujo y del fluido y la cantidad de principal in-
terés, k,.
Por mediodel método de Buckinghamde agrupamientb de variables,

que se estudió en el capítulo 11, se puede determinar que habrá tres grupos
adimensionales. Los tres grupos pi que se van a formar tomando a D A B , p, y
D como las variables principales, son:
Si se escribe 7~~ en forma dimensional,
y se igualan los exponentes de las dimensiones fundamentales de ambos lados

de la ecuación,se tiene para:
L: 0 = 2 ~ - 3 b + c + l
t: o=-u-1
Y
M : O=b
La solución a estas ecuaciones, correspondiente a los tres exponentes incóg-
nitos, produce el resultado:
U=-]
b=O
Y
c=l
así pues: 7 r 1 = k,D/DAB,que es el número de Nusselt de transferencia de masa
o bien, el número de Sherwood.Los otros dos grupospi se pueden determinar
en la misma forma, dando:
Dv
nz=-
DAB
Y
P -- sc
r3=--
@AB
que es el número de Schmidt. Si se divide IT^ entre r3 se obtiene:
EL
Análisis de la transferencia convectiva de masa 663
el número de Reynolds. El resultado del análisis dimensional de la transfe-

rencia de masa en condiciones de convección forzada en un conducto circular
indica que la correlación podría ser de la forma:
NuAB= f(Re, Sc) (28-8)
que es análoga a la correlación de transferencia de calor:
Nu =f(Re, Pr) (19-7)
TRANSFERENCIA A UNA FASE CUYO MOVIMIENTO

. SE DEBE A LA CONVECCION NATURAL
Las corrientes de convección natural se crearán si existe cualquier va-

riación en la densidad en una fase líquida o en una gaseosa. La variación en
la densidad se puede deber a diferencias de temperatura o a diferencias de
concentración, relativamente grandes.
En el caso en el que la convección natural incluye la transferencia demasa
desde una pared vertical plana hasta un fluido adyacente, las variables diferirán
de las utilizadas en el análisis de la convección forzada. Las variables impor-
tantes, sus símbolos y sus representaciones dimensionales, aparecenen la lista
siguiente:
Variable Sz'mbo,!o Dimensiones
longitud característica L L
difusividad del fluido DAB Lyt
densidad del fluido P MIL3
viscosidad del fluido P MILT
fuerza boyante g 'PA. MIL?'
coeficiente de transferencia de masa kc Lit.
Debido al teorema de Buckingham, puede predecirse que habrá tres grupos

adimensionales, los cuales serán, con DA.B,L, y p como variables principales:
Si se resuelven los tres grupos pi, se obtiene:

el número da Nusselt de transferencia de masa, o número de Sherwood:
el recíproco del número de Schmidt y
El resultado del análisis dimensional de la transferencia de masa por convec-

ción natural sugiere una interrelación de la forma:
SC)
NUAB= ~ ( G ~ A B, (28-9)
Las relaciones que se obtuvieron mediante el análisis dimensional, tanto
para la convección natural como para la forzada, sugieren que una correlación
de los datos experimentales puede estaren función de tresvariables, en vez de
las seis originales. Esta reducción en el número de variables ha ayudado a los
investigadores que han sugerido interrelaciones de esta forma para producir
muchas de las ecuaciones empíricas que aparecen en el capítulo 30.
28.4 A N A L I S I S E X A C T O D E L A C O N C E N T R A C I O N
LAMINAR DE LA CAPA LIMITE
Blasius obtuvo una solución exacta, correspondiente a la capa límite

hidrodinámica para un flujo paralelo a una superficie plana. Esta solución se
examinó en la sección 12.5. En la sección 19.4 se hizo una extensión a la
solución de Blasius para explicar la transferencia convectiva de calor. En una
forma absolutamente análoga se ampliará, ahora, la solución de Blasius, para
incluir la transferencia convectiva de masa para lamisma geometría yel mismo
flujo laminar.
Las ecuaciones correspondientes a la capa límite, estudiadas en la trans-
ferencia de momento en estado permanente, incluyeron la ecuación bidimen-
sional e incompresible, de continuidad:
au, au
-+".?=o (12-1 l b )
ax ay
Análisis de l a concentración laminar de la capa límite 665
y la ecuación de movimiento en la dirección de x, para una presión y una v

constantes,
(12-1la)
La ecuación que describe la transferencia de energía de un flujo isobárico en

estado permanente, incompresible y bidimensionatl, con una difusividad tér-
mica constante, correspondiente a la capa térmica límite, es:
(19-15)
Hay una ecuación que se aplica a la transferencia de masa dentro de una

capa límite de concentración cuando no hay producción de la componente
en difusión y cuando la segunda derivada de cA con respecto a x, d2C,4/aX2,
es mucho más pequeña en magnitud que la segunda derivada de cA con res-
pecto a y . Esta ecuación, escrita para un flujo permanente incompresible y
bidimensional, cuya difusividad de masa es constante, es:
(28-10)
La capa límite de concentración aparece en forma esquemática, en la figura

28.2. Las siguientes son las condiciones de frontera correspondientes a las
tres capas límite:
Orilla de
limite
la capa n
Figura 28.2 Capa límite de concentración de un flujo laminar alrededor de una

superficie plana.
o, como la velocidad en la dirección de x en la pared, vXTses igual a cero
T - T, T - T,
térmico -- -0 eny=O Y -- -1 eny=m
Tm - T s T m- T s
Yconcentración
La semejanza entre las tres ecuaciones diferenciales: (12-1la), (19-15) y

(28-10) así como las condiciones de frontera, sugieren que deben obtenerse
soluciones semejantes para los tres fenómenos de transferencia. Enel capítulo
19 se modificó la solución de Blasius correspondiente a la ecuación (12-1l a )
y se le aplicó con éxito en la explicación de la transferencia convectiva de
calor cuando la razón del momento a la difusividad térmica, v / a = Pr = 1 se
utiliza el mismo tipo de solución paradescribir la transferencia convectiva de
masa cuando la razón: momento-difusividad de masa, v / D = Sc = 1. Usando
la nomenclatura definida en el capítulo 12:
(28-12)
al aplicar el resultado de Blasius ala capa límite de

la concentración se obtiene:
La ecuación (28-13) se puede reordenar para obtener una expresión corres-

pondiente al gradiente de la concentración en la superficie,
(28-14)
Es importante recordar que la solución de Blasius correspondiente a la ecua-

ción (12-lla) no incluía una velocidad en la superficie en la direccibn de x ,
por lo tanto, la ecuación (28-14) hacc la suposición, muy importante, de que
La rapidez con la cual entra o sale la masa de la capa límite enla superficie, es
tan pequeíia que no altera el perfil de la velocidad predicho por la solución
de Blasius.
Cuando la velocidad en la dirección de y en la superficie, UW es igual a
cero, el término correspondiente a la contribución global, en la ecuación de
Fick, el flujo de masa en la dirección de y también es igual a cero. La trans-
ferencia de masa de la superficie pIana de la capa laminarlímite, está descrita
por la ecuación:
(28-15)
Análisis de la concentración de la capa límite 667
AI sustituir la ecuación (28-14) en la (28-15),se okltiene:
(28-16)
El flujo de masa de la componente en difusión se definió en función del coe-

ficiente de transferencia de masa por medio de:
Los lados derechos de las ecuaciones (28-16) y (28-4) se pueden igualar, ob-
teniéndose:
DAB
kc =“-[0.332
X
ó
-= NuAB = 0.332 Re,1’2 (28-17)
DAB
La ecuación (28-17) está restringida a sistemas cuycl número de Schmidt, Sc,

sea igual a uno y con una rapidez de transferencia de masa entre la placa plana
y la capa límite, muy baja,
En la figura 28.3 puede observarse una representación gráfica de la so-
lución a la ecuación de capa límite de concentración, ecuación (28-lo), ob-
tenida por Heartnett y Eckert.” En esta figura se muestran las curvas que
representan los valores positivos y negativos del parámetro límite de la super-
ficie, (uY,JudRex)”?Los valores positivos son válidos cuando la transferencia
de masa de la placa plana ocurre hacia la capa límite y los valores negativos
describen la transferencia de masa del fluido la placa.
a Cuando el parámetro de
la superficie de frontera tiende a cero, la rapidez de transferencia de masa
disminuye hasta que se considera que no tiene ningún efecto en el perfil de
velocidad. La pendiente de la recta cero, calculada eny = O es 0.332, tal como
lo predice la ecuación (28-13).
En la mayoría de las operaciones físicas que se relacionan con la trans-
ferencia de masa, el parámetro de la superficie límite es despreciable y se uti-
liza el tipo de solución de Blasius, de baja transferencia de masa, para definir
la transferencia hacia la capa laminar límite. La vaporización de un material
volátil hacia una corriente gaseosa que fluye a baja presión, es un caso en el
que no se puede hacer la suposición de baja transferencia de masa.
*J. P. Hartnett y E. G . Eckert, Trans. A .

R. S. M. E., 13, 247 (1957).
”. .. .. .” . c_x._ . . . .. . * ” ” .
O 2 4 6 8 10 12
11 = ;(Re,) lh
Figura 28.3 Perfiles de la concentración correspondientes a la transferencia de

masa en una capa laminar límite sobre una placa plana.
Para un fluio cuyo número de Schmidt sea diferente de la unidad, se

pueden definir curvas similares a las de la figura 28.3. La semejanza de las
ecuaciones diferenciales y de las condiciones de frontera sugiere el empleo
de un método análogoalasolución de Pohlhausen,correspondiente a la
transferencia de calor convectivo, para la transferencia convectiva de masa.
La capa límite de concentración está relacionada con la capa hidrodinámica
límite por medio dela ecuación:
(28-18)
donde 6 es el espesor de la capa hidrodinámica límite y 6, es el grosor de la

capa límite de concentración, por lo cual, el término de Blasius, q, se debe
multiplicar por S C ” ~ .En la figura 28.14 aparece una gráfica de concentración
adimensional contraqSc”’ correspondiente a uy.s =O. La variación en la con-
centración, expresada en esta forma conduce a una expresión correspondiente
Análisis de la concentración laminar de la capa límite 669
al coeficiente de transferencia convectiva de masa, semejante a la ecuación

(28-17). En y = O, el gradiente de la concentración es:
(28-19)
el cual, al emplearse con la ecuación (28-15)da origen a la ecuación:
kcx - Nu,,AB = 0.332

"
SC"~ (28-20)
DAB
El coeficiente medio de transferencia de masa, que se aplica sobre una
placa de ancho W y longitud L , se puede obtener por integración. En una placa
de estas dimensiones, se puede evaluar la rapidez, total de transferencia de
masa, W , , por medio de:
WA = k,A(cA,,-cA,d= JA
k,(cA,s-CA,m) d~
= k, WL (cA,, - C A . ~
- 0.332DABRe,"2Sc'/3 dA
- (C.4.S - cA,m X
Así pues:
Y
k,L
"
- NuL,,AB = 0.664Re,_1/2Sc1'3 (28-21)
DAB
El número local deNusselt a una distanciax , en. elsentido dela corriente,
se relaciona con el número medio de Nusselt correspondiente a la placa, por
medio de la expresión:
NUL,AB =2NUx,AB(x=L (28-22)

que es análoga al resultado de la transferencia convectiva de calor:
NUL = 2Nu,lx,L (19-27)
Las ecuaciones (28-20) y (28-21) se han verificado experimentalmente. Es
importante hacer notar que este análisis, totalmente diferente, ha producido
*W. J. Christian y S. P. Kezíos, A.1.Ch.E. J., 5 , 61 (1959).
resultados de la misma forma predicha en la sección 28.3 mediante el análisis

dimensional, para la transferencia forzada-convectiva de masa,
NuAB = f(Re, Sc) (28-8)

Si se recuerdan los perfiles adimensionales de concentración de Harnett
y Eckert, tal como aparecen en la figura 28.3, se podrá observar que cuando
se calcula la pendiente de cada una de las curvas en y = O, ésta disminuye
cuando el parámetro positivo de la superficie límite ( v /v,)(Re)”*, aumenta.
y:
Como la magnitud del coeficiente de transferencia esta directamente relacio-
con
nado ecuación:
lalapendiente
medio
depor .
(28-23)
la disminución de Ia pendiente indica que los sistemas cuyos valores del pará-
metro de la superficie límite son altos, tienen coeficientes más bajos de trans-
ferencia de masa.
Cuando la energía y la masa se transfieren a través de la capa laminar
límite, los perfiles adimensionales de la figura 28.3 también pueden representar
los perfiles adimensionales de temperatura cuando los números de Prandtl y
Schmidt son iguales a la unidad. En el párrafo anterior, se señaló que el coe-
ficiente de transferencia de masa disminuye en magnitud cuando se transfiere
masa de la superficie a la capa límite, por lo tanto, es de esperarse que el coe-
ficientedetransferenciadecalordisminuya al transferirse masa a la capa
límite. Esto puede realizarse haciendo pasar un fluido a través de una placa
porosa hacia afuera, llegando a la capa límite, o sublimando el material mismo
de la placa. Estos procesos simultáneos de transferencia de calor y masa, a los
cuales se llama a menudo enfriamiento por transpiración y ablación, respec-
tivamente, se utilizanparaayudarareducir los enormesefectosdelcalor
durante el regreso de un proyectil ala atmósfera terrestre.
EJEMPLO 1
~1 coeficiente de transferencia de masa correspondiente a una capa límite turbulenta

formada sobre una placa plana se ha correlacionado, en función del número local de N U d t , -
por medio de:
Nu,,,, = 0.0292Re~’5Sc’’’ (28-24)
donde X es la distancia en el sentido de la corriente, al borde de ataque de la placa Plana.

transición de flujo laminar a turbulento ocurre en Re, = 3 X 10’.
a ) Encuéntrese una expresión que corresponda al coeficiente medio de transferencia
de masa en una placa plana de longitud L , .
Por definición:
(28-25)
Análisis de la concentración’de la capa límite 671
donde L, es la distancia medida del borde de ataqueal punto de transición. kc,,,, está defi- .
nida por la ecuación (28-20):
DAB
k,,,, = 0.332---(Re,)’/2(Sc)’/3
X
kc,turbestá
definida
por ecuación
la (28-24): +
Das
kc.turt,= 0.0292-(Re,)“”(S~)‘/~
X
AI sustituir estas dos ecuaciones en la expresión correspondiente al coeficiente medio de

transferencia de masa, se obtiene:
donde L, es la distancia del borde de ataqueal punto de traasición, enel cualRe, = 3 x 10’.
?
1 I
a, i 4/5
0.664DAB(E)
1 /2
+ 0.0361DA,(
(SC)’/~J!,:/~ ): - &)“/’I
(SC)”~[(L)~”
- 0.664D,,(Re,)”2(Sc)”+0.0361D,,(Sc)”3[(Re,~)4’5 -(Re,)“/’]
kc = (28-26)
L
b ) Un vaso de alcohol etílico se vokó accidentalmentesobre un banco de laboratorio,
cubriendo la superficie lisa superior del mismo. El extractor de aire quese encuentra en el
. laboratorio produjo un flujo de aire de 6 m/seg, paralelo a h superficie que fluye a través
del banco de 1 m de ancho. El aire se mantuvo a 289 K y 1 atm. (1.013 x lo5 Pa). La
presión del vapor de alcohol etílico a 289 K es de 4,000 Pa. Determínese la cantidad de
alcohol que se evapora de un área de un metro cuadrado cadar 60 seg.
x
A 289 K, la viscosidad cinemática es de 1.48 10-5m2/seg. y la difusividad de la
masa deetanol en aire es de 1.26 x m2/seg. En este sistema:
sc==
- Y 1.48 X lo-’ m2/s
= 1.17
DAB 1.26 X lo-’ mz/s
El número de Reynolds correspondiente a toda la longitud de: 1 m. es:
Como esta cantidad es mayor que 3 x l o 5 , se observa que habrá un punto de transición
donde la capa límite cambia de flujo laminar a turbulento. El :punto de transición se puede
evaluar a partir del número de transición de Reynolds, Ret = 3 x 10’.
L Re,v (3 x 10’)(1.48 x
=-=
m2/s)
= 0.74 m
U 6 mfs
Se puede calcular el coeficiente medio de transferencia de masa usando la ecuación (28-26)

obtenida en la parte (a)
- 0.664(1.26x lo-’ mZ/s)(3x 105)1’2(1.17)’’3

kc =
lm
+ 0.0361(1.26~lO”m’/~)(l.17)~’~[(4.05~
l m
105)4/5-(3x
= 0.00483 + 0.00313 = 0.00796 m/s

Se puede calcular la concentración del alcohol etílico en el vapor inmediatamente encima
de la superficie del líquido, por medio de la relación:
PA 4000 P a
C& = -=
P a . .m3
RT (8.314-)(289 K mol K)
= 1.66 mo1/m3
La cantidad de alcohol que abandona la superficie es de k; ( c ~- ,c ~~ 6:

, ~ )
WA = (0.00796 m/s)(l m”(1.66 mol/m3)

= 1.32 x lo-’ mol/s = 0.793 mo1/60 S (0.104 lb mole/hr)
-
28.5 A N A L I S I S A P R O X I M A D O D E L A C A P A
LIMITE DE CONCENTRACION
Cuando el flujo es de un tipo diferente al laminar o cuando la configu-

ración es otra que no sea una placa plana, existen actualmente pocas soluciones
exactas correspondientes a la transferencia en una capa limitada. El método
aproximado, encontrado por von Karmán para describir la capa límite hidro-
dinámica, se puede utilizar en el an6lisis de la capa límite de concentración.
En los capítulos 1 2 y 19 se estudió el uso de este método.
Analicese, ahora, el volumen de control localizado en l a capa límite de
concentración, tal como aparece enla figura 28.5. Este volumen representado
por las líneas punteadas, tiene un anchoigual a Ax, una altura igual al espesor
!V
Figura 28.5 Volumen de control de la capa límite de concentración.

Análisis aproximado de la capa límite de concentración 673
de la capa límite de concentración, S, y profundidad unitaria. Si se hace un

balance de la masa molar en estado permanente sobre el volumen de control,
se tendrá la relación:
wA, -k WA, + WA, = wA, (28-27)
donde W, es la rapidez molar de transferencia de masa de la componente A .
En cada una de las superficies, la rapidez se expresa así:
Y
WA,= k c (CAS - C A . ~ A) X
En términos de rapidez molar, se puede reescribir la ecuación (28-27) en la
forma:
(28-28)
AI redordenar y dividir cada uno de los términos entre Ax y calcular el resul-
tado en el límite, cuando Ax tiende a cero, se obtiene:
ó
d
& j,,,(CA - cA,m) ux dy = k c ( C A S - CA,m) (28-29)
La ecuación (28-29) es análoga a las ecuaciones (12-37) y (19-30). Para

poder resolver la ecuación (28-29) deben conocerse los perfiles de la velocidad
y la concentración que usualmente se desconocen y se deben suponer. Las
siguientes son algunasde las condiciones de frontera que deben satisfacer
las condiciones supuestas de frontera, son:
1) u, = o eny = O
2) u, = urn eny = 6
avx -
3) --O eny = 6
JY
674 Transferencia convectiva demasa
Y
4) a2ux
2-0 eny = O
El perfil de concentración supuesto debe satisfacerlas condiciones correspon-

dientes de frontera en función delas concentraciones.
3) en y = S, (28-32)
4) eny=O (28-33)
Si se recuerda el flujo laminar paralelo a una superficie plana, se puede

usar la ecuación integral de von Kármán, que es la (28-29), para obtener una
soluciónaproximada. Los resultados se puedencompararconlasolución
exacta, ecuación (28-20) y verificar la exactitud de las suposiciones hechas
acerca de la velocidad y de los perfiles de concentración. Como primera aproxi-
mación, vamos a tomar una expresión en serie de potencias, para la variación
de la concentración con y.
2 3
~A-~A,~=a+by+cy+dy
Al aplicar las condiciones de fronterase obtendrá la expresión:
(28-34)
Si se supone que el perfil de la velocidad tiene la misma formaen serie de po-

tencias, entonces la expresión resultante, tal como se vio en el capítulo 1 2 es:
( 12-40)
Cuando se sustituyen las ecuaciones (28-34)y (12-40) en la expresión integral

(28-29) y se resuelve, se obtiene:
Nu,,AB = 0.36Re,”2S~”’ (28-35)
resultadoque se aproximabastantealasoluciónexactaexpresadaenla
ecuación (28-17).
Aunque esta no es la expresión correcta,es suficientemente aproximada
para indicar que se puede usar el método integral con bastante confianza en
Analogía de transferencia de la masa, energía y momento 675
otras situaciones en las que no se conozca una solución exacta. La exactitud

del método depende totalmente de la habilidad para suponer buenos perfiles
de velocidad y concentración.
La ecuación integral de von Kármán, ecuaciijn (28-29), también se ha
utilizado para obtener una solución aproximada de la capa límite turbulenta
sobre una placa plana. Si se toma el valor aproximado del perfil de la velo-
cidad igual a:
v,=(Y+py1"
y el valor aproximado de l a concentración, igual a:

1/ I
CA-CA,~=V+~Y
el número local de Nusselt de la capa turbulenta es:

= 0.0292 (28-24)
Se anima al lector a llevar a cabo la obtención de la ecuación (28-24) que se
solicita en uno de los problemas que aparecen al final del capítulo.
28.6 A N A L O G I A S D E T R A N S F E R E N C I A D E
MASA, ENERGIA Y MOMENTO
En el análisis que se hizo anteriormente de la transferencia convectiva

de masa, se,reconocieron las semejanzas que existen entre las ecuaciones di-
ferenciales de transferencia de momento, energía y rnasa y en las condiciones
de frontera, cuandose expresaron los gradientes de tr,ansferencia en función de
variables adimensionales. Estas semejanzas nos han permitido predecir solu-
ciones correspondientes a los procesos semejantes de transferencia. En esta
sección se estudiarán diversas analogías de entre lots fenómenos de transfe-
rencia, que se han propuesto debido a la semejanza de sus mecanismos. Las
analogías resultan de utilidad para comprender los fenómenos de transferencia
y también como un medio satisfactorio de predecir el comportamiento de
los sistemas para los cuales existen datos cuantitativos limitados disponibles.
La semejanza entre los fenómenos de transferencia y, por lo tanto, la
existencia de las analogías, requieren de la existencia de las siguientes cinco
condiciones dentro del sistema:
1. Que las propiedades físicas sean constantes.
2. Que no haya producción de masa o energía dentro del sistema. Esto,
desde luego, implica que no pueden ocurrir reaccione:: químicas homogéneas.
3. Que no exista emisión o absorción de energía radiante.
4. Que no haya disipación viscosa.
5. Que el perfil de la velocidad no esté afectado por la transferencia de
masa, por 10 cual habrá una transferencia lenta de masa.
ANALOGIA DE REYNOLDS
Reynolds* fue quien primero se dio cuenta del comportamiento análogo

de la transferencia de energía y momento. Aunque esta analogía es de zpli-
cación limitada, ha servido como catalizador para encontrar mejores analogías
y ha sido utilizada con éxito en elanálisisde los complejos fenómenos de
capa límite de la aerodinámica.
Reynolds postuló que los mecanismos de transferencia y energía y mo-
mento eran idénticos.En el estudiorealizadoanteriomente,acerca de las
capas laminares límite, se ha visto que esto se cumple cuando el número de
Prandtl, Pr es igual a la unidad. A partir del estudio que se hizo en la sección
28.4, se puedeampliar el postulado de Reynolds para incluir el mecanismo
correspondiente a la transferencia de masa cuando el número de Schmith, Sc,
también es igual a la unidad. Por ejemplo, si se analiza el flujo laminar sobre
una placa plana, donde S c = 1 , los perfiles de la velocidad y la concentración
dentro de las capas límite están relacionados por medio de:
a u,
"i
a y cc~A ,- c~A-, sc A , ~ ) l y = " = ~ ( ~ ) l y = " (28-36)
Recuérdese que en la frontera próxima a la placa, donde y = O , se puede ex-
presar el flujo de masa en función de la difusividad de la masa o del coeficiente
de transferencia de masa, por medio de la ecuación
a
NA,y = - D A R - ( C A - C A , s ) l y = O = kc(CA,s-CA,CC) (28-37)
ay
Se pueden combinarlas ecuaciones (28-36) y (28-37) para lograr una expresión
que relacione el coeficiente de transferencia de masa con elgradientede la
velocidad en la superficie:
(28-38)
El coeficiente de fricciónsuperficial se relacionó, en el capítulo 1 2 ,

con este mismo gradiente de velocidad, mediante la expresión:
( 12-2)
Si se hace uso de esta definición se puede reordenar la ecuación (28-38) para

obtener la analogía de Reynolds de transferencia de masa correspondiente a
los sistemas duyo número de Schmidt es igual a uno,
L C , "
(28-39)
Urn 2
*O. Reynolds, R o c . Manchester Lit. Phil. Soc., 8, (1874).

La ecuación (28-39) es análoga a la analogía de Reynolds de transferencia de

energía para los sistemas cuyo número de Prandtl es igual a uno. En el capí-
tulo 19 se analizó esta analogía, la cual puede expresarse así:
(19-36)
Los datos experimentales correspondientes a la transferencia de masa hacia

corrientes gaseosas, concuerdan, aproximadamente, con la ecuación (28-39)
cuando el sistema tiene un número de Schmidt cercano a 1 y si s u resistencia
al flujo se debe a la fricción superficial. Por lo tanto, no debe usarse la ecua-
ción (28-39)para describir situaciones en donde existaun arrastre de forma.
CONSIDERACIONES ACERCA DEL FLUJO TURBULENTO
En la mayor parte de las aplicaciones prácticas, el fl.ujo, en la corriente

principal, es turbulento en lugar de ser laminar. .Aunque muchos investiga-
dores han hecho grandes aportaciones a la comprensión del flujo turbulento,
hasta el momento nadie ha tenido éxito en la predicción de los coeficientes
de transferenciaconvectiva o factores defricción por medio del análisis directo,
lo cual no resulta demasiado sorprendente cuando se recuerda, en relación
con el estudio previo acerca del flujo turbulento, sección 13.1, que el flujo
en cualquier punto está sujeto a fluctuaciones irregulares en dirección y ve-
locidad. De acuerdo con esto, cualquiera de las partículas de fluido sufre una
serie de movimientos efectuados al azar, sobrepuestosal flujo principal. Estos
movimientos en forma de remolino, propician el mezclado en el centro tur-
bulento. A este proceso se le llama, a menudo, "difusión de remolino". El
valor de la difusividad de remolino de la masa seri mucho mayor que la di-
fusividad molecular que existe en el interior del núcleo turbulento.
En un esfuerzo para caracterizar este tipo de movimiento, Prandtl pro-
puso la hipótesis dela longitud de mezclado, tal como se estudió enel capítulo
13. En esta hipótesis, cualquier fluctuación de la velocidad, u%' se debe al
movimiento en la dirección de y de un remolino, en una distancia igual a la
longituddemezclado, L. El remolinodefluido,que posee una velocidad
media, üx(,,se desplaza hacia una corriente en la queel fluido adyacente tiene
una velocidad media, La fluctuaciónde la velocidad se relacionacon
el gradiente de la velocidad media, por medio de la ecuación
(13-10)
El esfuerzo cortante total se definió por medio dela expresión:
( 13-8)
Al sustituir la ecuación (13-10)en la (13-8) se obtiene:
(28-40)
6
d üx
T=p[v+€&f]- (28-4 1)
dY
donde = Lvy' se denomina difusividadderemolinode momento y es aná-
loga a la difusividad de momento,v.
Ahora se puede analizar, de manera semejante, la transferencia de masa
que ocurre en un flujo turbulento, ya que este mecanismo de transferencia
también se debe a la presencia de las fluctuaciones o remolinos. En la figura
28.6, la curva representa una porción del perfil turbulentolade concentración
con un flujo medio en la dirección de x. La rapidez instantánea de transfe-
rencia de la componenteA en la dirección de y es:
I '
NA.y = CA uy (28-42)
donde CA = FA + c A ' , o sea la suma de
las fluctuaciones instantáneasel promedio
temporal, de la concentración de la componente A . De nuevo puede utilizarse
el concepto dela longitud
Figura 28.6 Porción de la curva del perfil de la concentración turbulenta, mos-

trando la longitud del mezclado, de Prandtl.
de mezclado, para definir la fluctuación de la concentración por medio de la

relación:
(28-43)
AI insertar la ecuación (28-43) en la (28-42), ST obtiene una expresión que

corresponde a la transferencia turbulenta de masa por medio de remolinos.
La transferencia total de masa, normallaa dirección de flujo es:
6
(28-44)
"
donde ED = Luy' se denomina difusividad de remolino de la masa.

Por medio de un razonamiento semejante st: obtuvo, en el capítulo 19,
una expresión correspondiente a la transferencia convectiva de calor:
(19-49)
donde CY es la difusividad térmica molecular y e,,, , l a difusividad térmica de

remolino.
L a difusividad de remolino tiene un papel muy importante en muchos
procesos de transferencia de masa. Por ejemplo: existe una transferencia de
masa entre un fluido que circula alrededor de los sblidos de los reactores ca-
talíticos heterogéneos, los calefactores, los secadores, etc. Como resultadode
la difusión de remolino, la transferencia es rápida en el núcleo turbulento re-
duciendo la composición de cualquier gradiente. Al acercarse a la pared, la
turbulencia va disminuyendo hasta que, enla vecindad inmediata de la super-
ficie del sólido, desaparece esencialmente y la transferencia se realiza casi
completamente por difusión molecular. La mayor parte de la resistencia a la
transferencia ocurre en la capa límite próxima a. lasuperficieen lacual el
gradiente de la composición es más pronunciado.
LAS ANALOGIAS DE PRANDT Y VON KARMAN

En el capítulo 19 se obtuvo la analogía de Prandtl correspondiente a l a
transferencia de calorymomento,cuando se estudió el efecto del núcleo
turbulento y de l a subcapa laminar. Se puede utilizar el mismo razonamiento,
con respecto a la transferencia de masa y de momento, para desarrollar una
analogía semejante. En Ia subcapa laminar las difusividades de remolino de
momento y de masa son despreciables y en la superficie el esfuerzo cortante,
7s y el flujo de masa,NA,y,s, son constantes. La ecuación (28-41) se puede in-
tegrar sobre el espesor de la subcapa, resultando:
6
(28-45)
La ecuación (28-44) también se puede integrar sobre el espesor de la subcapa,

dando como resultado:
(28-46)
Si se elimina [ de estas dos ecuaciones, se obtendrá:
L a analogía de Reynolds, k c / v m= cf/2

= 7,/pUm2, se puede utilizar en el núcleo
turbulento, de y = .$ a y en las condiciones globales. El flujo de masa del
núcleo turbulento se transforma en:
(28-48)
Al eliminar a C A I se
~ obtiene, de las ecuaciones (28-46) y (28-48):
(28-49)
Cuando se sustituyen las ecuaciones que definen el proceso:
y también:
en la ecuación (28-49),se puede simplificar la relación hasta quedar así:
1 - 2
2[vcD+ v x p - 1)1
-"
kc cpm
o, escrita de manera ligeramente distinta:
(28-50)
Nótese que la ecuación (28-50) simplifica la analogía de Reynolds con la si-

guiente restricción: S c = l. En el capítulo 13 se definió la subcapa laminar
mediante la expresión: v+ = y + = 5, donde v+ = v x l t / ( v , q 2 ) ; por lo cual:
(28-5 1)
Cuando se sustituye el valor de Uxlt/vm en la ecuación (28-50) se obtiene una

analogía correspondiente a la transferencia convectiva de masa, semejante a
la analogía de Prandtl de la transferencia convectivade calor:
kc
-=
G /2
(28-52)
vm 1+ 5@(Sc - 1)
Al reordenar y multiplicar ambos lados de la ecuacilón (28-52) por vooL/D~*,

donde L es una longitud característica, se obtiene:
ó
(Cf/2jReSc
NUL,AB= (28-53)
1+5@(Sc- 1)
Las ecuaciones (28-52) y (28-53) son semejantes a la analogía de Prandtl de

transferencia de energía-momento, ecuación ( 19-57).
Von Kármán amplióla analogía de Prandtl, tomandoen cuenta la llamada
“capa amortiguadora” además de l a subcapa laminar y del núcleo turbulento.
Esto lo condujo a la formaciónde la analogía de uon Kcírmún:
(Cf/2jRePr (19-58)
Nu =
1 + 5JCf/2{Pr- 1 + In [( 1 + 5Pr)/6]}
que corresponde a la transferencia de energía y momento. El análisis de von

Kármán de la transferencia de masa, daorigen a las ecuaciones:
( Cf/2)ReSc
N u A ~= (28-54)
1+5~{S~-l+In[(1+5Sc)/6]}
ó
(28-55)
682 Transferencia convectiva demasa
Los resultados de la mayor parte de las analogías se pueden escribir en

la forma general, tal comola de las ecuaciones (28-52) y (28-55),en las que el
denominador del lado derecho forma un grupo complejo de términos, que
sirve de corrección a la analogía simple de Reynolds.
ANALOGIA DE CHILTON-COLBURN
Chilton y Colburn” descubrieron que la mejor correlación entre los datos

experimentales obtenidos en l a transferencia de masa, está dada por la ex-
presión:
(28-56)
Esta analogía es válida para h‘quidos y gases dentro delos valores 0.6<Sc<2500.
El término jD se denomina factor j de transferencia de masa y es análoga al
factor j de transferencia de calor, definido por la ecuación (19-39). Aunque
la ecuación (28-56) es una ecuacibn empírica basada en los datos recogidos
en los regimenes de flujo laminar y turbulento, se puede demostrar que sa-
tisface la solución exacta correspondiente a un flujo laminar sobre una placa
plana,
=0.332Re,”’Sc1’’ (28-20)
Si se dividen ambos lados de la ecuación entreRe,Sc”’, se obtiene:
(28-57)
Esta ecuación se reduce a la analogía de Chilton-Colburn cuando en la ex-

presión de arriba se sustituye la solución de Blasius correspondiente a la capa
laminar Iímite,
N ~ X , A BN k . 4 B S C 2 / 3 = C
1/3=-
Re, Sc Re, Sc 2
ó
*A. P. Colburn, Trans. A.I.Ch.E., 29, 174-210 (1933);T.H. Chilton y A. P. Colburn, Ind. Eng. Chem.,
26, 1183 (1934).
La analogía completa de Chilton-Colburn es:
(28-59)
la cual relaciona los tres tipos de transferencia en una sola expresión. La ecua-
ción (28-59) es exacta para placas planas y satisfactoria cuando se emplea en
sistemas cuya geometría sea diferente, a condici6n de que no exista ningún
arrastre de forma. Ensistemas enlos que este último tipo de arrastrese encuen-
tra presente, se ha encontrado que:
(28-61)
La ecuación (28-61) relaciona la transferencia convectiva de calor con la de

masa; permite la evaluación de un coeficiente desconocido de transferencia
a través de la información obtenida para otros fenómenos de transferencia. Es
válida para gases y líquidos, entre los valores: 0.6<Sc<2,500 y 0.6<Pr<100.
EJEMPLO 2
Fluye aire a presión atmosfkrica a través de un termómetro cuyo bulbose ha cubierto
con una tela húmeda. Este clásico termómetro de bulbo himedo indica que una pequeña
cantidad de líquido que se evapora en una gran cantidad de mezcla no saturada de vapor
y gas, alcanza U M temperatura de estado uniforme. La lectura del termómetro es de 65" F.
A esta temperatura se calcularon las propiedades siguientes:
presión de vapor del agua O .305 6 psi

densidad
aire del ,-
0.076 lb, /ft2
calor latente de vaporización del agua : 1057 Btu/lb,
número de Prandtl 0.72
número de Schmidt 0.61
calor específico, cp, del aire 0.24 Btu/lb," F
(Cuál es la tempertura del aire seco?

La ecuación (28-1) define el flujo molar del agua que se! evapora
NHzO = kc (cHKh -C H 2 0 , ~ )
La energía requerida para evaporar esta agua es suministrada por transferencia convectiva
de calor, o sea:
A
4=h(7',-T,)=AM~2JVH20
donde X es el calor latente de vaporización a la temperatura de la superficie. Esta ecuación

se puede resolver para la temperatura global:
Si se sustituye la ecuación ( 2 8 - 1 ) en esta ecuación, se obtiene:
kc
Tw = A M H ~ o ~ ( c -cHzo.m)
I - ~ ~ o ,+Ts
~
LOSfactores j de Chilton-Colburn nos dan una relación que corresponde a la razón k,lh:
. .
IH = I D
h
-((pr)Z/3 kc
= _(SC)*/'
PmCp
Cuando se sustituye esta expresión en la ecuación que se tiene, de la temperatura global,

se obtiene:
Las concentraciones son:
/0.3056 \
\ 14.7 atm) lb mole
cHzO.s = , atm . f t 3 = 5.42 X 10-~-
.73 )(525"R) ft'
lb mole OR
CH20.m =0
AI sustituir los valores conocidos, se obtiene:
(1057 Btu/lb,)(l8 lb,/lb mole) 0.72

T= (-)"'(5.42 x
(0.076 Ib,/ft3)(0.24 Btu/lb, OF) 0.61
T = 128.2"F (326 K)
28.7 M O D E L O S D E C O E F I C I E N T E S D E
T R A N S F E R E N C I A DE M A S A
Durante muchos años se han estado usando los coeficientes de transfe-

rencia convectiva de masa en e1 diseño de equipo de transferencia de masa.
Sin embargo, en muchos casos han sido coeficientes empíricos determinados
a partir de investigaciones experimentales. Una explicación teórica de los coe-
ficientes requeriría de una mejor comprensión de los mecanismos de la tur-
bulencia,yaqueestándirectamenterelacionadosconlascaracterísticas
dinámicas del flujo. En el capítulo 26, se introdujeron dos posibles modelos
para explicar la transferencia convectiva de masa. Tanto la teoría de película
como la de penetración se han aplicado ampliamente.
La teoría de película se basa en la presencia de una película ficticia de
fluido en flujo laminar, pr6xima a la frontera, que ofrece una resistencia a la
transferencia de masa, igual a la que existe verdaderamenre en todo el fluido
en movimiento. En otras palabras, se supone que toda la resistencia a la trans-
Modelos de coeficientes de transferncia de masa 685
ferencia existe en una película ficticia en la qule l a transferencia se realiza

totalmente por medio de difusión molecular. El espesor de la película, 6, debe
extenderse más allá de la subcapa laminar para illcluir la resistencia equiva-
lente que se encuentra cuando la concentración c,ambia dentro de la capa de
amortiguamiento y el núcleo turbulento. En cuanto a la difusión a través
de una capa que n o se está difundiendo, o fluido estancado, esta teoría pre-
dice que el coeficiente de transferencia de masa es:.
(26-9)
tal como sevio en el capítulo 26. El coeficiente de transferencia de masa,

correspondiente a la contradifusión equimolar, se expresó en la forma:
kco =-DAB (26-24)

6
En ambos casos, el coeficiente de transferencia convectiva de masa se rela-

ciona directamente con al difusividad de la masa molecular. Obviamente, el
espesor de la película ficticia 6, nunca puede med.irse puesto que no existe.
Por esto y porque aparentemente no es adecuada para explicar físicamente
la transferencia convectiva demasa, se han postulado otros modelos para
describir este fenómeno.
La teoría de la penetración fue expuestaoriginalmente,por Higbie*
para explicar la transferencia de masa en la fase líquida durante la absorción
del gas. Danckwertst laha aplicado al flujo turbulento y también lo han
hecho otros investigadores cuando la componente en difusión sólo penetra a
corta distancia enlafase de interés a causa de su desaparición rápida por
medio de unareacción química o del tiempo, re1ativa:mentecorto, de contacto.
Higbie supuso que la masa se transfería a la fase líquida por medio de la
transferencia molecular en estado no permanente. Con este concepto, se
expresó el flujo de masa en la interfase líquido-gas, en la forma:
(26-86)
Danckwerts aplicó este concepto de estado no permanente a la absorción de

la componente A en una corriente turbulenta de líquido. Su modelo supone
que el movimiento del líquido ocasiona constantemente remolinos de líquido
nuevo del interior hacia la superficie, donde desplazan a los elementos de
líquido que se encontraban anteriormente en la superficie. Mientras se en-
cuentra en la superficie cada uno de los elementos del líquido queda expuesto
a la segunda fase y se transfiere masa al líquido como si este estuviera estan-
cado y fuera infinitamente profundo. La rapidez dle transferencia depende
*R.Higbie, Trans. A. I. Ch. E., 31,368-389 (1935).

f P. V. Danckwerts, Ind. Eng. Chern., 43, 1460-67 (1951).
del tiempo de exposición. Se pueden hacer muchas suposiciones respecto a la

renovación de la superficie. Por ejemplo, cada uno de los elementos de lasu-
perficie puede tener el mismo tiempo de exposición antes deser reemplazado.
Esto permite suponer que la transferencia instantánea de masa ocurrirá de
acuerdo con la ecuación (26-86). La cantidad total de soluto que penetra en
el remolino en un tiempo de exposiciónigual a texp,es:
y la rapidez mediade transferencia durantela exposiciónse obtiene dividiendo

esta ecuación entre el tiempo de exposición:
(28-62)
Danckwerts modificó la suposición de un periodo constante de exposición,

proponiendo un conjunto “infinito” devalores de las edades de los elementos
de la superficie. Las funciones de distribución de edades en la superficie se
introdujeron para predecir la probabilidad de que un elemento de superficie
sea reemplazado por un remolino con elementos nuevos. La rapidez de reno-
vación superficial se creía constante en un cierto grado de turbulencia e igual
a un factor S de renovación superficial. La rapidez de transferencia de masa
con renovación superficial aleatoria es:
NA= G s ( c A , S -C A , ~ ) (28-63)
Actualmente los valores de S se obtienen por medio de investigaciones expe-
rimentales. El concepto de renovación superficial ha tenido mucho éxito en
la explicación y análisis de la transferencia convectiva de masa, particular-
mente cuando la transferencia de masa se presenta acompañada de reacciones
químicas en la fase líquida*?, pero se necesitan una considerable verificación
experimental y un enorme desarrollo para definir claramente este modelo.
En este texto nose hace un estudio detallado de los coeficientes de trans-
ferencia de masa correspondientes a los sistemas quereaccionanquímica-
mente. En el estudio que se realizó en la sección 26.2 acercade la transferencia
de masa molecular asociado con una reacción química, se demostró que Ia
transferencia de masadepende de la constantederapidezde la reacción
química. Es de esperarse que exista una dependencia similar correspondiente
al coeficiente de transferenciaconvectiva de masa. Existen estudios excelentes
acerca de este tema*$.
Toor y R4archello$ señalaron que es válido el concepto de penetración
de Danckwerts solamente cuando la renovación de la superficie es relativa-
*P. V. Danckwerts, Gas-Liquid Reactions, McGraw Hill Book Co., New York, 1970).
G. Astarita, Mas 7ransfer with Chemical Reaction, Elsevier Publishing Co., Amsterdam, 1967.
$ H. I,. Toor y J. M. Marchello,A. I. Ch. E. J. 1 , 97 (1958).
Conclusión 687
mente rápida, proveyendo, así, la a superficie de elementos jóvenes sobre una

base continua. Para los elementos de mayor edad de la superficie, se establece
un gradiente de concentración en estado permanente, tal como lo predice la
teoría de película. Por esto, la rapidez de transferencia de masa debe ser
directamente proporcional a la difusividad molecdar de masa. Para valores
pequeños del número de Schmidtel gradiente perm.anente de la concentración
se crea muy rápidamente en cualquier elemento nuevo de superficie de tal
modo que, a menos que la rapidez de renovaciónsea lo suficientemente grande
para hacer desaparecer la enorme fracción de los elementos de la superficie
antes de que se penetre en ellos, la mayor parte de los elementos de la superfi-
cie se comportarán como elementosviejos. AI aumentar el número de Schmidt,
el tiempo necesario para provocar un gradiente constante, aumenta rápida-
mente y por esto, es suficiente una rapidez de renovación superficial pequeña
para evitar que se penetre a la mayoría delos elementos, tanto jóvenes como
viejos, las características de transferencia son intermedias entre los modelos
de película y penetración. Los coeficientes de transferencia convectiva de
masa son proporcionales a una potencia de la difusividad molecular de lamasa
que se encuentra entre los valores 0.5 y 1.0. Estas conclusiones están basadas
en datos experimentales.
En ambos modelos: de película y de penetración, la transferencia de
masa comprende una interfase entre dos fluidos en movimiento. Cuando una
de las fases es un sólido, la velocidad del fluido paralela a la superficie en la
interfase debe ser igual a cero, por lo cual es de esperarse que exista la nece-
sidad de u r ~tercer modelo, el modelo de la capa límite para correlacionarlos
datos correspondientes a. un sólido que se sublima convirtiéndose en gas o un
sólido que se disuelve, transformándose en liquido. El coeficiente promedio
de transferencia de masa correspondiente a la difusión a través de una capa
laminar límite, es:
*"""
k, = 0.664-Re,"2Sc''3
DAB (28-21)
L
Esto demuestra que el coeficiente de transferencia de masa varía con
forma típica de los c~lculosde capa límite.
En este capítulo se estudiaron los principios de transferencia de masa

por convección forzada, los parámetros importantes que ayudan en la des-
cripción de la transferencia convectiva de masa y los rnodelos propuestos para
explicar el mecanismo de transferencia convectiva. Se ha visto que la trans-
ferencia de masa por convección está íntimamente relacionada con las carac-
terísticas dinámicas del fluido en movimiento, particularmente con el fluido
que se encuentra en la vecindad de la frontera. A causa de las grandes seme-
janzas entre los mecanismos de transferencia de momento, energía y masa,

que vimos que se pudieron utilizar los mismos cuatro métodos en la evaluación
de los coeficientes de transferencia convectiva de masa que se encontraron
originalmente con el fin de analizar los coeficientes de transferencia convec-
tiva de calor. En los cuatro análisis se correlacionó el coeficiente de trans-
ferencia convectiva de masa por medio de la ecuacióngeneral:
NuAB = f(Re, Sc)
La transferencia de masa hacia las corrientes turbulentas, tambiénse estudió
y se definió la difusividad de remolino de la masa. Se presentaron analogías
correspondientes a la transferencia convectiva de masa hacia las corrientes
turbulentas.
PROBLEMAS
28.1 En una columna de atomización de transferencia de masa, un líquido

se atomiza hacia una corriente de gas y hay un intercambio de masa
entre las fases líquida y gaseosa. La masa de las gotas formadas en el
atomizador se considera que es función del diámetro de la tobera, de
la aceleración de la gravedad, de la tensiónsuperficial del líquido contra
el gas, de la densidad del fluido, de la viscosidad y la velocidad y de la
viscosidad y la densidad del medio gaseoso. Ordene estas variables en
grupos adimensionales. iCree usted que debían haberse incluido algunas
otras variables?
28.2 Un cilindro largo de barro poroso, cuya concentración inicial de agua
era de cA0, se inserta repentinamente en una corriente de aire que tiene
un contenido de humedad de cA,m. Si el radio del cilindro es ro y el
coeficiente promedio de transferencia demasa del cilindro a la corriente
de aire es k c , demuestre, usando el análisis dimensional, que el perfil de
la concentración dentro del cilindro debe expresarse en función de los
parámetros:
CA(r)-CA,m - k t
Y 2
CAo- CA,m r0 r0
28.3 La transferencia de masa en una capa límite turbulenta formada sobre

una placa plana, está definida en función del número local de Nusselt,
por medio de:
N u A B , ~= 0.0292 Re,4'SS~1'3
Si la transición a flujo turbulento ocurre en x = 0.6 m y la placa es de
1.5 m de longitud y 1 m de ancho, obtenga las ecuaciones correspon-
dientes a:
u ) el coeficiente medio de transferencia de masa k;, correspondiente a
toda la placa,
b ) la rapidez de transferencia de masa, W, .

Lectura Complementaria 2 PDF

Cargado por

Copyright:

Formatos disponibles

Lectura Complementaria 2 PDF

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

Lectura Complementaria 2 PDF

Cargado por

Copyright:

Formatos disponibles

FUNDAMENTOS DE LA

LOS capítulos anteriores relacionados con los fenómenos de transferencia

La masa puede transferirse por medio del movimiento molecular fortuito

Ya en el aiio de 18 15 Parrot observi, cuantitativamentequecuando

Con la existencia del gradiente de concentración, habrá más moléculas

LA ECUACION DE RAPIDEZ DE FICK

Las leyes de transferencia demasa ponen de manifiesto la relacibn entre

Figura 24.1 Volumenelemental que contieneunamezcladecomponentesmúltiples.

donde n es el número de especies presentes en la mezcla.La fracción de masa,

Por definción, la suma de las fracciones de masa, debe ser uno.

La concentración molar de la especie A , C A , se define como el número

dende MA es el peso molecular de la especie A.

o, en una mezclagaseosa que obedezca la ley de los gases ideales:

Cuando una mezcla obedece la ley de

La ecuación (24-8) es una representación algebraica de la ley de Dalton

En la Tabla 24.1 aparece un sumario de los diferentes términos de con-

oxigeno, O,, yo, = 0.2 1

Determínese la fracción de masa, tanto de oxígeno como de nitrógeno y el peso

oxígeno presente= (1 mo1)(0.21) = 0.21 mol

TABLA 24.1 CONCENTRACIONES EN UNA MEZCLA BINARIADE A Y B

p = Densidad total demasa de la mezcla

wg = Fracción de masa de la especie B = P B / P

Rlezcla líquida o sólida Mezcla gaseosa

c=Densidad molar dela mezcla= n / V c=n/V=P/RT

nitrógeno presente = (1 mo1)(0.79)= 0.79 mol

= (0.7Y mol) kg) = 0.0221 kg

El volumen de la mezcla gaseosa a 298K3,es

- 0.79 mol = 32.3 mol N,

La densidad total, p , es:

y el peso molecular medio de la mezcla es:

M = -p= 1.180 kg/m' = 0.0288 kg/mol

velocidades. Enunsistemadecomponentesmúltiples, las diferentes

donde vi denota la velocidad absoluta de la especie i con relación a ejes esta-

Esta es l a velocidad que mediría un tubo de Pitot y la que se encontró

La velocidad de una especie particular con relación a la masa promedio

ui-u, que es la velocidad de difusihn de l a especie i con relación a la velo-

vi-V. que es l a velocidad de difusión de l a especie i con relación a l a velo-

Flujos. El flujo de masa ( o molar) de una especie dada es una cantidad

donde , / A , G es el flujo molar en la direccih de z relativa a l a velocidad molar

DA E , el factor de proporcionalidad, es la difusividad d e la masa o coeficiente

Como la concentración total, c es constante bajo condiciones isotérmicas

donde d w , / d z es el gradiente de concentración en funcihn de la fracción de

En un sistema binario con una velocidad media constante en la direc-

JA,z = C A (VA.=- Vz) (24-18)

Si se igualan las expresiones (24-16) y (24-18 ) , se obtiene:

AI sustituir esta expresión enla relación que se tenía, se obtiene:

Como las velocidades componentes, V A y U B , ~ son , velocidadesrelati-

AI sustituir estos símbolos en la ecuación (24-19), se obtiene una rela-

Esta relación se puede generalizar y escribir en forma vectorial de la ma-

N A =-cDAB V Y A + Y A (NA+ NB) (24-21)

Es importante notar que el flujo molar NA es la resultante de las dos

DAB VYA molar, JA , que

donde D AM es el coeficiente de difusión de A de la mezcla.

En condiciones isotérmicas e isobáricas, esta relaciónse reduce a:

'TABLA 24.2 FORMAS EQUIVALENTES DE LA ECUACION DE FLUJO

De acuerdo con la segunda ley de la termodinámica, los sistemas que nose

donde uA es la "movilidad" de la componente A o la velocidad resultante de

La ecuación (24-24) se puede usar para definir todos los fenómenos de

dende p o es una constante: el potencial químico del estado estándar. Cuando