Capítulo 7

FLUJO DE FLUIDOS POR LECHOS EMPACADOS

Diversas operaciones en la ingeniería requieren el flujo de fluidos a

través de lechos empacados con partículas de formas variadas. Es el caso
de la filtración, del intercambio iónico y las reacciones heterogéneas en
reactor catalítico. En otras ocasiones el flujo es bifásico como sucede en
la absorción líquido -gas. En este capítulo únicamente se estudia el flujo
de una sola fase fluida.

Los lechos empacados pueden ser fijos y fluidizados. En los primeros

el flujo es descendente o bien, si es ascendente su velocidad no es
suficientemente alta como para arrastrar el material de empaque. En los
lechos fluidizados en cambio, la velocidad del fluido ascendente es tal que
las partículas del lecho pierden contacto entre sí.

7.1. LECHOS FIJOS

El primer trabajo sobre lechos empacados fue realizado por Darcy en

Francia en 1830. Estableció que el gasto era proporcional a la caída de
presión e inversamente proporcional al espesor del lecho: Q (-P/L) lo que
se conoce como ecuación de Darcy. Esta relación lineal parece sugerir que
el flujo es laminar de acuerdo con la ecuación de Hagen -Poiseuille (2.19).

La ecuación de Fanning, aplicada a un ducto de sección no circular,

2 f  L
2

establece que − P = siendo: f=f(Re), Re=D e <>/ y D e el

De
diámetro equivalente, D e =4r h . <> es la velocidad media del fluido. El flujo
por un lecho empacado es mucho más complejo en cuanto que los pasajes
para el fluido son tortuosos, de diámetro no uniforme y no atra viesan todos
ellos el lecho. En el modelo del haz de tubos capilares la descripción
anterior se simplifica considerando el lecho como un haz de tubos rectos
que atraviesan todos ellos su espesor L. Se considera que el diámetro del
empaque es mucho menor qu e el de la columna que contiene el lecho y
que el diámetro de la columna es constante. Se trata ahora de poner la
ecuación de Fanning en función de los parámetros del lecho empacado.

Sea S c la sección transversal de la corriente, suma de las secciones

de todas las corrientes parciales a través del lecho. Sea b el perímetro
mojado por la misma. La definición de radio hidráulico es r h =S/b;
multiplicando y dividiendo por el espesor L del lecho se tiene:
 7-2

S  L  Volumen to tal de huecos

rh =  c   =
 b  L  área mojada de las partículas

Se supuso que todo el espacio entre las partículas está lleno de

líquido. Dividiendo ahora numerador y denominador por V (volumen total
del lecho) se obtiene
r h =/a s (7.1)

siendo a s el área mojada específica del empaque en metros cuadrados de

área mojada por metro cúbico de lecho empacado y  el grado de vacío o
porosidad, que es la fracción del volumen del lecho que no está ocupada
por las partículas. Este último p arámetro únicamente toma en
consideración lo s espacios interparticulares pero no los poros o cavidades
interiores de la partícula, cuando los hay.

Sea a e el área específica de partícula, igual a la relación entre su

área y su volumen: a e =A P /V P . Multiplicando y dividiendo esta relación por
el número total de partículas N se obtiene:

N  AP área superficia l partículas

ae = =
N  VP volumen partículas
área mojada partículas
as vol. lecho
La relación a s /a e es: = . Si se supone que toda
ae área sup. partículas
vol. partículas
el área superficial a s está mojada entonces la relación anterior es igual a
A
1- y por lo tanto, a s = P (1 −  ) . Introduciendo esta igualdad en la Ec. (7.1)
VP
 VP
se obtiene rh = y por lo tanto
1 −  AP
4 VP
De = (7.2)
(1 −  ) AP
Introduciendo la Ec. (7.2) en la ecuación de Fanning y en la definición
del reynolds se obtiene respectivamente:

f    2 L 1 −  AP
− P =   (7.3)
2  VP

 4VP   
Re =   (7.4)
 AP (1 −  ) 

 7-3

Por otro lado la velocidad medible en un lecho empacado no es < >

sino una velocidad media global llamada "superficial" definida por  s =Q/S.
Es conveniente entonces introducir esta va riable en las Ecs. (7.3, 7.4). Se
tiene por el balance de masa que: < >(sección vacía)=  s (sección total) y
(sección va cía )  L 
por lo tanto,  s =    =   o bien:
(sección to tal )  L 
s
 = (7.5)


En cuanto a la relación A P /V P , para partículas de forma geométrica

definida puede ser calculada directamente y usada en las Ecs. (7.3, 7.4).
As Ds 2 6
Así, para una esfera = = . No ocurre lo mismo con formas
Vs Ds / 6 Ds
3

irregulares. Para este caso se define, por semejanza con el de la esfera,

un diámetro de partícula que cumple la igualdad:

AP 6
= (7.6a)
VP DP

El diámetro D p puede ser expresado en función de la esfericida d. Se

define la esfericidad de una partícula,  s , como la relación entre el área de
una esfera y el de la partícula, ambas con el mismo volumen:
 s =(A s /V s )/(A P /V P ) de donde A P /V P =(A s /V s )/ s . De la igualdad anterior se
obtiene:

AP 6
= (7.6b)
V P Ds  s

Igualando las Ecs. (7.6a), (7.6b) se obtiene D P =D s  s siendo D s el

diámetro de la esfera del mismo volumen que la partícula o "diámetro de
la esfera equivalente"; es igualmente el diámetro nominal de la partícula.
Introduciendo las Ecs. (7.5, 7. 6a) en (7.3, 7.4) se obtienen las siguientes:

3 f s L(1 −  )
2

− P = (7.7)
DP  3

2 DP s 
Re = (7.8)
3 1−  

donde D P =6/a e =D s  s . Puesto que para una masa dada el área superficial
menor posible es la de una esfera, 0<  s < 1.0. En la Tabla 7.1 pueden verse
 7-4

algunos valores de la esfericidad. En cuanto a la porosidad del lecho

es muy difícil de predecir, sobre todo cuando se trata de una mezcla de
partículas de diversos tamaños; en cambio resulta muy fácil de determinar
experimentalmente. En la Fig. 7.1 se presenta una relación gráfica entre la
porosidad y la esfericidad. En la Tabla 7.2 se dan algunos valores de
porosidad.

Las condiciones extremas de aplicación de la Ec. (7.8) son el flujo

totalmente laminar y el flujo total mente turbulento.

Flujo laminar
La pérdida de presión debida a la fuerza de arrastre ejercida por el
fluido sobre las partículas, es en este caso directamente proporcional al
gasto (ecuación de Hagen -Poiseuille). En estas condiciones f=k 1 /Re.
Introduciendo en esta expresión el valor de Re dado por la Ec. (7.8) queda:
k (1 −  ) 
f = 2 (7.9)
DP s 

T abla 7.1. E s fer icid ad de d iver sas par tíc ulas ( 2 , 3 )

Mater ial esfer ici dad
pol vo de c arbó n, nat ural 0.65
pol vo de c arbó n, mol i d o 0.73
corcho 0.69
vi dri o mol i d o 0.65
hoj uel as d e mi ca 0.28
arena redo ndea da 0.83
arena angul ar 0.73
pol vo de t ungst eno 0.89
catal i zador de hi erro 0.58
carbón bi tumi n oso 0.625
Montur a B erl 0.30

Introduciendo ahora la Ec. (7.9) en (7.7) se obtiene:

− P k 3 s  (1 −  ) 2
=
L  3 DP 2
Utilizando múltiples datos experimentales Carman obtuvo: k 3 =150
quedando la ecuación:

− P 150 s  (1 −  ) 2
= (7.10)
L  3 DP 2

T abla 7.2. P or osid ad de le cho de diver sa s p ar tícula s( 1 , 4 )

Mater ial P or osida d
E sferas 0.412
C ubos 0.425
C i l i ndros, 3.1 75x3. 175 m m 0.401
 7-5

C i l i ndros, 3.1 75(di ám)x 6.35 mm 0.397

A ni l l os Lessi ng, 6 mm 0.87
A ni l l os R aschi g 1 cm 0.655
A ni l l os R aschi g, 2.4 c m 0.826
Montur as B erl , 1 cm 0.694
Montur as B erl , 1 .2 cm 0.71
P erdi gón de pl om o t amañ o uni for me,
0.635 c m 0.375- 0.42 1

La ecuación anterior es llamada de Carman -Koseny o Blake -Koseny.

D 
Es válida para Re = P s  10
(1 −  )

Figur a 7.1. R el aci ón en tre l a esferi ci d ad y l a p orosi dad ( 1 )

Flujo turbulento
La pérdida de presión en este caso es de bida principalmente a
efectos cinéticos asociados con remolinos, cambios de dirección,
contracciones y expansiones. Es por lo tanto proporcional al cuadrado del
gasto de acuerdo con la ecuación de Fanning (5.3). Para flujo altamente
turbulento, que sólo se puede dar en canales de mayor diámetro, el factor
de fricción, f T , depende principalmente de la rugosidad y no del número de
Reynolds, por lo que tenderá a ser constante como sucede en un tubo
(diagrama de Moody). Se puede sustituir 3f en la Ec. (7.7) por f T .
Experimentalmente se ha determinado que f T =1.75. La Ec. (7.7) queda para
el caso turbulento:

1.75 s L(1 −  )
2

− P = (7.11)
DP s  3
DP s 
que es la ecuación de Burke -Plummer, válida si:  1000 .
(1 −  )
 7-6

En realidad coexisten en diversas proporciones ambos flujos laminar

y turbulento. Si se supone que las pérdidas de presión por efectos viscosos
y cinéticos pueden sumarse se tiene una ecuación para el rango total de
flujo sumando las ecs(7.10, 7.11). Se obtiene así la llamada ecuación de
Ergun:

P 150 s  (1 −  ) 1.75  s (1 −  )
2 2

=− + (7.12a)
L DP  3
2
DP  3
o en forma adimensional:

− P DP  3 150
= + 1.75 (7.12b)
L  s 1 −  Re
2

siendo Re=D P  s /(1-) y D p =D s  s =6/a e . La ecuación de Ergun se ajusta

bien a los datos experimentales.

El producto D s  s es el único parámetro que caracteriza el tamaño y

la forma de una mezcla de partículas irregulares. Cuando las partículas
son redondeadas el diámetro D s puede ser determin ado por tamizado. Si
la forma de la partícula se aleja de la esférica el diámetro obtenido por
tamizado corresponderá a la segunda medida más grande de la misma. En
el caso de partículas aplanadas como hojuelas y discos el diámetro
obtenido será menor que D s mientras que si se trata de partículas
alargadas sucederá lo contrario.

Cuando el lecho está formado por partículas de diversos tamaños y

forma será necesario emplear en la ecuación de Ergun el valor medio D p =
Ds s . Para lechos de partículas de forma semejante y diversos tamaños Ds
1
puede ser calculado mediante la ecuación Ds = donde x i es la
( xi / d p ,i )
fracción en peso retenida en cada malla y d p =(d 1 +d 2 )/2, siendo d 1 , d 2 son
las aberturas de las mallas que limitan la fracción retenida. Si todas las
partículas quedan retenidas en la misma malla la ecuación anterior se
reduce a D s =(d 1 +d 2 )/2.

El diámetro de partícula puede también ser determinado

experimentalmente, a partir de mediciones de caída de presión en un lecho
empacado usando de preferencia un líquido, en condiciones en que Re<10,
por ser más precisa la Ec. (7.10).

Ejemplo 7.1. Una columna de 20 cm de diámetro se empaca con 4 kg de ácido

benzoico peletizado en partículas redondas de 6 mm de diámetro (Fig. 7.2). Se hace
pasar agua con objeto de obtener una solución del ácido. Si el gasto de agua es de 1.5
lps y el flujo es descendente a)¿Cuánto valdrá la caída de presión? b) Cuál será la lectura
 7-7

de un manómetro de mercurio? c) Calcular la potencia de bombeo. La densidad del

ácido es s = 1300 kg/m3. Considérese ’=1.0, =1 cp.

Solución:

0.0015
s = = 0.0477 m/s ; =0.412 (de Tabla 7.2)
(0.2)2 ( / 4)
 s =1.0; por la Ec. (7.12a):

− P 150 x0.0477 x10 −3 (1 − 0.412) 2 1.75 x1000 x(0.0477) 2 (1 − 0.412)

= + =6,562 Pa/m
L 0.006 2 x0.412 3 0.006 x0.412 3

Volumen sólido:4/1300=3.077x10 - 3 m 3
Volumen del lecho: (3.077x10 - 3 )/(1-0.412)=5.23x10 - 3 m 3 . S=0.0314
m2
a) L=5.23x10 - 3 /0.0314=0.17 m; -P=6562x0.17=1116 Pa

b) Por la Ec. (5.46): 1116=9.81h m (13590-1000); h m =0.009 m Hg.

c)Potencia consumida: Pot=1116x0.0015=1.674 W

Figur a 7.2. E j em pl o 7. 1.

Ejemplo 7.2 . A un lecho formado de partículas de catalizador de

fierro se introduce una mezcla gaseosa con peso molecular M= 36.5, a
193,000 Pa, 423K. A la salida M=38.2 y T=468K. La profundidad es de 2
m. Las partículas son semejantes en forma; su composición granulométrica
es dada en la Tabla 7.3, cols. 1,2. El flujo es de 2600 kg/m 2 h, descendente.
a) Calcular -P/L. b) Si el área transversal es de 0.15 m 2 ¿cuánto vale la
potencia disipada?. =3.76x10 - 5 Pas.

T abla 7.3. E jem p lo 7. 2

 7-8

A ber tur a de m alla, m m fr acció n r ete nida  d1 + d 2 

(d 1 , d 2 ) xi xi/  
 2 
0.62-0 .67 0.4 0.62
0.67-0 .72 0.4 0.57
0.72-0 .76 0.2 0.27
S uma: 1.46

Solución:
Ds =1/1.46=0.685 mm= 6.85x10 - 4 m
De la Tabla 7.1  s =0.58; de la Fig. 7.1 =0.56;
P 1 =193,000+101,300=294,300 Pa
 1 =(294300x36.5)/(8314x423)=3.05 kg/m 3 .  s =2600/(3600x3.05)=0.237
m/s. Por la Ec. (7.12a):

− P 150 x0.237 x(3.76 x10 −5 )(1 − 0.56) 2 1.75 x3.05 x0.237 2 (1 − 0.56)
= + =11,226 Pa/m.
L (6.85 x10 −4 ) 2 x0.58 2 x0.56 3 (6.85 x10 −4 ) x0.58 x(0.56) 3
a) -P=22,452 Pa. Siendo la caída de presión menor de 10% de la
original, la consideración de fluido incompresible está justificada.
b) Pot=22452x0.237x0.15=798 W

7.2. LECHOS FLUIDIZADOS

Si a través de un lecho empacado con partículas se hace pasar un fluido en forma

ascendente y con gasto creciente, al principio el lecho se comporta como fijo pero al
aumentar el gasto se expande, al igualarse la fuerza de arrastre con el peso aparente de
las partículas. Si el gasto continúa aumentando las partículas se separan; se dice que el
lecho está fluidizado. El lecho se comporta entonces como un verdadero fluido: su
superficie es horizontal, puede caer por un vertedero, salirse por una abertura y un
cuerpo ligero flota en él; los puntos situados en el mismo plano horizontal tienen la misma
presión de acuerdo con la ley hidrostática.

Al aumentar el flujo más allá del punto de fluidización, aumenta la porosidad del
lecho manteniéndose la caída de presión en un valor aproximadamente constante que
equilibra el peso aparente de las partículas según la ecuación:

− P = (1 −  )( s −  f )gL (7.13)

Si la velocidad continúa aumentando el comportamiento cambiará en una forma

que dependerá de la naturaleza del fluido. En el caso de un líquido el lecho continuará
siendo homogéneo al aumentar el flujo; se dice que es un lecho particulado o
simplemente fluidizado líquido. Si se trata de un gas, en cambio, el aumento del flujo
determina la formación de burbujas que se desplazan dentro de la emulsión sólido-gas,
la cual constituye la fase densa o continua en este sistema bifásico. La fase ligera o
discontinua es la gaseosa, que se segrega parcialmente en burbujas. Dentro de la fase
densa la velocidad relativa entre partículas y gas no aumenta prácticamente con el gasto,
 7-9

e incluso puede mantenerse la condición laminar. Este tipo de lecho fluidizado se

llama agregativo o burbujeante o simplemente lecho fluidizado gaseoso. En raras
ocasiones el lecho fluidizado líquido es agregativo o el gaseoso particulado.

La superficie del lecho fluidizado permanece definida en tanto que la velocidad

del fluido no supere la velocidad terminal de las partículas, caso en el cual su frontera se
hace indefinida y comienza el arrastre o transporte del sólido. Velocidad terminal es la
que alcanza una partícula que cae en un fluido quieto. No existen en la literatura abierta
correlaciones fiables para manejar la fluidización gaseosa. Para la fluidización líquida en
cambio, aun cuando también existe la tendencia a la formación de burbujas, se dispone de
ecuaciones empíricas.

Graficando -P v.s. s en escala doble logarítmica se obtiene la curva de la Fig. 7.3.
Se observa que -P llega a un valor en el que permanece constante. La variación de
condiciones representada en ella presenta un fenómeno de histéresis al reducirse
gradualmente el gasto. En efecto, la caída de presión resulta menor que en la curva de
gasto ascendente debido a que el movimiento del fluido ya ordenó las partículas de forma
tal que ofrezcan una menor resistencia al paso del fluido.

Entre las ventajas que el lecho fluidizado presenta con relación al fijo están la
facilidad de control de temperatura que permite mantener isotérmico el sistema, así como
las elevadas velocidades de transferencia de masa alcanzadas. Entre sus desventajas
están la pulverización de partículas friables y la ya mencionada dificultad en el tratamiento
matemático del caso gaseoso. El lecho fluidizado es adecuado para operaciones en gran
escala. Ha adquirido interés práctico en relación con los reactores catalíticos,
principalmente cuando la reacción es muy exotérmica. En este caso está la
descomposición térmica de hidrocarburos pesados, la síntesis de combustibles a partir
del gas natural y otras reacciones heterogéneas. La fluidización y el arrastre se usan
también en el transporte neumático de sólidos, en intercambio de calor, secado, adsorción,
etc.

Velocidad mínima de fluidización

Es el valor de la velocidad superficial a la cual comienza el fenómeno de fluidización.
Existen algunas discrepancias entre el valor calculado y el medido. Una ecuación
recomendada para calcular esta velocidad es (3)

(0.0408)Ds 3  f ( S −  f )g 
1
Ds mf  f 2

Re mf = = (33.7 ) +  − 33.7
2
(7.14)
  2 
6
Para velocidad terminal se tiene la ecuación ( ):

 s  ( s −  f )gDs
2

T = (8.43x10 −3
)  ln   (7.15)
−4
 6.5 x10  18
 7-10

Figura 7.3. Caída de presión en un lecho empacado con flujo ascendente

La relación mf/T tiene valores que van de 9 para sólidos gruesos (Re P>1000) a
90 para sólidos finos (ReP<0.4). Para la porosidad en función del gasto puede usarse la
ecuación de Richardson y Zaki (en ref. 6):

1
  n
 =  s  (7.16)
 T 
donde n es dada por las siguientes ecuaciones:

 D  −0.03
n =  4.35 + 17.5 s  Re si 0.2<Re
 DT 

 D  −0.1
n =  4.45 + 18 s  Re si 1<Re<200
 DT 
−0.1
n = 4.45 Re si 200<Re<500
n = 2.39 si 500<Re

Ds s  f
En las ecuaciones anteriores Re =

Ejemplo 7.3. Sea nuevamente el caso del ejemplo 7.1 . Si el gasto de agua es
ahora ascendente de 1.5 lps, ¿cuánto valdrá la caída de presión y a qué valor del gasto
se efectuará el arrastre del empaque?.

Solución:

1. Velocidad superficial. Ya fue calculada en el ejemplo 7.1: s=0.0477 m/s

2. Velocidad mínima de fluidización:

 7-11

1
 1300 − 1000  2
Re mf = 33.7 2 + 0.0408(0.006 ) x1000 9.81 − 33.7 = 130.83
3

 (
10 −3
2
) 
de donde

130.83x10 −3
 mf = = 0.0218 m/s; Puesto que s>mf el lecho está fluidizado
1000 x0.006

3. Velocidad terminal:

  (1300 − 1000 )x9.81(0.006 )

 2
1
 T = 8.43x10 ln 
−3
  = 0.364 m/s
 6.5 x10
−4
 
18 x0.001 
Puesto que  s   T el lecho no es arrastrado por el líquido.

4. Cálculo de la porosidad:

0.006 x0.0477 x1000

Re P = −3
= 286 ; n = 4.45 x 286 −0.1 = 2.527 . La porosidad resultante
10
1
 0.0477  2.527
es  =   = 0.447
 0.364 

5. Profundidad del lecho:

Masa 1 1
L= , donde S es el área seccional. Introduciendo valores:
s 1−  S
4 1 4
L= = 0.177 m
1300 1 − 0.447  (0.2 )2
6. Caída de presión:

-P=(1300-1000)(1-0.447)x9.81x0.177=288 Pa

El lecho será arrastrado cuando s>0.364 m/s o sea Q>0.364x0.0314=0.0114

3
m /s

Ejercicios (ver respuestas en el Apéndice D)

1.-En un tubo vertical de 8" std hay un lecho de partículas cuya forma es cercana a la cúbica, con
2 mm de lado. El lecho ocupa 30 cm del tubo. Por el mismo pasa agua a temperatura ambiente a razón
de 0.9 lps. La densidad del sólido es '=1.4. Calcular la caída de presión cuando el flujo es descendente.
2. Se quiere determinar el diámetro de partícula de un sólido que va a ser utilizado en un lecho
fluidizado. Para ello se empaca un tubo que tiene 25 cm de diámetro interior con el material sólido,
obteniéndose un lecho con 3 m de profundidad. Por el mismo se hace pasar agua a la temperatura
ambiente a razón de 1.2x10-4 m3/s. La caída de presión producida es de 2670 Pa. =0.38. Calcular D p y
comprobar que Re<10.
 7-12

3. En un tubo con 15 cm de diámetro interior hay 6 kg de catalizador sólido con la siguiente

granulometría: 53% se retiene entre los tamices con abertura 2.362-1.651 mm, y el resto entre los
tamices 1.651-1.168 mm Calcular el flujo mínimo de líquido para que el lecho se fluidice. Las propiedades
del líquido son: ’=0.89, =0.7 cp. La densidad relativa del sólido es 2.7.

Bibliografía

1. Brown,G.G. et al, Unit, Operations, John Wiley and Sons, New York, 1950
2. Perry R.,Green,D.,Perry's Chemical Engineers' Handbook, 6 th Ed, McGraw-Hill International Editions,
1984
3. Kunii, D., Levenspiel, O. Fluidization Engineering. Robert E.Krieger Publishing Co. Inc. Florida, 1962
4. Coulson,J.M.,Richardson,J.M.,Chemical Engineering, Pergamon Press, 1978.
5. Bird,R.B.,Stewart,W.,Lightfoot,E.,Transport Phenomena. John Willey,1960.
6. Leva,M., Fluidization, McGraw-Hill Book Company, New York, 1959

Notación

Ap área superficial de 1 partícula, m2

ae área específica de partícula, ae=Ap/Vp, m-1
as área específica del lecho, m2/m3 de lecho
d1, d2 abertura de malla, m
D diámetro de clilindro, m
dp diámetro de partícula definido por d p =(d 2 +d 2 )/2, m
Dp diámetro de partícula definido como DP=6/ae=DSs, m
Ds diámetro de la esfera equivalente; diámetro nominal de la partícula, m
DT diámetro del lecho empacado (sección circular), m
Frmf froude mínimo de fluidización = mf2/DPg, adimensional
h altura del cilindro, m
Remf reynolds mínimo de fluidización (adim)
S área seccional del lecho, m2
Sc área seccional total de las corrientes de fluido en el lecho, m2
mf velocidad mínima de fluidización, m/s
Vp volumen de partícula, m3
s velocidad superficial ms-1
T velocidad terminal, ms-1
<> velocidad media, ms-1
x fracción (en peso) retenida en una malla, adimensional.

 porosidad
s esfericidad

subíndices
f fluido
mf mínimo def luidización
s sólido; superficial; esfera equivalente
p partícula
7-13