LECTURA COMPLEMENTARIA 2 Transferencia de Materia

FUNDAMENTOS DE LA
TRANSFERENCIA DE MASA
LOS captulos anteriores relacionados con los fenmenos 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 ms componentes cuyas concentra-
ciones varan de un punto a otro, presenta una tendencia natural a transferir
la masa, haciendo mnimas las diferencias de concentracin dentro del siste-
ma.
La transferencia de un constituyente de una regin de alta concentraci6n
a una de baja concentracin se llama transferencia d e masa.
Muchas de nuestrasexperienciascotidianasestnrelacionadascon la
transferencia de masa. Un terrn de azcar en una taza de caf negro se di-
suelve finalmente y despus 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 atmsfera 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 difusin de neutrones den-
tro delos reactores nucleares, la difusin de sustancias quelos poros del carbbn
activado absorben, la rapidez de las reacciones qumicas catalizadas y biolgi-
cas as como el acondicionamiento del aire, son ejemplos tpicos.
Si se piensa en el terrn de azcar puesto dcntro de l a taza de caf ne-
gro, la experiencia nos ensea queel intervalo de tiempo que se requiere para
destribuir el azcar depende desi el lquido 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 dinmica 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 caractersticas dinmicas del flujo.
Estos dos modos diferentes de transferencia de masa: molecular y con-
vectiva, son anlogos a la conduccibn calorfica 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 actan, a menudo simultnea-
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 ms especies moleculares cuyas concentra-
ciones relativas varan de un punto a otro, resulta un proceso, aparentemente
natural, que tiende a disminuir cualesquiera desigualdades de composicibn.
Estatransferenciamacroscpicademasa,independientedecualquier
conveccin que se lleve a cabo dentro del sistema, se define con el nombre de
difusin melocular.
En el caso especfico delas mezclas gaseosas se puede deducir una expli-
caci6n l6gica de este fenmeno de transferencia a partir de la teora cintica
de los gases.
A temperaturas superiores al cero absoluto, las molculas individuales se
encuentran en un estado de movimiento continuo, aunque fortuito. Dentro
de las mezclas de gases diluidos, cada una de las molculas de soluto se com-
portaenformaindependientede las otrasmolculasdesoluto,yaque
rara vez se topa con ellas. Estin ocurriendo contlnuamente colisiones entre
el solvente y el soluto. Como resultado de estas colisiones, las molculas del
soluto describen trayectorias en zigzag, a veces hacia una regibn de mayor
concentracin, a veces hacia una de concentracin m i s baja.
Examinemos una seccibn hipotktica que pase en forma normal al gra-
diente de concentracin dentro de una mezcla gaseosa isobrica e istrmica
que contenga molculas de solvente y soluto.
I,os dos elementos delgados e iguales de volumen que se encuentran so-
bre y por debajo de la seccin, contienen el mismo nmero de molculas, tal
Como lo estipula la ley de Avogadro.
Aunque n o es posible establecer la direccin especfica en la queviajar
una mol6cula particular en un intervalo dado de tiempo, puede decirse que
un nmero definido de molculas que se encuentren enel elemento inferior de
volumen cruzar la seccin hipottica desde abajo y el mismo nimero de mo-
lculas abandonar6 el elemento superior y atravesar la seccin desde arriba.
Transferencia demasa molecular 535
Con la existencia del gradiente de concentracin, habr ms molculas

de soluto en uno delos elementos de volumen que en el otro; as pues, resul-
tar una transferencia total neta de unaregin de concentracihn nlayor a una
de cocentracin menor.
El flujo de cadaunade las especiesmolescularesocurreen la direc-
cin del gradiente negativo de la concentracibn.
Como se seal en los Captulos 7 y 15, la transferencia molecular de
momento y la transferencia de energa por conduccin se deben tambin al
movimiento al azar de las molculas. De manera que es de esperarse que los
tres fenbmenos de transferencia dependan de muchas propiedades y caracte-
rsticas iguales, tales como la de trayectoria media libre y el hecho de que los
anlisis tehricos de los tres fenmenos tienen muchoen comn.
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 concentracin
responsable de esta transferencia de masa. Desafortunadamente, la descrip-
cin cuantitativa de la difusin molecular es considerablemerlte mis comple-
ja que las dcscripciones anlogas correspondientes a la transferencia molecular
de momento y energa que tienen lugar en una fase de una componente.
Como la transferencia de masa o difusin, como se le llama tambin,
ocurre solamente en mezclas, su evaluacin debe inclulr un examen del efec-
to de todaslas componentes. Por ejemplo, a menudo se desea conocer la rapi-
dez de difusin de una componente especfica 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 comn 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 mltiple de componentes la concentra-
cin 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 molculas
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 contieneunamezcladecomponentesmltiples.
donde n es el nmero de especies presentes en la mezcla.La fraccin de masa,

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

n
w,=1 (24-3)
t=l
La concentracin molar de la especie A , C A , se define como el nmero

de moles de A,presentes por unidad de volumen dela mezcla. Por definicin,
un mol de cualquier especie contiene una masa equivalente a su peso mo-
lecular.. Los trminos de l a concentracin de masa y de l a concentracih molar
estn relacionados por medio de la siguiente expresin:
(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 trminos de las presiones parciales. Bajo las condiciones
en las cuales es vlida la ley de los gases ideales, PA V=nA RT, la concentra-
cin molar es:
nA PA
C A =-=- (24-5)
V RT
donde PA es la presin parcial de la especie A en la mezcla, n A es el nmero
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 concentracin molar total,c , es el nmero 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 presin total. La fraccin molar correspon-
diente a las mezclas de lquidos o slidos, 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 = - (lquidos y slidos)
C
(24-7)
CA
yA = - (gases)
C
Cuando una mezcla obedece la ley de

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

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

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

nitrogeno, N2, yN2= 0.79
Determnese la fraccin de masa, tanto de oxgeno como de nitrgeno 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 oxgeno es de 0.032 kg/mol y el del nitrgeno es de 0.028 kg/mol.
Como base para nuestros clculos, tmese 1 mol de mezcla gaseosa;
oxgeno 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
Concentracin 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, = Fraccin de masa de la especie A = p a / p
wg = Fraccin de masa de la especie B = P B / P

P=PA+PB
1=6J,+w,
Concentraciones molares
Rlezcla lquida o slida 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)
nitrgeno 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 nme-
ros redondos, como 0.029 kglmol.
Este problema podra, tambin resolverse utilizando la ley de los gases ideales PV=
nRT. En condicionesideales, OC o 273 K y 1 atm o 1.O13X1O5 Pa de presin 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. Enunsistemadecomponentesmltiples, las diferentes

especies se movern 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 mltiples 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 relacin a ejes esta-

cionarios de coordenadas.
540 Fundamentos de la transferencia de masa
Esta es l a velocidad que medira 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 mltiples, se define
en funcin de las concentraciones molares de todas las componentes, por me-
dio de la expresi6n:
(24- 14)
La velocidad de una especie particular con relacin a la masa promedio

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

cidad media de l a masa.
vi-V. que es l a velocidad de difusin de l a especie i con relacin 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
concentracin.
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 refirindose a las coordenadas que permanecen
fijas en el espacio, a las coordenadas que se estn moviendo con la velocidad
promedio de l a masa o con la velocidad molar promedio.
La relacibn bsica correspondiente a l a difusin molecular define el flu-
jo molar relativo a l a velocidad molar media, la. l"ick* fue quien primero pos-
tul una relacin emprica 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 concentracin 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 difusin correspondiente a una componenteA quese difunde a travs de la
componente B.
Groot* propuso una relacin ms general de flujo que no est restringi-
da a sistemas isotrmicoso isobricos.
Groot escribi:
concentracin difusin
) ( gradiente de )
o
Como la concentracin total, c es constante bajo condiciones isotrmicas

e isobricas, la ecuacin (24-15) es una forma especial de la relacin (24-16),
que es ms general. Una expresin equivalente, que corresponde ajA .z,quees
el flujo d e masa en la direccin de z , relativo a la velocidad promedio de la
masa, es:
(24-17)
donde d w , / d z es el gradiente de concentracin en funcihn de la fraccin de

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

cin de z , el flujo molar en la direccin de z , relativo a la velocidad molar
media tambin 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 ecuacin
(24-14), en la forma:
AI sustituir esta expresin enla relacin que se tena, 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 relacin 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 smbolos en la ecuacin (24-19), se obtiene una rela-

cin que corresponde al flujo dela componente A , relativa al eje z .
(24-20)
Esta relacin 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
concentracin. Este trmino se llama contribucin .
del gradiente dela concentracin.
Y
Y A ( N A +N,)=,, V El flujomolarqueresultacuando la componente
molar A circula con el flujo global. Este trmino
delflujo se llama contribucibndelmovimiento
global.
Cualquiera de estas cantidades puede ser una parte importante del flu-
j o molar total,N, . Cuando la ecuacin (24-21) se usa para describir una difu-
sin molar, la naturaleza vectorial de los flujos individuales N, y N,, se debe
analizar y despus debe evaluarse l a direccin de cada una de las dos cantida-
des vectoriales.
Si la especie A se estuviera difundiendo en una mezcla de componentes
mltiples, la expresin equivalente a la ecuacin (24-21), sera:
donde D AM es el coeficiente de difusin 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 funcin de la densidad de masa y
d i la fraccin 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 isotrmicas e isobricas, esta relacinse reduce a:
Las cuatro ecuaciones que definen los flujos JA ,j,, N, y nA son enun-
ciados equivalentes de la ecuacin de Fick de la rapidez.
El coeficiente de difusin DAB ,es idntico en todas.las ecuaciones, cual-
quiera de las cuales es adecuada para describir la difusin molecular; sin em-
bargo, ciertos flujos son ms fciles de utilizar en casos especficos. Los flujos
de masa, n, y j, se usan cuando tambin se requiere que las ecuaciones de
Navier Stokes describan el proceso. Ya que las reacciones qumicas se descri-
ben en funcin 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 qumicas. 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 difusin em-
pleadas para medir el coeficiente de difusin.
La Tabla 24.2 resume las formas equivalentes de la ecuacin 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 termodinmica, los sistemas que nose

encuentran en equilibrio tienden a ste despus de algn tiempo. Existe una
fuerza motriz generalizada, que, en trminos de qumica termodinmica es
-dpc/dz donde ,uc es el potencial qumico. La velocidad de difusin molar de
la componente A se define en funcin del potencial qumico por medio de:
(24-23)
donde uA es la "movilidad" de la componente A o la velocidad resultante de

la molcula mientras se encuentra bajo la influencia de una fuerza motriz
unitaria.
La ecuacin (24-23) se conoce como la relacin de Nerst-Einstein. El
flujo molar deA se transforma en:
(24-24)
La ecuacin (24-24) se puede usar para definir todos los fenmenos de

transferencia de masa molecular. Como ejemplo, estdiense las condiciones es-
pecificadas para la ecuacin (24-15); el potencial qumico de una componen-
te en una solucin ideal homognea a temperatura y presin constantes, se
define por medio de:
@, = E*.'+R T In cA (24-25)
dende p o es una constante: el potencial qumico del estado estndar. Cuando

se sustituye esta relacin en la ecuacin (24-24), se obtiene la ecuacin de
Fick de la rapidez correspondiente a una fase homognea:
(24-15)
Hay otrascondicionesfsicasadems de las diferenciasdeconcentra-

cin,queproducenungradiente de potencialqumico: las diferenciasde
temperatura y de presin as como las diferencias en las fuerzas creadas por
campos externos, tales como el gravitacional, el magntico y el elctrico. Por
ejemplo, se puede obtener la transferencia de masa aplicando un gradiente de
temperatura a un sistema de componentes mltiples;este fenmeno de trans-
ferencia, el efecto Soret de difusin trmica, aunque usualmente pequeo en
relacin con otros efectos de difusin, se usa con xito en la separacin de
istopos. Las componentes de una mezcla de lquidos se pueden separar me-
diente una fuerza centrfuga pordifusin de presin..
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-
cinporsedimentacinbajo la influencia de la gravedad,laprecipitacin
electroltica debida a un campo electrosttico de fuerzas y la separacin mag-
ntica de mezclas minerales a travs de la accin de un campo magntico de
fuerzas. Aunque estos fenmenos de transferencia de masa son importantes,
son procesos muy especficos.
La transferenciademasamolecularqueresulta de las diferencias de
concentracin y a la cual describe la ley de Fick, es el resultado del movi-
miento molecular fortuito en pequeas trayectorias medias libres, indepen-
dientes de las paredes del recipiente. La difusin de los neutrones rpidos y
de las molculas en poros pequefisimoso a una densidad muy baja degas no
pueden describirse mediante esta relacin.
Los neutrones producidos en un proceso de fisin nuclear poseen, ini-
cialmente, energas cinticas altas yse llaman neutrones rpidos a causa de SUS
altas velocidades, hasta de 15 millones de metros por segundo. Con altas ve-
locidades, los neutrones pasan a travs de las cubiertas electrnicas de otros
tomos o molculas con muy poco problema. Para deflectarse, los neutrones
rpidosdebenhacercolisinconunncleo,elcualrepresentaunblanco
muy pequeo comparado con el volumen de la mayora de los tomos y mo-
lculas. La trayectoria media libre de los neutrones rpidos es aproximada-
mente un milln de veces mayor que las trayectorias medias libres de 10s ga-
ses a presiones ordinarias. Despus de que los neutrones rpidos reducen su
velocidad por medio de colisiones elsticas de dispersin entre 10s neutrones
Y el ncleo del moderador del reactor. Estos neutrones ms lentos, quese Ila-
man neutrones trmicos, emigran de las posiciones de alta concentracin a
posiciones de ms baja concentracin y la ley de Fick de la difusin describe
este movimiento.
Si la cantidad de gas es baja, o si los poros a travs de los que viaja el gas
son muy pequeos, las molculas chocarn con las paredes con mayor fre-
cuencia que entre s. Esto se conoce con el nombre de flujo de Knudsen o di-
fusin d e Knudsen. Despuks dechocarcontralapared, las molculas son
absorbidas, momentneamente y despus 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 difusin de Knudsen en los gases median-
te el uso de la teora cintica. La constante de proporcionalidad que relaciona
el flujo de Knudsen con el gradiente de concentracin es independiente de la
presin. Se demostrar que la constante de proporcionalidad correspondiente
a la difusin de gases, tal como la describe la ley de Fick, es inversamente
proporcional a la presin.
La transferencia de masa descrita por la ley de Fick, que resulta de las
diferencias de concentracin es el proceso principal con el que se topan los
ingenieros as como el fenmeno de transferencia de masa que se estudiar
en captulos siguientes. Es importante, sin embargo, percatarse de que estos
otros fenmenos 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 difusin. Sus dimensiones fundamentales, que pueden obte-
nerse a partir de la ecuacin(24-15),
sonidnticasa las dimensiones fundamentales de las otras propiedades de

transferencia: la viscosidad cinemtica, v y la difusividad trmica, 01, o su ra-
zn 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-
gls, 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 difusin depende de la presin, de la temperatura y de

la composicin del sistema.En el Apndice de Tablas, Tablas J 1,J 2 y 53,pue-
den consultase los valores experimentales correspondientes a las difusivida-
El coeficiente de difusin 547
des de los gases lquidos y slidos, respectivamente. Como es de esperarse, de

acuerdo con la movilidad de las molculas, los coeficientes de difusin son
generalmente mayores en relacin con los gases (entre los valores de 5X lop6
y 1X m' /seg.), queenrelacincon los lquidos (entre losvalores: 10"'
y lo-' m' /seg.), que son mayores a los valores obtenidos en relacin con los
slidos (entre y .lo"' m' /seg.).
En ausencia de datos experimentales, se han obtenido expresionessemi-
tericas que aportan aproximaciones cuya validez es tan buena como la de
los valores experimentales debido a las dificultades que existen para la medi-
cin de estas ltimas.
DIFUSIVIDAD DE LA MASA GASEOSA
Las expresiones toricas correspondientes al coeficiente de difusin en

las mezclas gaseosas de baja densidad, en funcin de las propiedades molecu-
laresdelsistema, fueronobtenidasporJeans",Chapman?ySutherlandS
usando la teora cintica de los gases.Si se utiliza el razonamiento de estos
cientficos, para explicar los fenmenos de transferencia, se puede examinar
el movimiento de las molculas de gas, tal como se hizo en las secciones 7.3 y
15.2, y despus obtener una expresin que relacione el coeficiente de difu-
sin con las propiedades de los sistemas gaseosos.
Vase ahora el volumen de control de la figura 24.2. Si se especifica
que el gas es esttico o que fluye en forma laminar en la direccin de x, se
puede considerar que la transferencia de masa de la especie A en la direccin
de y se lleva a cabo solamente aescala molecular. Si se aplica la ecuacin(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 travs de la cara superior del elemento,
Esta ecuacibn establece sencillamente, que el flujo ascendente de masa

debe ser igual al flujo descendente.
Como primera aproximacin estudiemos un sistema que contenga mol-
culas de igual tamao y masa, que posea velocidades mediasiguales. Slo una
mezcla gaseosa formada por istopos del mismo elemento se parecera apro-
ximadamente a este sistema. Si se vuelve a examinar la ecuacin obtenida so-
bre una vase microscpica, se puede concluir que el nmero de moIculas que
atraviesa la cara superior desde abajo deber ser igual al nmero de molculas
que atraviesa esta cara desde arriba. Ya que existe una concentracin dela es-
pecie A , tal como aparece en la figura 24.2, se transportarn ms molculas a
travgs de la superficie de control desde arriba, que desde abajo. Esto da como
resultado un flujo neto de molculas deA en la direccin de y .
Tal como se hizo anteriormente en los Captulos 7 y 15, se usarn las si-
guientes ecuaciones, obtenidas de la teora cintica de los gases de baja densi-
dad:
donde? es la velocidad molecular fortuita media,?/4 es la velocidad de una

molcula individual mientras pasa a travs del rea A x Az, A es el paso me-
dio libre, K es la constante de Boltzman, m es la masa de una molcula, d es
el dimetro de las molculas esfricas, N es la concentra.cin molecular y Z es
la frecuencia con la que llegar una molcula al rea Ax Az. La ecuacibn de
continuidad, descrita en funcibn delas molculas en movimiento, es:
c c
m,- Ax Azl,--C m,- Ax Azly+= O
n 4 n 4
o, al sumar las N molculas en un volumen unitario:

Si van a tomarse en cuenta solamente las mcdculas de A que atraviesan

esta superficie, la ecuacin debe ser vlida para un flujo neto de masa en la
direccin de y
(24-27)
De nuevo, como en los Captulos 7 y 15, se supondr que el perfil de la

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

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

mezcla de molculas semejantes, o sea, A y el istopo A * , es:
1 - (2430)
DAA* =3 C
A
Este coeficiente de difusin se denomina, a menudo, coeficiente de au-
todifusibn y se le utiliza para explicar la difusin de las molculas trazadoras.
La sustitucin de los resultados cinticos correspondientes a y A, en la c
ecuacin (24-30), da:
(24-31)
En un gas ideal se puede reemplazar a N , usando la relacin
NKT=cRT= P
y obtenindose:
La ecuacin (24-32) demuestra que el coeficiente de difusin se puede

expresar totalmente en funcin delas propiedades del gas. De manera distinta
a lo que ocurre en otros coeficientes de transferencia de dos molculas,l a vis-
cosidad y la conductividad trmica, el coeficiente de difusin depende de la
presin as como de un orden ms 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 teora cintica han intentado tomar en
cuenta las fuerzas deatraccin y repulsinexistentesentre las molculas.
Nirschfelder, Bird y Spotz", utilizando el potencial de Lennard Jones para
evaluar la influencia de las fuerzas intermoleculares, encontraron una ecua-
cin adecuada al coeficiente de difusin correspondiente a parejas gaseosas
de molculas no polares? y no reactivas.
(24-33)
/"donde DA, es la difusividad de la masa de A , que se difunde a travs de B ,

,/ en /seg., T es latemperaturaabsolutaen K;M,, M, son los pesos mo-
leculares de A y B , respectivamente; P es la presin absoluta en atmsferas;
*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, colsltense F. Daniels y R. A. Al-
berty, Physical Chemistry, Wiley, Nueva York, 1955.
IEl coeficiente de difusin 551
u A B es el dimetro de colisin, que es un parlmetro de Lennard-Jones en

Angstroms; R, es la integral de colisin correspondiente a la difusin mo-
lecular, que es una funcin adimensional de la temperatura y del campo po-
tencial intermolecular correspondiente a una molcula de A y una de B. En
la Tabla K 1 del apndice, aparece una lista de valores de a, en funcin de
; k es la constante de Boltzmann,igual a 1.38 X 10 ergios /I( y E
es la energa de la interaccin molecular que corresponde al sistema binarlo
4B
AB,un parmetro de Lennard-Jones en ergios; ver ecuacin (24-40). La ecua-
cin (24-33) tambin seala la dependencia de la! presin y la temperatura,
tal como se demostr en la ecuacin (24-32).Como la ecuacin (24-33) defi-
ne la propiedad de transferencia de una mezcla gaseosa binaria, no es sor-
prendente encontrar que el coeficiente de difusin tambin 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 lacomposicinen la ecuacin(2.4-32) ni en las ecuaciones
similares correspondientesa la viscosidad y a laconductividadtrmica.
Los parmetrosdeLennard-Jones, u y , se obtienenusualmentede
datos acerca de la viscosidad. Desafortunadamente esta informacin puede
obtenerse solamente para unos cuantos gases puros. La Tabla K2 Apndice
tabula estos valores. En la ausencia de datos experimentales,los valores de las
componentes puras se pueden calcular a partir de las siguientesrelaciones
empricas :
u = 1.18 vp3 (24-34)

u = 0.841 Vc13 (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 normaldeebullicin,en

(cm)3 /g mol (esto se calcula utilizando la Tabla 214.3); V, es el volumen mo-
/ g T, es la temperaturacrticaen
lecular crticoen ( ~ m ) ~mol; K; es la
temperatura normal de ebullicin en K y P, es a l presin crtica, en atm-s-
feras.
Para un sistema binario compuesto de pareja:s moleculares no polares, se
pueden combinar los parmetros de Lennar-Jones de la componente pura,
emperricamente, por medio de las siguientes relaciones:
(24-39)
Y
EAR = JEAER (24-40)
Estas relaciones debern modificarse para parejas moleculares: polar-po-

lar y polar-no polar. Hirschfelder, Curtiss y Bierd* analizaron las modifica-
ciones propuestas.
La ecuacin de Hirschfelder, Bierd y Spotz, ecuacin (24-33), parece
ser la mejor ecuacin de que se dispone actualmente para el clculo 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 mximo de un 6% con
relacin a los valores medidos.
Como los coeficientes de difusin son difciles de medir con exactitud,
esta desviacin se puede deber, parcialmente, a los errores experimentales.
Otras ecuaciones empricas sugeridas por Arnoldt y Gillilandt y por Slattery
y Bird 5 , no concuerdan tanto.
De acuerdo conlo anterior, se recomienda utilizar la ecuacin deHirsch-
felder cuando no pueda hacerse uso de los valores experimentales.
La ecuacin de Hirschfelder, que es la (24-33), se utiliza a menudo para
extrapolar los datos experimentales. Para valores moderados de la presin,
hasta de 25 atmsferas, el coeficiente de difusin vara inversamente a la pre-
sin. Las presiones ms altas requieren, aparentemente, correcciones corres-
pondientes alos gases, densos; desafortunadamente no existe ninguna relacin
satisfactoria para las presiones altas. La ecuacion (24-33) tambin establece
que el coeficientededifusinvaracon la temperatura,deacuerdocon
l a relacin: ? I 2 laD.Si se simplifica la ecuacin (24-33), se puede predecir
el coeficiente de difusin a cualquier temperatura y a cualquier presin me-
nor de 25 atmsferas, a partir de un valor experimental conocido, por medio
de :
(24-41)
En la Tabla J 1 del Apndice, 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 eucacin (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)
Calclese el coeficiente de difusin del bixido de carbono en aire a 20C y presin

atmosfrica. Comprese este valor con el valor experimental que aparece en la Tabla J 1
del Apndice. Los valores de (T y E/K se obtienen de la Tabla. K2 del Apndice.
u, inA EA/K, in K
Bixido de carbono 3.996 190
aire 3.167 97
Los diversos parmetros que se necesitan para resolver la ecuacin (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 ecuacin, 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 atmsfera, $;etiene:
2
DAB= O. 136 cm/s

La ecuacin (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 (clculo anterior)

El valor corregido correspondiente a l coeficiente de difusin 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 pequea. Por esto, la mayora de los valores de las difusivi-
dades con relacin a la temperatura slo incluyen la razn ( T I / T),3 2 .
La transferencia de masa en las mezclas gaseosas de varias componentes
puede describirse por medio de ecuaciones tericas que incluyen los coefi-
cientes de difusin Correspondientes a las diversas parejas binarias que for-
man la mezcla. Hirschfelder, Curtiss y Bird*, elaboraron una expresin en su
forma ms general. Wilkei simplific la teora y demostr que una buena
aproximacin a la formacorrecta, es la dada por la relacin:
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-
cin 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
Determnese la difusividad del monxido de carbono en una mezcla de oxgeno 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 atmsferas de presin total.

En la Tabla J 1 del apndice 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 coeficientededifusin 555
Los dos coeficientes binarios de difusin se pueden corregir. para diferencias de presin y
temperatura utlilizando la ecuacin (24-41),
D A B condicin 1 /
I
DA, condicin 2
<\ 52\
Para 298 K y 2 atm, se tiene:
Las composiciones del oxgeno del nitrgeno 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 ecuacin (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 teora cintica

avanzada para explicar el movimiento molecular,las teoras de que se dispone
para explicar la estructura de los lquidos y sus caractersticas de transferen-
cia an son inadecuadas para permitir un tratamiento riguroso.La inspeccin
de los valores experimentales publicados en el Apndice 52,relativos a los
coeficientes de difusin de los lquidos, revelan que existen diversos rdenes
de magnitud menores que los coeficientes de difusin de los gases y que depen-
den de l a concentracin debida a los cambios de viscosidad con l a concentra-
cin as como a modificaciones en el grado en que la solucin sea ideal.
Ciertas molculas se difunden como molculas, en tanto que otras, que
se llaman electrlitos, 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 elctr~ica de la solucin indica que los iones deben difundirse en alguna
proporcin; por lo tanto, es posible hablar de un coeficiente de difusin de

10s electrlitos tales como el DjaC1. Sin embargo, si se encuentran presentes
varios iones, se debe tener en cuenta la rapidez de difusin de los cationes y
aniones individuales y los coeficientes de difusin molecular no tienen signifi-
cado. No es necesario decir que se necesitarn correlaciones separadas que se
utilizan en la prediccin de l a relacin entre las difusividades de la masa l-
quida, para los electrcilitos, as corno para 10s no electr6litos.
Se han postulado dos teoras, l a del hoyo de Eyring y la hidrodinmi-
ca, como posibles explicaciones al fenmeno de difusin de solutos no elec-
trolticos en soluciones de baja concentracin. En el concepto de Eyring el
lquido ideal se trata como un modelo de capas casi cristalinas interespacia-
das con hoyos. Entonces se describe el fenmeno de transferencia por medio
de un proceso que incluyeel salto de las molculas de soluto en los hoyos del
modelo de capas. Estos saltos se relacionan empricamente con la teora de
Eyring de rapidez de reaccin*. La teora hidrodinmica establece que el coefi-
ciente de difusin de los lquidos se relaciona con la movilidad de las molcu-
las de soluto; esto es, con l a velocidad neta de la molcula mientras sta se
encuentra bajo l a influencia de una fuerza motriz unitaria. Las leyes de l a hi-
drodinmica establecen relaciones entre l a fuerza y l a velocidad. Una de las
ecuaciones desarrolladas a partir de la teora hidrodinmica es l a ecuacin de
Stokes-Einstein:
donde D A B es l a difusividad de A en una solucin diluida de B,k es la conso-

nante de Boltzmann, T es la temperatura absoluta,7 es el radio de l a partcula
de soluto y pB es l a viscosidad del solvente. Esta ecuacin ha sido bastante
til en l a descripcin de la difusibn de las patculas coloidales de molculas
redondas y grandes, en un solvente que se comporta como un continuo enre-
lacin con l a especie en difusin.
Los resultados de ambas teorasse pueden reordenar en la forma general:
(24-44)
en la cual f(u) es una funcin del volumen molecular del soluto en difusin.
Se han obtenido correlaciones empricas en un intento de predecir el coefi-
ciente de d i f u s i h de los lquidos en funcin de las propiedades del soluto y
del solvente. Wilke y Changt propusieron l a siguiente correlacin, que es la
mejor de que se dispone para los no electrlitos en una solucin diluida
* S. Glasstone, K. J. Laider y H. Eyring, Theory o f Rate Processes, McGraw-Hill Book Company, Nue-
va York, 1941, Captulo 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 solucin, 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'ullicin, en cm3/g mol y
aBes el parmetro de "asociacin" correspondiente al solvente B.
Los volmenes mokculares a los puntos normalesdeebullicin, v6 ,
aparecen tabulados en la Tabla 24.3. Para otros compuestos, los volmenes
atmicos de cada uno de los elementos presentes se suman mediantelas frmu-
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 configuracin especfica de anillo. Se recomiendan las siguientes correc-
ciones:
para el anillo triple, como el xido de etileno restar 6
para el anillo cudruple 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
Hidrgeno H, 14.3 Oxido
ntrico NO 23.6
Oxgeno, O, 25.6 Oxido nitroso!NzO 36.4
Nitrgeno N, 31.2 N H,, 25.8
Aire 29.9 d;bo $ 5 18.9
Monxido decarbono Co 30.7
Bixido de carbono CO, 34.0
hidrgeno
Bromuro Br,
& 32.9
53.2
Sulfuro decarbonilocos 51.5 Cloruro Cb 48.4
Bioxido
sulfrico SO, 44.8 Yodo I, 71.5
LOS valores recomendados del parimetro de asociacin, aBpara unos Po-

COS solventes se proporcionan en seguida:
E.JEhlPL0 4
Calcidese el coeficiente de difusindel lquido, 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
atmico
Volumen
atmico
Volumen
Elemento
cm3
en /g Elemento
mol en cm3 /g mol
Bromuro 27.0 Oxgeno,
excepto en la 7.4
Carbono 14.8 forma anotadaabajo
Cloro 21.6 Oxgeno, en los steres 9.1
Hidrgeno 3.7 metlicos
Yodo 37.0 Oxgeno en los teres 9.9
Nitrgeno, doble 15.6 metlicos
enlace Oxgeno enteres
los 11.0
Nitrgeno en los 10.5 altos y otros steres
aminos primarios Oxgeno, en los 12.0
Nitrgeno 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 10C, la viscosidad de una solucin que contiene 0.05 mol de alcohol porcada li-
tro de agua es de 1.45 centipoises; los parmetros restantes quevan a usarse, son:
F283K
correspondiente al agua=2.6
para agua=l8
si se constituyen estos valores en la ecuacin (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 Apndice J.
Laecuacin(24-45) se recomiendasolamenteensoluciones diluidas

de solutonodisociativos.Para estas soluciones,predice valores cuyomar-
gen de error es de 13 % con relacin a los datos experimentales*. La ecua-
cin ( 2 4 - 2 5 ) n o es aplicable a la difusin de las molculas 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 molculas grandes de soluto.
Las propiedades de las soluciones conductoras de la electricidad se han
estudiado ya en forma intensiva durante ms de 75 aos. An as, las relacio-
nes conocidasentrelaconductividadelctricayelcoeficiente de difusin
del lquido son vlidas nicamente para soluciones de sales diluldas en agua.
El coeficiente de difusin de una sal equivalente en solucin dilulda, est da-
da por la ecuacinde Nernst:
(24-47)
donde DA, es el coeficiente de difusin basado en la concentracin 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
inicas lmite (concentracin cero) en (amp/cm2) (volt/cm) (g equivalente/
cm3 ) y es la constante de Faraday, de 96,500 coulombs/g equivalente. Esta
ecuacin se ha ampliado a iones polivalentes reemplazando la constante nu-
mrica2por (l/n' !-1/n-)donde n y n- son las valencias delcatiny del
anin, respectivamente
DIFUSIVIDAD DE LA MASA SOLIDA
En cualquier estudiodel movimiento molecular en el estado slimdo, la ex-

plicacin de la transferencia de masa se divide autornticamente en dos cam-
pos mayores de inters, la difusin de los gases o lquidos en los poros del
slidoylainterdifusin de los constituyentes slidos por medio del movi-
miento atmico. La primera clase de difusin tiene un papel preponderante
en la catlisis y es importante para el ingeniero qu.mico. Los metalurgistas
son los principales investigadores de la difusin de los tomos en los slidos.
La difusin en los poros se puede llevar a cabo por medio de tres o ms
mecanismos: difusin de Fick, difusin de Knudsen y difusin superficial. Si
los poros son grandes y el gas relativamente denso, la transferencia de masa
se llevar a cabo por medio de una difusin de Fick. Dentro del catalizador,
las trayactorias de difusin son como canales de forma irregular; por lo tanto,
el flujo es menor de lo que sera en poros uniformes de la misma longitud y

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

fase de difusin, tales como la temperatura y la presin, as como de las pro-
piedades del catalizador, tales como espacio fracciona1 vaco,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 relacin 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 tamao de los poros es del or-
den del de la trayectoria media libre dela molcula en difusibn. Se ha encon-
trado una relacin que describe la difusin de Knudsen y su coeficiente de
difusin, usando la teora cintica 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

partcula catalizadora y 111 el peso molecular del gas en difusin
La difusin superficial tiene lugar cuando las molculas que se han ab-
sorbido son transportadas a lo largo de la superficie como resultado de un
gradiente bidimensional de concentracin superificial. f Usualmente la difu-
sin superficial tiene un papel menor en la difusin en un slido catalizador a
menos que exista unagran cantidad de absorcin.
La introduccibn de un tomo extrao puede alterar las propiedades im-
portantes de un metal.
El edurecimiento del acero es un proceso industrial basado en la difu-
sin del carbono y otros elementos en el hierro. Se han producido semicon-
ductores por medio de la difusin de tomos "dopados" dentro de cristales.
Los cuatro mecanismos que se postulan para la difusin en los slidos son: la
difusin de vacancia, la intersticial, la de intersticialidad y la difusin por in-
tercambio directo de molculas 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 trmico a tempe-

raturas superiores al cero absoluto, existen sitios desocupados o vacantes en
las redes. Un tomo puede saltar de su posicin 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)Difusinlugar
de vacante (b) Difusin intersticial
O.... ...e.
0 .A.
0 . 0 . .
0 . 0 . .
0
(C) Difusin de intersticialidad

A.
....e
(d) Intercambio directo
e
Figura 24.3 Difusin en un cristal.
El tomo en difusin contina movindose a travs 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 distorsin de la red. Este
mecanismo se ha descrito matemticamente suponiendo un proceso de rapidez
unimolecular y aplicandoel concepto de Eyring de estado activado, tal como
se estudi en la teora de los hoyos correspondiente a la difusin de los l-
quidos. La ecuacin resultante es complicada y no se estudiar aqu. Rela-
ciona la difusividad con la relacin geomtrica entre las posiciones en la red,
la longitud de la trayectoria de salto, la frecuencia de vibracin del tomo
en la direccin del salto y la energa de activacin asociada al salto. La difu-
si& de vacunciu es el mecanismodominante de difusin en los metales y
aleaciomes cbicos centrados en las caras. El mecanismo se ha usado para ex-
plicar la difusin en materiales cbicos centrados en el cuerpo as como en
compuestos inicos 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 dilatacin o distorsin de la red; as pues, el mtodo matemtico
que incluye la teora de Eyring de rapidez unimolecular tambin se usa pa-
ra definir la difusividad intersticial. Este coeficiente de difusin se define en
funcin de las mismas variables que en el caso de la. difusin de vacancia.
Cuando el tomo intersticial en difusin posee un tamao comparable al
de tomos de la red, tiene lugar un proceso diferente del de difusin 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 autodifusin en
los metales y aleaciones incluye el intercambio directo de ms tomos por
medio de una rotacin en forma de anillo. Los tres tomos giran como en un
tiovivo y despus de dos intercambios, el tomo A se corre hasta una posi-
cicjn nueva, A . Despus efecta una rotacin con los otros dos tomos, Ile-
gando finalmente a A . El mecanismo de anillo no se sabe que ocurra en
ningim metal o aleacin, pero se ha sugerido como un mecanismo que puede
explicaralgunasanomalasaparentesen los coeficientesdedifusin de los
metales con redes centradas en el cuerpo.
Existen referencias excelentes para un estudio ms detallado de la trans-
ferencia de masa de los tomos a travs 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 caractersticas dinmi-
cas del fluido que est fluyendo.
Se debe lhacer una distincin 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-
veccin forzada. Si el movimiento del fluido se debe a una diferencia de den-
sidades, que puede haber surgido como resultado de una concentracin o de
una diferencia de temperatura, el procesose llama conveccin libre o natural.
La ecuacin de rapidez correspondiente a l a transferencia convectiva de
masa, generalizada en forma anloga a la ley de Newton del enfriamiento,
ecuacin (15-1l ) , es:
NA kc ACA (24-52)
donde NA es l a transferencia de masa molar de la especie A , medida con re-

lacin a coordenadas especiales fijas; AcA es la diferencia entre l a concentra-
cin de l a superficie lmite y la concentracin media de l a corriente de fluido
de la especie A e n difusin 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.
Conclusin 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 direccin de
una concentracin decreciente. La ecuacin (24-52) define el coeficiente Kc
en funcin del flujo de masa y la diferencia de concentracin desde el princi-
pio hasta elfinalde la trayectoria de difusin. Por lo tanto, el coeficiente
incluye las caractersticas de las regiones de flujos laminar y turbulento del
fluido, en cualesquiera proporciones en las que se encuentren. En los captu-
los 28 y 29 se estudiarn los mtodos de determinacin de este coeficiente.
En general es una funcin de: la geometra del sistema, las propiedades del
fluido del flujo, y la diferencia de concentracin 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 partculas prximas a la frontera slida 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 travs de
las capas estancaday laminar de fluido. L a resistencia que controla la transfe-
renciaconvectiva de masaes,amenudo,elresultado de esta "pelcula"
de fluido y el coeficiente Kc se llama, de acuerdo con esto, coeficiente de la
pelcula.
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 tcnicas desarrolladas para la evalucin del
coeficiente de transferencia convectiva de calor s e puede aplicar a la transfe-
rencia convectiva de masa. En el Captulo 28 se har un estudio completo de
los coeficientes de transferencia convectiva demasa as como de su evalua-
cin.
~~
24.4 C O N C L U S I O N
En este captulo se han introducido los dos modos de transferencia de

masa, molecular y convectivo. Como la difusin de masa incluye una mezcla
mltiple 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 difusin o difusividad de la masa en sistemas gaseosos,
lquidos y slidos 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 lquido, GNL, de la Pennsula de Kenai, en

Alaska, a un lugar de la Baha Yakina, en Oregon, por medio de un
transportemartimo. La composicinmolardelGNLcomercial es:
metano Ch, 94.9%

etano C, H, 4 .O%
propano C, H, 0.6%
bixidodecarbono CO, 0.5%
Determinese:
(a) la fraccin 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 presin parcial del metano cuando la presin total del sistema es
de 1.O13X1O5 Pa.
(e) la fraccin de masa del propano en partes por milln.
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) fraccin molar de C, I & ,
(b) fraccin de peso de N,
(c) peso molecular medio de la mezcla de gas natural
(d) presin parcial de cada una de las componentes si la presin 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) fraccin molar de O,,
(b) fraccin volumtrica de O,,
(c) peso de la mezcla
25
ECUACIONES DIFERENCIALES DE
LA TRANSFERENCIA DIE MASA
En el captulo 9 se obtuvieron las ecuaciones diferenciales generales de

transferenciademomentoutilizando el conceptodevolumendiferencial
de control. En el captulo 16 se generaron en forma ;anlogalas 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, tambin, la ecuacin de continuidad de una especie dada haciendo
un balance de la masa en un volumen diferencial de control.
Se obtendrn ecuaciones diferenciales adiciona.les al incluir en la ecua-
cin de continuidad, las relaciones correspondientes al flujo de masa encon-
tradas en el captulo 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 travsdelcualest

fluyendo una mezcla, incluyendo la componente A , tal como se aprecia en
la figura 25.1. La expresin de volumen de control clue corresponde a la con-
servacin de la masa, es:
lo cual puede decirse con palabras, de la siguiente manera:
Rapideznetadeflujoderapideznetadeacumulacin
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 relacin debe

incluir un trmino que corresponda a la produccihn o desaparicin de A por
medio de una reaccin qumica que ocurra dentro del volumen. La relacin
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
[acumulacin
volumen
de
control
volumen
de I[de A dentro del
duccin
control
- qumica
(25-1)
1
Se evaluarn los trminos correspondientes a la especieA y se explicarn
sus significados mis adelante.
I,a rapideznetadeflujo de masadel volumendecontrol se puede
calcular tomando en cuenta la masa transferida a travs de la superficie de
control. Por ejemplo, la masa de A transferida a travs 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 , sera
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 direccin 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 ecuacin diferencial dt?transferencia de masa 573
Si A se produce dentro del volumen de control por medio de una reaccin

qumica a una rapidez r A , donde r A tiene las unidatdes: (masa producida de
A)/(volumen)(tiempo), la rapidez de produccin de .4 es:
rA AX Ay AZ
Este trmino, correspondiente a la produccin, es anlogo al trmino que

corresponde a la generacin de energa, que apareceen la ecuacin diferencial
de la transferencia de energa, tal comose estudi en el captulo 16.
Si se sustituye cada uno de los trminos dela ecuacin (25-l), se obtiene:
(25.2)
Si se evala esta expresin en el lmite, cuando Ax,, Ay y AZ tienden a cero,

se obtiene:
(25-4)
La ecuacin (25-4) es la ecuacin 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 ecuacin (25-4) se puede escribir de la siguiente forma:
V.n,+--rA=O
aPA (25-5)
at
Se puede obtener una ecuacin 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 producirB dentrodel volumen d t control

por medio de una reaccin qumica. 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 conservacin de la masa y sustituyendo estas relaciones

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

(25-9) es idntica a la ecuacin de continuidad (9-2), correspondiente a un
fluido homogneo.
La ecuacin de continuidad de la mezcla y de unaespecie dada, se pueden
escribir en funcin de la derivada sustancial. La ecuacin de continuidad de
la mezcla se puede reordenar y escribir de la siguiente manera, tal como se
demostr en el captulo 9:
2 + p v .v = o (9-5)
Dt
A travs de pasos matemticos semejantes se puede obtener la ecuacin de

continuidad dela especie A en funcin dela derivada sustancial.La ecuacines:
r, =o (25-10)
El mismo mtodo se puede emplear en funcin de las unidades molares.

Si RA representa la rapidez de produccin molar de B por unidad de volumen,
las ecuaciones molares equivalentes, son:
Componente A :
(25-11)
Formas especiales de la ecuacin diferencial 575
Componente B:
(25-12)
Mezcla :
(25-13 )
Sin embargo slo cuando la estequiometra de la reaccin

es:
lo cual establece que se produce una molcula de B por cada mol de A que
desaparece y se puede decir que R, = -RB. En general, la ecuacin 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 ecuacin de contilnuidad vlidas en las situa-

ciones comunes, son las que aparecen a continuacin. Para poder usar las
ecuaciones en la evaluacin de los perfiles de concentracin, se reemplazan
los flujos, nA y N, por las ecuaciones apropiarlas obtenidas en el captulo 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 ecuacin (24-22) en la (25-5) se obtiene:
(25-15)
y al sustituir la ecuacin (24-21) en la (25-11) se obtiene:
Cualquiera de las dos ecuaciones: (25-15) o (25- 16) se pueden usar en la

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

constante, la ecuacin (25-15) se trhnsforma en:
Si se divide cada uno de los trminos entre el peso molecular de A y se reor-

dena la expresin, se obtiene:
v VCA + 3CA
~
at
+ RA
= DABV2cA (25-17)
I ii)Si no hay trmino de produccin, RA = O y si

coeficientesde
los
J
densidad y difusin se suponen constantes, la ecuacin (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 ecuacin (25-18), se obtiene:
DCA
"
- D,,v'c* (25-19)
D
t
Formasespecialesde la ecuacin diferencial 577
que es una ecuacin anloga a la (16-14) de transferencia de calor,
(16-14)
DT
--=aV2T
Dt
donde 01 es la difusividad trmica. L a semejanza entre estas dos ecuaciones

sirve de base para las analogas sealadas entre las transferencias de calor y
de masa.
iii) En una situacin en la cual no exista movimiento del fluido, v = O,

ni trmino de produccin, R, = O , ni variacin alguna en la difusividad o en
la densidad, la ecuacin (25-18) se reduce a:
(25-20)
L a ecuacin (25-20) se denomina segunda "ley" de Fick de la difusin. La

suposicin de que no hay movimiento del fluido restringe su aplicabilidad a
la difusin en slidos o lquidos estacionarios y a sistemas binarios de gases
o lquidos, donde N, es igual en magnitud pero acta en direccin opuesta
a NB ; esto es, en el caso de la contradifusin equimolar. La ecuacin (25-20)
es anloga a la segunda "ley" de Fourier de la conduccin de calor:
aT
-= aV2T (16-18)
at
iv) Las ecuaciones (25-17), (25-18)y (25-20) se pueden simplificar an

ms cuando el proceso que se va a definires un proceso en estado permanente;
esto es: =O. La ecuacincorrespondiente
una
a densidad constantey
una difusin tambin constante, es:
v VCA= DABV2cA+ RA (25-21)

Cuando la densidad es constante, lo mismo que la difusividad y no hay pro-
duccin qumica, RA = O , se obtiene:
v .VCA = DABV 2 C A (25-22)
Si adems, v = O , la ecuacin se reduce a
V2CA= o (25-23)
L a ecuacin (25-23) es la ecuacin de Laplace en funcin de la concentracin
molar.
Todas las ecuaciones, desde la (25-15) hasta la (25-23) estn escritas en

forma vectorial, por lo cual, todas son vilidas en cualquier sistema ortogonal
de coordenadas. La transformacin de la ecuacin 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 difusin, en coordenadas rectangulares es:
(25-24)
y en coordenadas cilndricas
(25-25)
y en coordenadas esfricas es:
L a ecuacin diferencial general de transferencia de masa de la componente A ,

o la ecuacin de continuidad de A , en coordenadas rectangulares, es:
(25-27)
en coordenadas cilndricas es:
(25-28)
y en coordenadas esfricas
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
determinacin de las constantes d e integracin. 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 seccin 16.3 para la transferencia de energa. Tal vez sea de
utilidad al estudiante repasar, en esa seccin, lo referente a condiciones ini-
ciales y de frontera.
La condicin inicial en los procesos de transferencia de masa es la con-
centracin de la especie en difusin al principio del intervalo de tiempo bajo
estudio, expresada ya sea en unidades molares o en unidades de concentracin
de masa. La concentracin 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 expresin ms complicada, si la ,distribucin de concentra-
ciones es una funcin de las variables al iniciar la medicin del tiempo.
Las condicionesdefrontera que generalmente se encuentran, son las
siguientes:
1. Se puede especificar la concentracin de u:na superficie. Esta concen-
tracin puede estar en funcin de la concentracin molar, cA = cA1de frac-
ciones molares, yA = yAl , en gases o x, = x A 1 ,en lquidos y slidos; de
concentracin de masa, p A = p A l o de fraccin de masa, wA = w A 1 . Cuando
el sistema es un gas, la concentracin se relaciona con la presin parcial por
medio de la ley de Dalton;asimismo, la concentracin puede estaren funcin
de la presin parcial, pA =p = y A l P. En el caso especfico de la difusin de
un lquido a una fase gaseosa, la condicin de frontera en la superficie del l-
quido est definida para una solucin lquida ideal, por medio de la ley de
Raoult, de la siguiente manerap, = x A PA , donde x, es la fraccin molar en
el lquido y PA es la presin del vapor de la especie A evaluada a la tempera-
tura del lquido.
2. El flujo de masa en la superficie puede estar especificado; por ejemplo:
.
j, = j, o NA = NA1 Los casos de inters en inge:niera 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 . Larapidezdeunareaccinqumicapuedeestar especificada, por

ejemplo, si A desaparece en la frontera debido a una reaccin qumica de
primer orden, se puede escribir N A I = k 1 C ~ 1 ,donde kl es la rapidez constante
en una reaccin de primer orden. Cuando laespecie en difusin desaparece en
la frontera debidoa una reaccin instantnea, la concentracin de esa especie
se supone, a menudo,igual a cero.
4. Cuando un fluido fluye en la fase para la cual est escrita la ecuacin
diferencial de masa, la especie puede perderse de la fase de inters por medio de
la transferencia convectiva de masa. El flujo de masa en la frontera, se define,
entonces, as:
donde CA, es la concentracin de la corriente de fluido adyacente a l a super-

ficie y kc esel coeficiente detransferenciade masa deconveccin que se
defini en la seccin 24.3.
Se aplicarn estas condiciones iniciales y de frontera en la solucin de
problemas de difusin molecular en los captulos 26 y 27. En seguida se pre-
sentar un ejemplo de reduccin en la ecuacin diferencial general de trans-
ferencia de masa de manera que solamente quedenlos trminos relevantes.
EJEMPLO 1
En una barra cilndrica de combustible nuclear que contiene material fisionable, la

rapidez de produccin de neutrones es proporcional a la produccin de los mismos. Use
una de las ecuaciones diferenciales generales de transferencia de masa para escribir la ecua-
cin 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 ecuacin (25-1 1) estable lo siguiente:
ac,
V.N,+--R,=O
at
Como la rapidez de produccin es proporcional a la concentracin neutrnica R , = kc,.

En los casos de difusin en slidos, en los que la contribucin del movimiento global es
cero,
N, = -DABVc,
Al sustituir estas relaciones en la ecuacin ( 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~
cin no vara con el ngulo, 8, ?cA/d8z = O ; la ecuacin se reduce a:
La nica condicin de frontera que resulta obvia es:
la cual requiere que la concentracin de la especie en difusin sea finita en el centro del
cilindro.
Problemas 581
L a ecuacin diferencial general de transferencia de masa se ha obtenido

en este captulo para la componente en difusin de una mezcla. Se presentaron
las formas de la ecuacin general, vlidas en situaciones especficas. Tambin
se listaron las condiciones de frontera encontradas con mayor frecuenciaen los
procesos de transferencia de masa.
PROBLEMAS
25.1 Obtengalaecuacin (25-11) correspondientealacomponente A en

funcin de lasunidades molares. empezartdo con la expresin de vo-
lumen de control correspondiente a la conservacin de la masa.
25.2 Obtenga la ecuacin de continuidad de la. componente A , utilizando
un volumen de control cilndrico e infinitamentegrande y coordenadas
cilndricas.
25.3 Transforme la ecuacin (25-27) de coordenadas rectangulares:
a una ecuacin equivalenteen coordenadas cilndricas.

"$5.4 Una celda de difusin de Arnold es un apara.to sencillo usado para medir
coeficientes de difusin de gases. Un estanque de lquido se mantiene
en el fondo de un tubo de dimetro pequeo. Un gas, insoluble en el
lquido, fluye a travs de la boda del tubo extrayendo los vapores de A
que se difunden a travsdel espacio de gas que se encuentra sobre el
estanque de lquido. En condiciones isotrmicas e isobricas, la evapo-
racin de A es un proceso en estado permanente.
Reduzca la ecuacin diferencial general de transferencia de masa para
escribir la ecuacin diferencial especfica que describa este proceso de
transferencia de masa.
iCul sera la forma de la ley de Fick, correspondiente a la especie A
que se substituira en la ecuacin diferenc:ial? Proporcione dos condi-
ciones de frontera que podran ser utilizadas para resolver la ecuacin
diferencial resultante.
-Flow of gas B
Lquido A puro
28
TRANSFERENCIA CONVECTIVA
DE MASA
La transferencia de masa por conveccin es la que se realiza entre una

superficie lmite y un fluido en movimiento o entre dos fluidos en movimiento,
relativamente no miscibles. La ecuacin 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 direccin en la que dismi-
nuye la concentracin. Esta sencilla ecuacin es la relacin que define a kc,
que es el coeficiente de transferencia conuectiua de masa. Es anlogo ala
ecuacin 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 captulo 15, no fue una
empresa fcil.Amboscoeficientes de transporte estn relacionados con las
propiedades del fluido, con sus caractersticas dinkmicas y con la geometra
del sistema especfico 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 captulo 19 se pueda aplicar al coeficiente de transferen-
cia de masa y esto es, precisamente, lo que se har en los anlisis siguientes.
Se har un uso considerable de las secuelas y conceptos de los captulos 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 conveccin. 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 fsico de paquetes de materia alo largo de las lneas 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 distincin entre los flujos
laminar y turbulento ser de vital importancia en cualquier situacin convec-
tiva.
La capa lmite hidrodinmica que se analiz en el captulo 12, tiene un
papel muy importante en la transferencia convectiva de masa. Tambin se
definir y analizar una capa lmite de concentracin que resultar vital para
el anlisis del proceso de transferencia convectiva de masa. Esta capaes similar,
pero no necesariamente del mismo espesor que la capa trmica lmite que se
estudi en el captulo 19.
Cuando la transferencia de masa incluye a un soluto que se disuelve con
rapidez constante desde una superficie slida y despus se difunde a un fluido
en movimiento, el coeficiente de transferenciaconvectiva de masa se define as:
En esta ecuacin, 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-
posicin del soluto en el fluido, en la intercara, c ~ ,es~la, composicin del
fluido en equilibrio con el slido a la temperatura y presin del sistema. La
cantidad cA , representa la composicinen algn punto dentro dela fase fluida.
Cuando se define la capa lmite de la concentracin, se puede escoger a cA
como la concentracin de la componente A en el borde de la capa lmite y
sele puede llamar cA-. Si el flujo tuviera lugar en un conducto cerrado, la
composicin, c A , podra serla concentracin global o concentracin de
nzezcla homognea. La composicin de mezcla homognea es la concentra-
cin que se medira si se recogiera y se mezclara perfectamente el flujo en
un plano, esto es,se tendra una composicin promedio del flujo global.
Existen cuatro mtodos de evaluacin de los coeficientes de transferencia
convectiva de masa que se estudiarn en este captulo, que son:
1. anlisis dimensional acoplado con experimentos

2. anlisis exacto .de capa lmite
3. anlisis aproximado de capa lmite
4. analoga entre la transferencia de momento, energa
y masa.
En las secciones siguientes se estudiar cada uno de estos mtodos.

Parmetros 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 parmetros adimensionales para relacionarlos da-

tos relativos a la transferencia convectiva. En la transferencia de momento,
nos hemos topadoya con los nmeros de Reynolds, y de Euler. En la relacin
de datos de transferenciaconvectiva de calor, son importantes los nmeros de
Prandtl y Nusselt. Algunos de los mismos parmetrosjunto con ciertas razones
adimensionales recin definidas, resultarn de utilidad en la correlacin de
los datos correspondientes a la transferencia de masa. En esta seccin se es-
tudiar la interpretacin fsicad 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 trmica =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 razn de cualesquiera dos de ellas debe ser
adimensional. La razn de la difusividad molecular de momento a la difusi-
vidad molecular de la masa se denomina nmero de Schmidt:
difusividad de momento
-Sc="=-
- v El. (28-2)
dqusividad de masa DAB PDAB
El nmero de Schmidt tiene una importancia en la transferencia convectiva

de masa, anloga a la del nmero de Prandtl, enla transferencia convectiva de
calor. La razn de la difusividad molecular del calor a la de la masa, se deno-
mina nmero de Lewis:
difusividad trmica = = k (28-3)

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

vectiva simultnea de masa y energa. Los nmero,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 endifusin.
660 Transferencia convectiva de masa
Pensemos ahora en la transferencia demasa del soluto A desde un slido

hasta un fluido que pasa alrededor de la superficie del slido. En la figura
28.1 aparece el perfil de la concentracin.
Figura 28.1 Perfiles de la velocidad y la Concentracin correspondientes a un

fluido que circula alrededor de una superficie slida.
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 difusin mo-

lecular, la transferencia de masa tambin se puede escribir por medio de:
cuando la concentracin de frontera,CA,,,es constante, esta ecuacin se sim-

plifica, reducindose a:
(28-5)
Las ecuacines (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 relacin:
(28-6)
Si se multiplican ambos lados de la ecuacin (28-6) por una longitud signi-

ficativa, L , se obtiene la siguiente expresin adimensional:
(28-7)
Anlisis de la transferencia convectiva de masa 661
El lado derecho de la ecuacin (28-7) es la relacin entre el gradiente de la

concentracin en la superficie y un gradiente globad o de referencia de la con-
centracin. Por lo tanto, se puede considerar como la relacin 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 ecuacin (28-7) es
semejante a la de la ecuazin (19-5) del nmero cle Nusselt encontrado en la
transierencia convectiva de calor, la raznk,L/Dp,B se denomina nmero de
Nusselt de transfrencia de masa, NuAB.Tambin se le conoce como nmero
de Sherwood, Sh.
Estos tres pargmetros Sc, Nu,, o Sh y Le, aparecern en los anlisis 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 parmetros adimensionales

que resultan tiles en la correlacin de los datos experimentales. Hay dos
procesos importantes de transferencia de masa que se estudiarn: la trans-
ferencia de masa hacia una corriente que fluye en condiciones de conveccin
forzada y de transferencia demasa hacia unafase (quese mueve en condiciones
de conveccin 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 travs del conducto. La trans-
ferencia es el resultado de la fuerza impulsora de la concentracin, c ~ -cA.
, , ~
Las variables importantes, sus smbolos y sus represent?ciones dimensionales,
son los que aparecen a continuacin:
Vmiab le Smbolo Dimensiones
dimetro 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 trminos descriptivos de la geometra del

sistema, las propiedades del flujo y del fluido y la cantidad de principal in-
ters, k,.
Por mediodel mtodo de Buckinghamde agrupamientb de variables,

que se estudi en el captulo 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 ecuacin,se tiene para:
L: 0 = 2 ~ - 3 b + c + l
t: o=-u-1
Y
M : O=b
La solucin a estas ecuaciones, correspondiente a los tres exponentes incg-
nitos, produce el resultado:
U=-]
b=O
Y
c=l
as pues: 7 r 1 = k,D/DAB,que es el nmero de Nusselt de transferencia de masa
o bien, el nmero 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 nmero de Schmidt. Si se divide IT^ entre r3 se obtiene:
EL
Anlisis de la transferencia convectiva de masa 663
el nmero de Reynolds. El resultado del anlisis dimensional de la transfe-

rencia de masa en condiciones de conveccin forzada en un conducto circular
indica que la correlacin podra ser de la forma:
NuAB= f(Re, Sc) (28-8)
que es anloga a la correlacin 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 conveccin natural se crearn si existe cualquier va-

riacin en la densidad en una fase lquida o en una gaseosa. La variacin en
la densidad se puede deber a diferencias de temperatura o a diferencias de
concentracin, relativamente grandes.
En el caso en el que la conveccin natural incluye la transferencia demasa
desde una pared vertical plana hasta un fluido adyacente, las variables diferirn
de las utilizadas en el anlisis de la conveccin forzada. Las variables impor-
tantes, sus smbolos y sus representaciones dimensionales, aparecenen la lista
siguiente:
Variable Sz'mbo,!o Dimensiones
longitud caracterstica 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 sern, con DA.B,L, y p como variables principales:
Si se resuelven los tres grupos pi, se obtiene:

el nmero da Nusselt de transferencia de masa, o nmero de Sherwood:
el recproco del nmero de Schmidt y
El resultado del anlisis dimensional de la transferencia de masa por convec-

cin natural sugiere una interrelacin de la forma:
SC)
NUAB= ~ ( G ~ A B, (28-9)
Las relaciones que se obtuvieron mediante el anlisis dimensional, tanto
para la conveccin natural como para la forzada, sugieren que una correlacin
de los datos experimentales puede estaren funcin de tresvariables, en vez de
las seis originales. Esta reduccin en el nmero de variables ha ayudado a los
investigadores que han sugerido interrelaciones de esta forma para producir
muchas de las ecuaciones empricas que aparecen en el captulo 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 solucin exacta, correspondiente a la capa lmite

hidrodinmica para un flujo paralelo a una superficie plana. Esta solucin se
examin en la seccin 12.5. En la seccin 19.4 se hizo una extensin a la
solucin de Blasius para explicar la transferencia convectiva de calor. En una
forma absolutamente anloga se ampliar, ahora, la solucin de Blasius, para
incluir la transferencia convectiva de masa para lamisma geometra yel mismo
flujo laminar.
Las ecuaciones correspondientes a la capa lmite, estudiadas en la trans-
ferencia de momento en estado permanente, incluyeron la ecuacin bidimen-
sional e incompresible, de continuidad:
au, au
-+".?=o (12-1 l b )
ax ay
Anlisis de l a concentracin laminar de la capa lmite 665
y la ecuacin de movimiento en la direccin de x, para una presin y una v

constantes,
(12-1la)
La ecuacin que describe la transferencia de energa de un flujo isobrico en

estado permanente, incompresible y bidimensionatl, con una difusividad tr-
mica constante, correspondiente a la capa trmica lmite, es:
(19-15)
Hay una ecuacin que se aplica a la transferencia de masa dentro de una

capa lmite de concentracin cuando no hay produccin de la componente
en difusin y cuando la segunda derivada de cA con respecto a x, d2C,4/aX2,
es mucho ms pequea en magnitud que la segunda derivada de cA con res-
pecto a y . Esta ecuacin, escrita para un flujo permanente incompresible y
bidimensional, cuya difusividad de masa es constante, es:
(28-10)
La capa lmite de concentracin aparece en forma esquemtica, en la figura

28.2. Las siguientes son las condiciones de frontera correspondientes a las
tres capas lmite:
Orilla de
limite
la capa n
Figura 28.2 Capa lmite de concentracin de un flujo laminar alrededor de una

superficie plana.
o, como la velocidad en la direccin de x en la pared, vXTses igual a cero
T - T, T - T,
trmico -- -0 eny=O Y -- -1 eny=m
Tm - T s T m- T s
Yconcentracin
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 fenmenos de transferencia. Enel captulo
19 se modific la solucin de Blasius correspondiente a la ecuacin (12-1l a )
y se le aplic con xito en la explicacin de la transferencia convectiva de
calor cuando la razn del momento a la difusividad trmica, v / a = Pr = 1 se
utiliza el mismo tipo de solucin paradescribir la transferencia convectiva de
masa cuando la razn: momento-difusividad de masa, v / D = Sc = 1. Usando
la nomenclatura definida en el captulo 12:
(28-12)
al aplicar el resultado de Blasius ala capa lmite de

la concentracin se obtiene:
La ecuacin (28-13) se puede reordenar para obtener una expresin corres-

pondiente al gradiente de la concentracin en la superficie,
(28-14)
Es importante recordar que la solucin de Blasius correspondiente a la ecua-

cin (12-lla) no inclua una velocidad en la superficie en la direccibn de x ,
por lo tanto, la ecuacin (28-14) hacc la suposicin, muy importante, de que
La rapidez con la cual entra o sale la masa de la capa lmite enla superficie, es
tan pequeia que no altera el perfil de la velocidad predicho por la solucin
de Blasius.
Cuando la velocidad en la direccin de y en la superficie, UW es igual a
cero, el trmino correspondiente a la contribucin global, en la ecuacin de
Fick, el flujo de masa en la direccin de y tambin es igual a cero. La trans-
ferencia de masa de la superficie pIana de la capa laminarlmite, est descrita
por la ecuacin:
(28-15)
Anlisis de la concentracin de la capa lmite 667
AI sustituir la ecuacin (28-14) en la (28-15),se okltiene:
(28-16)
El flujo de masa de la componente en difusin se defini en funcin 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-
tenindose:
DAB
kc =-[0.332
X

-= NuAB = 0.332 Re,12 (28-17)
DAB
La ecuacin (28-17) est restringida a sistemas cuycl nmero de Schmidt, Sc,

sea igual a uno y con una rapidez de transferencia de masa entre la placa plana
y la capa lmite, muy baja,
En la figura 28.3 puede observarse una representacin grfica de la so-
lucin a la ecuacin de capa lmite de concentracin, ecuacin (28-lo), ob-
tenida por Heartnett y Eckert. En esta figura se muestran las curvas que
representan los valores positivos y negativos del parmetro lmite de la super-
ficie, (uY,JudRex)?Los valores positivos son vlidos cuando la transferencia
de masa de la placa plana ocurre hacia la capa lmite y los valores negativos
describen la transferencia de masa del fluido la placa.
a Cuando el parmetro de
la superficie de frontera tiende a cero, la rapidez de transferencia de masa
disminuye hasta que se considera que no tiene ningn efecto en el perfil de
velocidad. La pendiente de la recta cero, calculada eny = O es 0.332, tal como
lo predice la ecuacin (28-13).
En la mayora de las operaciones fsicas que se relacionan con la trans-
ferencia de masa, el parmetro de la superficie lmite es despreciable y se uti-
liza el tipo de solucin de Blasius, de baja transferencia de masa, para definir
la transferencia hacia la capa laminar lmite. La vaporizacin de un material
voltil hacia una corriente gaseosa que fluye a baja presin, es un caso en el
que no se puede hacer la suposicin 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 concentracin correspondientes a la transferencia de

masa en una capa laminar lmite sobre una placa plana.
Para un fluio cuyo nmero 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 mtodo anlogoalasolucin de Pohlhausen,correspondiente a la
transferencia de calor convectivo, para la transferencia convectiva de masa.
La capa lmite de concentracin est relacionada con la capa hidrodinmica
lmite por medio dela ecuacin:
(28-18)
donde 6 es el espesor de la capa hidrodinmica lmite y 6, es el grosor de la

capa lmite de concentracin, por lo cual, el trmino de Blasius, q, se debe
multiplicar por S C ~ .En la figura 28.14 aparece una grfica de concentracin
adimensional contraqSc correspondiente a uy.s =O. La variacin en la con-
centracin, expresada en esta forma conduce a una expresin correspondiente
Anlisis de la concentracin laminar de la capa lmite 669
al coeficiente de transferencia convectiva de masa, semejante a la ecuacin

(28-17). En y = O, el gradiente de la concentracin es:
(28-19)
el cual, al emplearse con la ecuacin (28-15)da origen a la ecuacin:
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 integracin. 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 nmero local deNusselt a una distanciax , en. elsentido dela corriente,
se relaciona con el nmero medio de Nusselt correspondiente a la placa, por
medio de la expresin:
NUL,AB =2NUx,AB(x=L (28-22)

que es anloga 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 anlisis, totalmente diferente, ha producido
*W. J. Christian y S. P. Kezos, A.1.Ch.E. J., 5 , 61 (1959).
resultados de la misma forma predicha en la seccin 28.3 mediante el anlisis

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

Si se recuerdan los perfiles adimensionales de concentracin 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 parmetro positivo de la superficie lmite ( v /v,)(Re)*, aumenta.
y:
Como la magnitud del coeficiente de transferencia esta directamente relacio-
con
nado ecuacin:
lalapendiente
medio
depor .
(28-23)
la disminucin de Ia pendiente indica que los sistemas cuyos valores del par-
metro de la superficie lmite son altos, tienen coeficientes ms bajos de trans-
ferencia de masa.
Cuando la energa y la masa se transfieren a travs de la capa laminar
lmite, los perfiles adimensionales de la figura 28.3 tambin pueden representar
los perfiles adimensionales de temperatura cuando los nmeros de Prandtl y
Schmidt son iguales a la unidad. En el prrafo anterior, se seal que el coe-
ficiente de transferencia de masa disminuye en magnitud cuando se transfiere
masa de la superficie a la capa lmite, por lo tanto, es de esperarse que el coe-
ficientedetransferenciadecalordisminuya al transferirse masa a la capa
lmite. Esto puede realizarse haciendo pasar un fluido a travs de una placa
porosa hacia afuera, llegando a la capa lmite, o sublimando el material mismo
de la placa. Estos procesos simultneos de transferencia de calor y masa, a los
cuales se llama a menudo enfriamiento por transpiracin y ablacin, respec-
tivamente, se utilizanparaayudarareducir los enormesefectosdelcalor
durante el regreso de un proyectil ala atmsfera terrestre.
EJEMPLO 1
~1 coeficiente de transferencia de masa correspondiente a una capa lmite turbulenta

formada sobre una placa plana se ha correlacionado, en funcin del nmero 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.

transicin de flujo laminar a turbulento ocurre en Re, = 3 X 10.
a ) Encuntrese una expresin que corresponda al coeficiente medio de transferencia
de masa en una placa plana de longitud L , .
Por definicin:
(28-25)
Anlisis de la concentracinde la capa lmite 671
donde L, es la distancia medida del borde de ataqueal punto de transicin. kc,,,, est defi- .
nida por la ecuacin (28-20):
DAB
k,,,, = 0.332---(Re,)/2(Sc)/3
X
kc,turbest
definida
por ecuacin
la (28-24): +
Das
kc.turt,= 0.0292-(Re,)(S~)/~
X
AI sustituir estas dos ecuaciones en la expresin correspondiente al coeficiente medio de

transferencia de masa, se obtiene:
donde L, es la distancia del borde de ataqueal punto de traasicin, 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,~)45 -(Re,)/]
kc = (28-26)
L
b ) Un vaso de alcohol etlico 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 travs
del banco de 1 m de ancho. El aire se mantuvo a 289 K y 1 atm. (1.013 x lo5 Pa). La
presin del vapor de alcohol etlico a 289 K es de 4,000 Pa. Determnese la cantidad de
alcohol que se evapora de un rea de un metro cuadrado cadar 60 seg.
x
A 289 K, la viscosidad cinemtica 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 nmero 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 transicin
donde la capa lmite cambia de flujo laminar a turbulento. El :punto de transicin se puede
evaluar a partir del nmero de transicin 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 ecuacin (28-26)

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

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

Se puede calcular la concentracin del alcohol etlico en el vapor inmediatamente encima
de la superficie del lquido, por medio de la relacin:
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-

racin es otra que no sea una placa plana, existen actualmente pocas soluciones
exactas correspondientes a la transferencia en una capa limitada. El mtodo
aproximado, encontrado por von Karmn para describir la capa lmite hidro-
dinmica, se puede utilizar en el an6lisis de la capa lmite de concentracin.
En los captulos 1 2 y 19 se estudi el uso de este mtodo.
Analicese, ahora, el volumen de control localizado en l a capa lmite de
concentracin, tal como aparece enla figura 28.5. Este volumen representado
por las lneas punteadas, tiene un anchoigual a Ax, una altura igual al espesor
!V
Figura 28.5 Volumen de control de la capa lmite de concentracin.

Anlisis aproximado de la capa lmite de concentracin 673
de la capa lmite de concentracin, S, y profundidad unitaria. Si se hace un

balance de la masa molar en estado permanente sobre el volumen de control,
se tendr la relacin:
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 trminos de rapidez molar, se puede reescribir la ecuacin (28-27) en la
forma:
(28-28)
AI redordenar y dividir cada uno de los trminos entre Ax y calcular el resul-
tado en el lmite, cuando Ax tiende a cero, se obtiene:
d
& j,,,(CA - cA,m) ux dy = k c ( C A S - CA,m) (28-29)
La ecuacin (28-29) es anloga a las ecuaciones (12-37) y (19-30). Para

poder resolver la ecuacin (28-29) deben conocerse los perfiles de la velocidad
y la concentracin 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 concentracin supuesto debe satisfacerlas condiciones correspon-

dientes de frontera en funcin 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 ecuacin integral de von Krmn, que es la (28-29), para obtener una
solucinaproximada. Los resultados se puedencompararconlasolucin
exacta, ecuacin (28-20) y verificar la exactitud de las suposiciones hechas
acerca de la velocidad y de los perfiles de concentracin. Como primera aproxi-
macin, vamos a tomar una expresin en serie de potencias, para la variacin
de la concentracin con y.
2 3
~A-~A,~=a+by+cy+dy
Al aplicar las condiciones de fronterase obtendr la expresin:
(28-34)
Si se supone que el perfil de la velocidad tiene la misma formaen serie de po-

tencias, entonces la expresin resultante, tal como se vio en el captulo 1 2 es:
( 12-40)
Cuando se sustituyen las ecuaciones (28-34)y (12-40) en la expresin integral

(28-29) y se resuelve, se obtiene:
Nu,,AB = 0.36Re,2S~ (28-35)
resultadoque se aproximabastantealasolucinexactaexpresadaenla
ecuacin (28-17).
Aunque esta no es la expresin correcta,es suficientemente aproximada
para indicar que se puede usar el mtodo integral con bastante confianza en
Analoga de transferencia de la masa, energa y momento 675
otras situaciones en las que no se conozca una solucin exacta. La exactitud

del mtodo depende totalmente de la habilidad para suponer buenos perfiles
de velocidad y concentracin.
La ecuacin integral de von Krmn, ecuaciijn (28-29), tambin se ha
utilizado para obtener una solucin aproximada de la capa lmite 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 concentracin, igual a:

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

= 0.0292 (28-24)
Se anima al lector a llevar a cabo la obtencin de la ecuacin (28-24) que se
solicita en uno de los problemas que aparecen al final del captulo.
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 anlisis 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, energa y rnasa y en las condiciones
de frontera, cuandose expresaron los gradientes de tr,ansferencia en funcin de
variables adimensionales. Estas semejanzas nos han permitido predecir solu-
ciones correspondientes a los procesos semejantes de transferencia. En esta
seccin se estudiarn diversas analogas de entre lots fenmenos de transfe-
rencia, que se han propuesto debido a la semejanza de sus mecanismos. Las
analogas resultan de utilidad para comprender los fenmenos de transferencia
y tambin como un medio satisfactorio de predecir el comportamiento de
los sistemas para los cuales existen datos cuantitativos limitados disponibles.
La semejanza entre los fenmenos de transferencia y, por lo tanto, la
existencia de las analogas, requieren de la existencia de las siguientes cinco
condiciones dentro del sistema:
1. Que las propiedades fsicas sean constantes.
2. Que no haya produccin de masa o energa dentro del sistema. Esto,
desde luego, implica que no pueden ocurrir reaccione:: qumicas homogneas.
3. Que no exista emisin o absorcin de energa radiante.
4. Que no haya disipacin 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 anlogo

de la transferencia de energa y momento. Aunque esta analoga es de zpli-
cacin limitada, ha servido como catalizador para encontrar mejores analogas
y ha sido utilizada con xito en elanlisisde los complejos fenmenos de
capa lmite de la aerodinmica.
Reynolds postul que los mecanismos de transferencia y energa y mo-
mento eran idnticos.En el estudiorealizadoanteriomente,acerca de las
capas laminares lmite, se ha visto que esto se cumple cuando el nmero de
Prandtl, Pr es igual a la unidad. A partir del estudio que se hizo en la seccin
28.4, se puedeampliar el postulado de Reynolds para incluir el mecanismo
correspondiente a la transferencia de masa cuando el nmero de Schmith, Sc,
tambin 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 concentracin
dentro de las capas lmite estn relacionados por medio de:
a u,
"i
a y cc~A ,- c~A-, sc A , ~ ) l y = " = ~ ( ~ ) l y = " (28-36)
Recurdese que en la frontera prxima a la placa, donde y = O , se puede ex-
presar el flujo de masa en funcin de la difusividad de la masa o del coeficiente
de transferencia de masa, por medio de la ecuacin
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 expresin
que relacione el coeficiente de transferencia de masa con elgradientede la
velocidad en la superficie:
(28-38)
El coeficiente de friccinsuperficial se relacion, en el captulo 1 2 ,

con este mismo gradiente de velocidad, mediante la expresin:
( 12-2)
Si se hace uso de esta definicin se puede reordenar la ecuacin (28-38) para

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

La ecuacin (28-39) es anloga a la analoga de Reynolds de transferencia de

energa para los sistemas cuyo nmero de Prandtl es igual a uno. En el cap-
tulo 19 se analiz esta analoga, la cual puede expresarse as:
(19-36)
Los datos experimentales correspondientes a la transferencia de masa hacia

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

principal, es turbulento en lugar de ser laminar. .Aunque muchos investiga-
dores han hecho grandes aportaciones a la comprensin del flujo turbulento,
hasta el momento nadie ha tenido xito en la prediccin de los coeficientes
de transferenciaconvectiva o factores defriccin por medio del anlisis directo,
lo cual no resulta demasiado sorprendente cuando se recuerda, en relacin
con el estudio previo acerca del flujo turbulento, seccin 13.1, que el flujo
en cualquier punto est sujeto a fluctuaciones irregulares en direccin y ve-
locidad. De acuerdo con esto, cualquiera de las partculas 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, "difusin 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 ncleo turbulento.
En un esfuerzo para caracterizar este tipo de movimiento, Prandtl pro-
puso la hiptesis dela longitud de mezclado, tal como se estudi enel captulo
13. En esta hiptesis, cualquier fluctuacin de la velocidad, u%' se debe al
movimiento en la direccin 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 fluctuacinde la velocidad se relacionacon
el gradiente de la velocidad media, por medio de la ecuacin
(13-10)
El esfuerzo cortante total se defini por medio dela expresin:
( 13-8)
Al sustituir la ecuacin (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
tambin se debe a la presencia de las fluctuaciones o remolinos. En la figura
28.6, la curva representa una porcin del perfil turbulentolade concentracin
con un flujo medio en la direccin de x. La rapidez instantnea de transfe-
rencia de la componenteA en la direccin de y es:
I '
NA.y = CA uy (28-42)
donde CA = FA + c A ' , o sea la suma de
las fluctuaciones instantneasel promedio
temporal, de la concentracin de la componente A . De nuevo puede utilizarse
el concepto dela longitud
Figura 28.6 Porcin de la curva del perfil de la concentracin turbulenta, mos-

trando la longitud del mezclado, de Prandtl.
de mezclado, para definir la fluctuacin de la concentracin por medio de la

relacin:
(28-43)
AI insertar la ecuacin (28-43) en la (28-42), ST obtiene una expresin que

corresponde a la transferencia turbulenta de masa por medio de remolinos.
La transferencia total de masa, normallaa direccin 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 captulo 19,
una expresin correspondiente a la transferencia convectiva de calor:
(19-49)
donde CY es la difusividad trmica molecular y e,,, , l a difusividad trmica 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-
talticos heterogneos, los calefactores, los secadores, etc. Como resultadode
la difusin de remolino, la transferencia es rpida en el ncleo turbulento re-
duciendo la composicin de cualquier gradiente. Al acercarse a la pared, la
turbulencia va disminuyendo hasta que, enla vecindad inmediata de la super-
ficie del slido, desaparece esencialmente y la transferencia se realiza casi
completamente por difusin molecular. La mayor parte de la resistencia a la
transferencia ocurre en la capa lmite prxima a. lasuperficieen lacual el
gradiente de la composicin es ms pronunciado.
LAS ANALOGIAS DE PRANDT Y VON KARMAN

En el captulo 19 se obtuvo la analoga de Prandtl correspondiente a l a
transferencia de calorymomento,cuando se estudi el efecto del ncleo
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
analoga 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 ecuacin (28-41) se puede in-
tegrar sobre el espesor de la subcapa, resultando:
6
(28-45)
La ecuacin (28-44) tambin 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 analoga de Reynolds, k c / v m= cf/2

= 7,/pUm2, se puede utilizar en el ncleo
turbulento, de y = .$ a y en las condiciones globales. El flujo de masa del
ncleo 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 tambin:
en la ecuacin (28-49),se puede simplificar la relacin hasta quedar as:
1 - 2
2[vcD+ v x p - 1)1
-"
kc cpm
o, escrita de manera ligeramente distinta:
(28-50)
Ntese que la ecuacin (28-50) simplifica la analoga de Reynolds con la si-

guiente restriccin: S c = l. En el captulo 13 se defini la subcapa laminar
mediante la expresin: 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 ecuacin (28-50) se obtiene una

analoga correspondiente a la transferencia convectiva de masa, semejante a
la analoga 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 ecuaciln (28-52) por vooL/D~*,

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

transferencia de energa-momento, ecuacin ( 19-57).
Von Krmn amplila analoga de Prandtl, tomandoen cuenta la llamada
capa amortiguadora adems de l a subcapa laminar y del ncleo turbulento.
Esto lo condujo a la formacinde la analoga de uon Kcrmn:
(Cf/2jRePr (19-58)
Nu =
1 + 5JCf/2{Pr- 1 + In [( 1 + 5Pr)/6]}
que corresponde a la transferencia de energa y momento. El anlisis de von

Krmn 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 analogas 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 trminos, que
sirve de correccin a la analoga simple de Reynolds.
ANALOGIA DE CHILTON-COLBURN
Chilton y Colburn descubrieron que la mejor correlacin entre los datos

experimentales obtenidos en l a transferencia de masa, est dada por la ex-
presin:
(28-56)
Esta analoga es vlida para hquidos y gases dentro delos valores 0.6<Sc<2500.
El trmino jD se denomina factor j de transferencia de masa y es anloga al
factor j de transferencia de calor, definido por la ecuacin (19-39). Aunque
la ecuacin (28-56) es una ecuacibn emprica basada en los datos recogidos
en los regimenes de flujo laminar y turbulento, se puede demostrar que sa-
tisface la solucin exacta correspondiente a un flujo laminar sobre una placa
plana,
=0.332Re,Sc1 (28-20)
Si se dividen ambos lados de la ecuacin entreRe,Sc, se obtiene:
(28-57)
Esta ecuacin se reduce a la analoga de Chilton-Colburn cuando en la ex-

presin de arriba se sustituye la solucin de Blasius correspondiente a la capa
laminar Imite,
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 analoga completa de Chilton-Colburn es:
(28-59)
la cual relaciona los tres tipos de transferencia en una sola expresin. La ecua-
cin (28-59) es exacta para placas planas y satisfactoria cuando se emplea en
sistemas cuya geometra sea diferente, a condici6n de que no exista ningn
arrastre de forma. Ensistemas enlos que este ltimo tipo de arrastrese encuen-
tra presente, se ha encontrado que:
(28-61)
La ecuacin (28-61) relaciona la transferencia convectiva de calor con la de

masa; permite la evaluacin de un coeficiente desconocido de transferencia
a travs de la informacin obtenida para otros fenmenos de transferencia. Es
vlida para gases y lquidos, entre los valores: 0.6<Sc<2,500 y 0.6<Pr<100.
EJEMPLO 2
Fluye aire a presin atmosfkrica a travs de un termmetro cuyo bulbose ha cubierto
con una tela hmeda. Este clsico termmetro de bulbo himedo indica que una pequea
cantidad de lquido 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 termmetro es de 65" F.
A esta temperatura se calcularon las propiedades siguientes:
presin de vapor del agua O .305 6 psi

densidad
aire del ,-
0.076 lb, /ft2
calor latente de vaporizacin del agua : 1057 Btu/lb,
nmero de Prandtl 0.72
nmero de Schmidt 0.61
calor especfico, cp, del aire 0.24 Btu/lb," F
(Cul es la tempertura del aire seco?

La ecuacin (28-1) define el flujo molar del agua que se! evapora
NHzO = kc (cHKh -C H 2 0 , ~ )
La energa 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 vaporizacin a la temperatura de la superficie. Esta ecuacin

se puede resolver para la temperatura global:
Si se sustituye la ecuacin ( 2 8 - 1 ) en esta ecuacin, se obtiene:
kc
Tw = A M H ~ o ~ ( c -cHzo.m)
I - ~ ~ o ,+Ts
~
LOSfactores j de Chilton-Colburn nos dan una relacin que corresponde a la razn k,lh:
. .
IH = I D
h
-((pr)Z/3 kc
= _(SC)*/'
PmCp
Cuando se sustituye esta expresin en la ecuacin 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 aos se han estado usando los coeficientes de transfe-

rencia convectiva de masa en e1 diseo de equipo de transferencia de masa.
Sin embargo, en muchos casos han sido coeficientes empricos determinados
a partir de investigaciones experimentales. Una explicacin terica de los coe-
ficientes requerira de una mejor comprensin de los mecanismos de la tur-
bulencia,yaqueestndirectamenterelacionadosconlascaractersticas
dinmicas del flujo. En el captulo 26, se introdujeron dos posibles modelos
para explicar la transferencia convectiva de masa. Tanto la teora de pelcula
como la de penetracin se han aplicado ampliamente.
La teora de pelcula se basa en la presencia de una pelcula 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 pelcula ficticia en la qule l a transferencia se realiza

totalmente por medio de difusin molecular. El espesor de la pelcula, 6, debe
extenderse ms all de la subcapa laminar para illcluir la resistencia equiva-
lente que se encuentra cuando la concentracin c,ambia dentro de la capa de
amortiguamiento y el ncleo turbulento. En cuanto a la difusin a travs
de una capa que n o se est difundiendo, o fluido estancado, esta teora pre-
dice que el coeficiente de transferencia de masa es:.
(26-9)
tal como sevio en el captulo 26. El coeficiente de transferencia de masa,

correspondiente a la contradifusin 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 pelcula ficticia 6, nunca puede med.irse puesto que no existe.
Por esto y porque aparentemente no es adecuada para explicar fsicamente
la transferencia convectiva demasa, se han postulado otros modelos para
describir este fenmeno.
La teora de la penetracin fue expuestaoriginalmente,por Higbie*
para explicar la transferencia de masa en la fase lquida durante la absorcin
del gas. Danckwertst laha aplicado al flujo turbulento y tambin lo han
hecho otros investigadores cuando la componente en difusin slo penetra a
corta distancia enlafase de inters a causa de su desaparicin rpida por
medio de unareaccin qumica o del tiempo, re1ativa:mentecorto, de contacto.
Higbie supuso que la masa se transfera a la fase lquida por medio de la
transferencia molecular en estado no permanente. Con este concepto, se
expres el flujo de masa en la interfase lquido-gas, en la forma:
(26-86)
Danckwerts aplic este concepto de estado no permanente a la absorcin de

la componente A en una corriente turbulenta de lquido. Su modelo supone
que el movimiento del lquido ocasiona constantemente remolinos de lquido
nuevo del interior hacia la superficie, donde desplazan a los elementos de
lquido que se encontraban anteriormente en la superficie. Mientras se en-
cuentra en la superficie cada uno de los elementos del lquido queda expuesto
a la segunda fase y se transfiere masa al lquido 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 exposicin. Se pueden hacer muchas suposiciones respecto a la

renovacin de la superficie. Por ejemplo, cada uno de los elementos de lasu-
perficie puede tener el mismo tiempo de exposicin antes deser reemplazado.
Esto permite suponer que la transferencia instantnea de masa ocurrir de
acuerdo con la ecuacin (26-86). La cantidad total de soluto que penetra en
el remolino en un tiempo de exposicinigual a texp,es:
y la rapidez mediade transferencia durantela exposicinse obtiene dividiendo

esta ecuacin entre el tiempo de exposicin:
(28-62)
Danckwerts modific la suposicin de un periodo constante de exposicin,

proponiendo un conjunto infinito devalores de las edades de los elementos
de la superficie. Las funciones de distribucin 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-
vacin superficial se crea constante en un cierto grado de turbulencia e igual
a un factor S de renovacin superficial. La rapidez de transferencia de masa
con renovacin 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 renovacin superficial ha tenido mucho xito en
la explicacin y anlisis de la transferencia convectiva de masa, particular-
mente cuando la transferencia de masa se presenta acompaada de reacciones
qumicas en la fase lquida*?, pero se necesitan una considerable verificacin
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 quereaccionanqumica-
mente. En el estudio que se realiz en la seccin 26.2 acercade la transferencia
de masa molecular asociado con una reaccin qumica, se demostr que Ia
transferencia de masadepende de la constantederapidezde la reaccin
qumica. 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$ sealaron que es vlido el concepto de penetracin
de Danckwerts solamente cuando la renovacin 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).
Conclusin 687
mente rpida, proveyendo, as, la a superficie de elementos jvenes sobre una

base continua. Para los elementos de mayor edad de la superficie, se establece
un gradiente de concentracin en estado permanente, tal como lo predice la
teora de pelcula. Por esto, la rapidez de transferencia de masa debe ser
directamente proporcional a la difusividad molecdar de masa. Para valores
pequeos del nmero de Schmidtel gradiente perm.anente de la concentracin
se crea muy rpidamente en cualquier elemento nuevo de superficie de tal
modo que, a menos que la rapidez de renovacinsea lo suficientemente grande
para hacer desaparecer la enorme fraccin de los elementos de la superficie
antes de que se penetre en ellos, la mayor parte de los elementos de la superfi-
cie se comportarn como elementosviejos. AI aumentar el nmero de Schmidt,
el tiempo necesario para provocar un gradiente constante, aumenta rpida-
mente y por esto, es suficiente una rapidez de renovacin superficial pequea
para evitar que se penetre a la mayora delos elementos, tanto jvenes como
viejos, las caractersticas de transferencia son intermedias entre los modelos
de pelcula y penetracin. 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 estn basadas
en datos experimentales.
En ambos modelos: de pelcula y de penetracin, la transferencia de
masa comprende una interfase entre dos fluidos en movimiento. Cuando una
de las fases es un slido, 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 lmite para correlacionarlos
datos correspondientes a. un slido que se sublima convirtindose en gas o un
slido que se disuelve, transformndose en liquido. El coeficiente promedio
de transferencia de masa correspondiente a la difusin a travs de una capa
laminar lmite, es:
*"""
k, = 0.664-Re,"2Sc''3
DAB (28-21)
L
Esto demuestra que el coeficiente de transferencia de masa vara con
forma tpica de los c~lculosde capa lmite.
En este captulo se estudiaron los principios de transferencia de masa

por conveccin forzada, los parmetros importantes que ayudan en la des-
cripcin 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 conveccin est ntimamente relacionada con las carac-
tersticas dinmicas 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, energa y masa,

que vimos que se pudieron utilizar los mismos cuatro mtodos en la evaluacin
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 anlisis se correlacion el coeficiente de trans-
ferencia convectiva de masa por medio de la ecuacingeneral:
NuAB = f(Re, Sc)
La transferencia de masa hacia las corrientes turbulentas, tambinse estudi
y se defini la difusividad de remolino de la masa. Se presentaron analogas
correspondientes a la transferencia convectiva de masa hacia las corrientes
turbulentas.
PROBLEMAS
28.1 En una columna de atomizacin de transferencia de masa, un lquido

se atomiza hacia una corriente de gas y hay un intercambio de masa
entre las fases lquida y gaseosa. La masa de las gotas formadas en el
atomizador se considera que es funcin del dimetro de la tobera, de
la aceleracin de la gravedad, de la tensinsuperficial del lquido 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 deban haberse incluido algunas
otras variables?
28.2 Un cilindro largo de barro poroso, cuya concentracin 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 anlisis dimensional, que el perfil de
la concentracin dentro del cilindro debe expresarse en funcin de los
parmetros:
CA(r)-CA,m - k t
Y 2
CAo- CA,m r0 r0
28.3 La transferencia de masa en una capa lmite turbulenta formada sobre

una placa plana, est definida en funcin del nmero local de Nusselt,
por medio de:
N u A B , ~= 0.0292 Re,4'SS~1'3
Si la transicin 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 Transferencia de Materia

Cargado por

Copyright:

Formatos disponibles

LECTURA COMPLEMENTARIA 2 Transferencia de Materia

Cargado por

Información del documento

Descripción 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 Transferencia de Materia

Cargado por

Copyright:

Formatos disponibles

FUNDAMENTOS DE LA

LOS captulos anteriores relacionados con los fenmenos 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 concentracin, habr ms molculas

LA ECUACION DE RAPIDEZ DE FICK

Las leyes de transferencia demasa ponen de manifiesto la relacibn entre

Figura 24.1 Volumenelemental que contieneunamezcladecomponentesmltiples.

donde n es el nmero de especies presentes en la mezcla.La fraccin de masa,

Por defincin, la suma de las fracciones de masa, debe ser uno.

La concentracin molar de la especie A , C A , se define como el nmero

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 ecuacin (24-8) es una representacin algebraica de la ley de Dalton

En la Tabla 24.1 aparece un sumario de los diferentes trminos de con-

oxigeno, O,, yo, = 0.2 1

Determnese la fraccin de masa, tanto de oxgeno como de nitrgeno y el peso

oxgeno 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 = Fraccin de masa de la especie B = P B / P

Rlezcla lquida o slida Mezcla gaseosa

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

nitrgeno 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. Enunsistemadecomponentesmltiples, las diferentes

donde vi denota la velocidad absoluta de la especie i con relacin a ejes esta-

Esta es l a velocidad que medira un tubo de Pitot y la que se encontr

La velocidad de una especie particular con relacin a la masa promedio

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

vi-V. que es l a velocidad de difusin de l a especie i con relacin 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 concentracin total, c es constante bajo condiciones isotrmicas

donde d w , / d z es el gradiente de concentracin en funcihn de la fraccin 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 expresin enla relacin que se tena, se obtiene:

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

AI sustituir estos smbolos en la ecuacin (24-19), se obtiene una rela-

Esta relacin 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 difusin de A de la mezcla.

En condiciones isotrmicas e isobricas, esta relacinse reduce a:

'TABLA 24.2 FORMAS EQUIVALENTES DE LA ECUACION DE FLUJO

De acuerdo con la segunda ley de la termodinmica, los sistemas que nose

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

La ecuacin (24-24) se puede usar para definir todos los fenmenos de

dende p o es una constante: el potencial qumico del estado estndar. Cuando