Circulación de fluidos a través de

lechos porosos
Ing. Jesús Alexander Sánchez González

El flujo a través de lechos de

partículas sólidas aparece en
numerosos procesos técnicos:
Filtración, flujo a través de
columnas de relleno
(Destilación, Absorción,
Adsorción e Intercambio Iónico)

Ingeniería
de
Alimentos
II
 Permeabilidad (Ley de Darcy)
 Darcy demostró que la velocidad
media, en una sección del lecho, era
Ing. Jesús Alexander Sánchez González

directamente proporcional a la caída

de presión que experimenta el fluido al
atravesar dicho lecho e inversamente
al espesor del mismo.

 (-∆p)= caída de presión a través del

lecho
 L= Espesor o altura del lecho
 K= Constante de proporcionalidad
(depende de las propiedades físicas del lecho y
del fluido que circula a través del lecho)

Ingeniería
de
Alimentos
II
 …Permeabilidad (Ley de Darcy)
 Suponiendo un régimen laminar (ya que la
velocidad a través de los intersticios del lecho granular es
baja y además la sección de cada uno de ellos es pequeña: por
Ing. Jesús Alexander Sánchez González

lo que el valor del módulo de Reynolds no es elevado)

la
relación entre la velocidad de circulación
del fluido y la caída de presión que
experimenta es lineal.

 Todo ello presupone que la resistencia

que ofrece el lecho al flujo de fluido es
debida principalmente a rozamientos
viscosos

Ingeniería
de
Alimentos
II
 …permeabilidad
La viscosidad (μ)
del fluido muestra
Ing. Jesús Alexander Sánchez González

la resistencia
intrínseca del
fluido
La constante
(1/α) recibe el La unidad de
nombre de permeabilidad es
coeficiente de el Darcy
permeabilidad

Ingeniería
de
Alimentos
II
 Superficie específica
También denominada área superficial específica, puede ser referida a
todo el lecho o a la partícula
Ing. Jesús Alexander Sánchez González

Superficie específica del lecho (aS)

Área superficial específica de las

partículas (aS0)

En el caso de una partícula esférica de diámetro “dr”, la superficie

específica es aS0=6/dr

Cuando las partículas no son esféricas, se define un diámetro equivalente de partícula (dp),
como aquel diámetro que poseería una esfera cuya relación área superficial a su volumen
fuese la misma que la que posee una partícula Ingeniería
de
Alimentos
II
 Forma de la partícula Esfericidad (Γ)

Esfericidad de Esfera 1

partículas Cubo 0.81

Cilindros
•h= d 0.87
Ing. Jesús Alexander Sánchez González

•h=5d 0.70
•h=10d 0.58
Discos
•h= d/3 0.76
•h=d/6 0.60
•h=d/10 0.47
Arena de playa Tan alta como 0.86

Arena de río Tan baja como 0.53

Distintos tipos de arena 0.75

Sólidos triturados 0.5 – 0.7

Partículas granulares 0.7 – 0.8

Trigo 0.85

Anillos Rasching 0.26 – 0.53

Sillas Berl 0.30 – 0.37 Ingeniería

de
Alimentos
II
 Porosidad (ε)
Porosidad o fracción de huecos
• El lecho de partículas no es compacto
Ing. Jesús Alexander Sánchez González

• Debida a la porosidad del lecho, la superficie específica del lecho y de

partícula no coinciden, sino que están relacionadas según la ecuación:

Es fácil observar que cuando la porosidad aumenta, el flujo a

través del lecho es mucho mejor, con lo que la permeabilidad
aumenta

En el estudio de las características del lecho, se va a suponer

que:
• Está compuesto por partículas dispuestas al azar, de forma que en el
interior de ellos se forman canales por los que circula el fluido,
• La longitud de estos canales (L’) es la misma para todos ellos, poseyendo el
mismo diámetro equivalente (De) Ingeniería
de
Alimentos
II
 Variables características de un lecho

n’ • Número de canales por metro cuadrado de

sección transversal del lecho
Ing. Jesús Alexander Sánchez González

• Número total de canales del lecho (n)

• Área interfacial de un canal

• Área interfacial del lecho

• Volumen del lecho (VL)

• Volumen de lecho ocupado por las partículas

Ingeniería
de
Alimentos
II
 Valores de las superficies específicas
Ing. Jesús Alexander Sánchez González

 Correlación entre superficie específica del lecho

y la de la partícula

Ingeniería
de
Alimentos
II
 Características de un canal
La velocidad de un canal es función de la velocidad global y de las
longitudes del lecho y del canal, así como de la porosidad
Ing. Jesús Alexander Sánchez González

• Diámetro equivalente de un canal (De)

• Expresado en términos de
lecho poroso

• El diámetro equivalente de un canal en

función de ε y de aS
• Ec. de la continuidad; donde SC es la sección de
paso de todos los canales, vC es la velocidad de
circulación del fluido a través de un canal

Ingeniería
de
Alimentos
II
 Ecuaciones para el flujo a través de lechos porosos

• Régimen Laminar
Ing. Jesús Alexander Sánchez González

• Ecuación de Kozeny-Carman

• Régimen Turbulento
• Ecuación de Burke-Plummer

• Flujo global laminar-

turbulento
• Ecuación de Ergun
Ingeniería
de
Alimentos
II
 Flujo en Régimen Laminar
 Si el flujo del fluido a través de los canales es laminar se podrá
aplicar la ecuación de Fanning para cada uno de ellos:
Ing. Jesús Alexander Sánchez González

 Pero si se tiene presente que el módulo de Reynolds es:

 Y que la velocidad de circulación del fluido a través de una canal vC

viene dada por la siguiente expresión

 Se obtiene
Ingeniería
de
Alimentos
II
  Se observa que la pérdida de presión depende, entre otros, de las
longitudes de cada canal y del lecho.
Ing. Jesús Alexander Sánchez González

 La longitud de cada canal L’ es superior a la del lecho.

 Si se supone que dichas longitudes son proporcionales L’=K’L, y
definiendo una constante K’’=2(K’)2, resulta:

 De la que se obtiene la Ec. de Kozeny-Carman

Ingeniería
de
Alimentos
II
  Al comparar la ecuación de Kozeny-Carman con la de Darcy, se
obtiene que la permeabilidad será:
Ing. Jesús Alexander Sánchez González

 La constante K’’ se denomina constante de Kozeny.

 En lechos en los que la porosidad y superficie específica no varían
con el espesor del lecho, se han encontrado experimentalmente
que esta constante posee un valor de 5±0,5.
 En realidad su valor depende del tipo de relleno y de la porosidad,
tomando distintos valores según sea la forma de las partículas y
porosidad del lecho.
 En el caso que las partículas sean de forma esférica, el valor de esta
constante es de 4,8±0,3. Ingeniería
de
Alimentos
II
 La constante de Kozeny (K’’)

 El valor de la constante de Kozeny no es el mismo para todos los

tipos de relleno, sino que depende de la relación L’/L
Ing. Jesús Alexander Sánchez González

 Carman ha demostrado que:

◦ La relación L’/L recibe el nombre de tortuosidad
◦ K0 es un factor que depende de la sección transversal del canal.

 A pesar de que la tortuosidad y el factor K0 pueden variar, esta

variación es tal que cuando una aumenta la otra disminuye, y
viceversa, de forma que su producto posee valores próximos a 5

Ingeniería
de
Alimentos
II
 Correcciones por efectos de la pared
 En los lechos de partículas, las que se hallan en contacto con la
pared están menos compactadas, lo que implica que la resistencia
Ing. Jesús Alexander Sánchez González

que ofrece el lecho al flujo es menor que el dado por la ecuación

de Kozeny-Carman
 Coulson ha obtenido un factor de corrección KP de modo que se
tenga presente este defecto.

◦ En la que AP es la superficie de la pared de la columna que contiene el

lecho por unidad de volumen de dicho lecho
 Para el cálculo de la caída real de presión, deberá multiplicarse la
calculada a partir de la ecuación de Kozeny-Carman por el valor de
este factor KP

Ingeniería
de
Alimentos
II
 Flujo en Régimen Turbulento
 Si el régimen de circulación es turbulento puede pensarse en
aplicar la ecuación de Fanning a la circulación del fluido por un
Ing. Jesús Alexander Sánchez González

canal, así:

 Teniendo en cuenta las expresiones del diámetro equivalente y de

la velocidad por un canal:

 Se obtiene lo siguiente

Ingeniería
de
Alimentos
II
  Al igual que antes, si se supone que L’=K’L, y que la superficie
específica esta relacionada con el diámetro de la misma por la
Ing. Jesús Alexander Sánchez González

expresión aS0=6/dP se obtiene:

 Si se define un factor de fricción modificado f’=(K’)3f se obtiene la

ecuación de Burke-Plummer:

 El valor de se obtiene a partir de experimentación, dependiendo

del número de Reynolds
Ingeniería
de
Alimentos
II
 Flujo Global Laminar-Turbulento
 Observando las Ecuaciones de Kozeny-Carman y de Burke-
Plummer, puede pensarse que la pérdida de presión por unidad de
Ing. Jesús Alexander Sánchez González

longitud del lecho será una expresión del tipo:

 O bien una combinación lineal de dichas ecuaciones:

 Para la obtención de las distintas constantes α’ y β’ que aparecen

en esta ecuación se recurre a datos obtenidos por
experimentación, ajustándolos a ella
Ingeniería
de
Alimentos
II
 Número de Reynolds
 En la circulación de fluidos por el exterior de sólidos, no sólo se
utiliza el módulo de Reynolds para un canal, sino que en algunos
casos se suele modificar, o bien se utiliza el módulo de Reynolds de
Ing. Jesús Alexander Sánchez González

partícula
◦ Reynolds para 1 canal:

◦ Reynolds modificado:

◦ Reynolds de partícula:

 Régimen de circulación en función del ReP

◦ (ReP<40) el régimen de circulación se considera laminar
◦ (ReP>40) el régimen se considera turbulento. Ingeniería
de
Alimentos
II
 Factor de fricción (f ’)

 De la Ecuación
Ing. Jesús Alexander Sánchez González

 Es fácil obtener la expresión para el factor de fricción modificado:

 A partir de una serie de datos experimentales, de valores del factor

de fricción modificado [3f’ ], calculados a partir de la ecuación de
Burke-Plummer, se obtiene los valores correspondientes de sus
Reynolds de partícula.
 Se busca una función que correlacione 3f’ con ReP
 Del ajuste de los datos experimentales la ecuación obtenida es la
siguiente:

Ingeniería
de
Alimentos
II
  Si en esta ecuación se sustituyen los valores de 3f’, dados por las
ecuaciones que expresan al factor de fricción y la correspondiente
Ing. Jesús Alexander Sánchez González

al módulo de Reynolds de partícula,

 se puede obtener la siguiente expresión denominada ecuación de

Ergun, y puede utilizarse para el cálculo de la pérdida de presión
que experimenta el fluido al atravesar un lecho relleno,
independientemente de cual sea el tipo de régimen de circulación

Ingeniería
de
Alimentos
II
 Comparaciones finales
Ing. Jesús Alexander Sánchez González

 Al comparar la Ecuación de Ergun con la linealizada inicial

 se observa lo siguiente:

Ingeniería
de
Alimentos
II
 CASO ESPECIAL: Flujo en contracorriente de gas y líquido

 Cuando en una columna rellena de partículas circulan en

contracorriente un gas con un líquido, es conveniente utilizar la
denominada Ecuación de Chilton-Colburn, que es una ecuación
Ing. Jesús Alexander Sánchez González

empírica basada en la de Fanning

 El factor de fricción modificado (f’) puede calcularse a partir de las

graficas que relacionan

 o bien pueden utilizarse las siguientes expresiones:

◦ Para régimen Laminar (ReP<40):

◦ Para régimen turbulento (ReP>40):

Ingeniería
de
Alimentos
II
 Condiciones para el uso de la Ec. Chilton-Colburn
 Se puede utilizar cuando el lecho está formado por partículas macizas
 Si las partículas son huecas debe multiplicarse el segundo miembro
Ing. Jesús Alexander Sánchez González

de la Ecuación por un factor Kr que viene dado por la ecuación:

◦ dP es el diámetro nominal de las partículas expresado en pulgadas.

 Los huecos creados en las inmediaciones de la pared pueden afectar
la caída de presión, para lo que se introduce el factor de corrección
de pared KP
◦ KP: depende del régimen de circulación y de la relación diámetro partícula a
diámetro del lecho (dP/D)
◦ En el caso de que: dP/D < 1,6 no se considera el efecto de la pared.
 En general la caída real que experimenta el fluido se obtendrá multiplicando la caída de presión
calculada a partir de la ecuación de Chilton-Colburn por estos dos factores; es decir:

Ingeniería
de
Alimentos
II
 FLUIDIZACIÓN
Ing. Jesús Alexander Sánchez González

 Cuando un fluido circula a través de un lecho

de partículas:
◦ si su velocidad de circulación es baja, el lecho
permanece estático,
◦ pero si su velocidad aumenta puede ocurrir que el
lecho se expanda, existiendo una reordenación de las
partículas, con el consiguiente aumento de la
porosidad del lecho

Ingeniería
de
Alimentos
II
 Si se supone un lecho de partículas, a través de las cuales circula un fluido
al que se le va aumentando paulatinamente la velocidad de circulación.

 A velocidades bajas el lecho permanece estático, sin aumento de su altura,

por lo que al representar en coordenadas doble logarítmicas la caída de
Ing. Jesús Alexander Sánchez González

presión frente a velocidad de circulación se obtiene una recta tal como lo

muestra el segmento 0A de la gráfica

log(-∆p)

D
B
A
C

Fluidización Fluidización
discontinua continua

0 vmf va log (v) Ingeniería

de
Alimentos
II
  Si se aumenta más la velocidad, las partículas empiezan a separarse, pero todavía
permanecen en contacto, aunque la relación entre la caída de presión y la velocidad
continua siendo lineal, pero la pendiente es menor (segmento AB).
 En el punto B en el que las partículas no se hallan en contacto se dice que el lecho es
fluidizado.
 A partir de este punto, si se aumenta la velocidad puede haber una pequeña caída de
presión (segmento BC), pero si se sigue aumentando la velocidad, la caída de presión
Ing. Jesús Alexander Sánchez González

aumentara linealmente con la velocidad, esta vez con menor pendiente, hasta que la
velocidad es lo suficiente mente elevada como para arrastrar las partículas, cosa que
ocurre a partir del punto D

log(-∆p)

D
B
A
C

Fluidización Fluidización
discontinua continua

Ingeniería
de
0 vmf va log (v) Alimentos
II
 Características de la fluidización
 La velocidad a la cual ocurre la fluidización se conoce como:
mínima de fluidización;
Ing. Jesús Alexander Sánchez González

 La velocidad correspondiente al punto de arrastre es la velocidad

de arrastre.
 En el tramo en que el lecho esta fluidizado, la velocidad del fluido
en contacto con las partículas es superior a cuando las ha
superado, ello hace que las partículas no queden arrastradas, sino
que vuelvan a caer al lecho.
 El conjunto de partículas poseen un movimiento desordenado,
presentando un aspecto como si las partículas estuvieran en
ebullición, denominándose lechos hirvientes o fluidización
discontinua a este tipo de proceso.
 Cuando las partículas son arrastradas, el tipo de fluidización se
denomina continua, y es un tipo de circulación de dos fases, que es
la base del transporte neumático Ingeniería
de
Alimentos
II
 VELOCIDAD MINIMA DE FLUIDIZACION

 La velocidad a la que empieza la fluidización del lecho se conoce

como velocidad mínima de fluidización. En este punto existe un
Ing. Jesús Alexander Sánchez González

equilibrio dinámico entre la fuerza que el campo gravitatorio y el

fluido ejercen sobre las partículas
 La fuerza gravitatoria ejercida sobre las partículas viene dad por la
expresión:

◦ En la que ρP y ρ son las densidades de las partículas y del fluido, respectivamente;

◦ S la sección de paso de la columna que contiene las partículas;
◦ L la altura del lecho
◦ ε su porosidad
◦ g la constante de gravedad

Ingeniería
de
Alimentos
II
 Relación de fuerzas

 La fuerza que el fluido ejerce sobre el lecho de partículas es la de

presión, y que se calcula por la expresión
Ing. Jesús Alexander Sánchez González

◦ En la que (-Δp) es la caída de la presión que experimenta el fluido al atravesar el

lecho, y cuya expresión depende del régimen de circulación del fluido
 Por tanto, para el cálculo de la velocidad mínima de fluidización
deben igualarse las so fuerzas, la de gravedad y la de presión

Ingeniería
de
Alimentos
II
 La expresión de la caída de presión viene dad por
la ecuación de Ergun
Ing. Jesús Alexander Sánchez González

 A pesar que esta es la ecuación general para el cálculo de la caída

de presión, según sea el régimen de circulación del fluido, podrá
simplificarse.
Régimen
laminar

Régimen
Turbulento

Régimen de
Transición Ingeniería
de
Alimentos
II
 POROSIDAD MÍNIMA DE FLUIDIZACIÓN
Ing. Jesús Alexander Sánchez González

 Lo mismo que ocurría con la caída de presión del fluido, que

variaba con la velocidad de circulación, la porosidad del lecho
también experimenta variación.
 Para velocidades bajas, la fuerza de presión no es suficiente para
variar la estructura del lecho, y su porosidad no varia, pero a
velocidades más altas, el lecho se expansiona aumentando el
volumen de huecos.

Ingeniería
de
Alimentos
II
 Si se representa en coordenadas doble logarítmica las variaciones de la
porosidad con la velocidad lineal de circulación se obtiene una gráfica
como por ejemplo

log(ε)
Ing. Jesús Alexander Sánchez González

εmf

log (v)
0 vmf va
Ingeniería
de
Alimentos
II
 Cálculo de la porosidad

 A velocidades bajas la porosidad no varía, pero llega un momento

en que aumenta con la velocidad, siendo esta variación tipo lineal.
Ing. Jesús Alexander Sánchez González

La porosidad en el momento que empieza la fluidización no se

corresponde al punto en que la porosidad empieza aumentar, sino
que ocurre a una velocidad más alta, que es la mínima de
fluidización.

◦ Ecuación en la que el diámetro de partículas debe expresarse en micrones (10 -

6m) y es válida para valores de d comprendidos entre 50 y 500 micrones
P

Ingeniería
de
Alimentos
II
 ALTURA DEL LECHO
 Evidentemente, al aumentar la velocidad de circulación del fluido a través del
lecho poroso, no sólo aumenta la caída de presión, sino que el aumento de la
porosidad es debido a que la altura del lecho también aumenta.
Ing. Jesús Alexander Sánchez González

 Este aumento de altura de lecho esta íntimamente ligado al aumento de

porosidad, de forma que las porosidades correspondientes a dos alturas
cualesquiera están relacionadas por la expresión

 De forma particular, si para el lecho fijo le corresponde una porosidad ε0 y una

altura de lecho L0, la altura y porosidad para otro instante están relacionadas
con ésta por la ecuación:

 Si el lecho fuese compacto, las partículas ocuparían todo el lecho y no existirían

huecos, por lo que la expresión que correlaciona las alturas se simplificaría, en
la que es la altura del lecho compacto:
Ingeniería
de
Alimentos
II
 Ing. Jesús Alexander Sánchez González

II
de
Alimentos
Ingeniería