FLUIDODINÂMICA EM

SISTEMAS PARTICULADOS

Giulio Massarani
Programa de Engenharia Química
COPPE/Universidade Federal do Rio de Janeiro

Versão Preliminar da 2º Edição

2001
 Ao amigo José Teixeira Freire

2
 SUMÁRIO
Prefácio 5
Capítulo 1 Fluidodinâmica da Partícula Sólida 6
1. Equação do Movimento da Partícula 6
2. A Força Resistiva Fluido-Partícula 11
Efeito da presença de fronteiras rígidas 17
Influência da concentração de partículas 20
3. O Movimento Acelerado da Partícula 26
4. Dinâmica da Partícula em Fluido Não-Newtoniano 28
Problemas 31
Bibliografia 36
Capítulo 2 A Decantação 38
1. A Trajetória da Partícula 38
2. Separação Sólido-Fluido na Fenda de Seção Retangular 42
3. O Conceito Sigma e a Especificação de Centrífugas 44
4. Ciclones a Gás e Hidrociclones 45
Problemas 53
Bibliografia 60
Capítulo 3 Escoamento de Fluidos em Meios Porosos 62
1. Equações da Continuidade e de Movimento para o Fluido 62
A força resistiva m 64
A tensão extra τ 65
A equação de Darcy 65
2. Propriedades Estruturais da Matriz Porosa 66
A determinação experimental de parâmetros estruturais 66
O modelo capilar 67
3. Escoamento em Meios Porosos: Aplicações Clássicas 73
A perda de carga no meio poroso 73
O escoamento compressível 75
O escoamento transiente 75

3
 4. O Escoamento Bifásico em Meios Porosos 75
Equação de Darcy-Buckingham 76
Generalização da Forma Quadrática de Forchheimer 77
Problemas 80
Bibliografia 95
Capítulo 4 Fluidodinâmica em Sistemas particulados Expandidos 97
1. Equações da Continuidade e do Movimento 97
2. Caracterização dos Meios Expandidos 101
3. O Elo entre a Fluidodinâmica de Partículas e a Teoria de Misturas 103
4. Transporte Hidráulico e Pneumático de Partículas 105
Transporte vertical homogêneo: partículas "grandes" 106
Transporte hidráulico homogêneo 108
Problemas 109
Bibliografia 117
Capítulo 5 Escoamento em Meios Porosos Deformáveis 118
1. Equações da Continuidade e do Movimento 118
2. Teoria da Filtração com Formação de Torta 122
Equacionamento da filtração plana com formação de torta 123
A teoria simplificada da filtração 126
3. A Sedimentação Contínua 128
Problemas 131
Bibliografia 145
Índice Onomástico 146

4
 PREFÁCIO
Primeira Edição (Ed. UFRJ, 1997)

Entre as múltiplas facetas que os Fenômenos de Transporte em Sistemas Particulados oferecem, tanto
do ponto de vista científico como numa larga gama de aplicações tecnológicas, este livro trata apenas dos
aspectos fluidodinâmicos da questão.
Inicialmente, nos primeiros capítulos, os sistemas em que a fase dispersa é diluída são analisados a
partir da fluidodinâmica da partícula isolada; efeitos como aqueles causados pela interação entre partículas são
levados em conta através de modificações do problema inicial.
Para contornar a dificuldade aparentemente intransponível na descrição geométrica do conjunto de
partículas que compõe o sistema denso, os capítulos seguintes utilizam uma Teoria de Misturas com base na
Mecânica do Contínuo. A formulação é estabelecida a partir das leis de conservação aplicadas às fases fluida e
particulada, e mais um conjunto de informações que caracterizam o sistema, as denominadas equações
constitutivas.
A poderosa formulação via Teoria de Misturas, com os seus teoremas, acarreta, no primeiro impacto, o
desconforto causado pela perda do referencial “partícula” na “estrutura amorfa do contínuo”. No cálculo da queda de
pressão no escoamento em duto, problema clássico na Mecânica dos Fluidos, leva-se em conta, por acaso, a estrutura
molecular da matéria? Da mesma forma, na Teoria de Misturas os detalhes da estrutura do Sistema Particulado
escapam pela luneta usada ao revés; as propriedades do sistema são medidas em experiências simples e os resultados
expressos de modo generalizado através das equações constitutivas, tal como na Mecânica dos Fuidos o escoamento
laminar em tubo capilar fornece informações sobre a reologia do fluido.
Não há como negar, o desafio em ministrar por uma centena de vezes a disciplina de Sistemas
Particulados, quer na forma de Operações Unitárias para os estudantes da graduação ou no enfoque de
Fenômenos de Transporte para os pós-graduados, foi sempre a busca de uma teoria que procura amalgamar e
correlacionar os diferentes temas. Assim, por exemplo, o escoamento em meios porosos, a filtração com
formação de torta e o espessamento, guardadas algumas poucas peculiaridades, podem e devem ser tratados
dentro de um mesmo arcabouço; os resultados alcançados na fluidização homogênea levam à reologia da
suspensão e ao projeto das linhas de tranporte hidráulico; a dinâmica da partícula no campo centrífugo permite
analisar o desempenho de ciclcones e de centrífugas.
A cena repete-se anualmente desde 1973, sempre em outubro, na atmosfera acolhedora do anfiteatro
universitário. Entre os veteranos circulam os debutantes tensos. O evento nasceu Encontro sobre o Escoamento
em Meios Porosos (ENEMP) e só recentemente, a partir da 23ª versão, passou a ser Congresso Brasileiro em
Sistemas Particulados. Pois é sobretudo neste foro que os últimos resultados são disseminados entre os grupos
participantes; esta Fluidodinâmica procura respeitosamente preservar e ordenar um pouco da memória dos
Encontros.

Rio de Janeiro, Outubro de 1996

Giulio Massarani

Versão Preliminar da Segunda Edição

A realização desta Versão Preliminar foi concretizada graças ao incentivo e ao apoio desta generosa
população que trabalha no Laboratório de Sistemas Particulados: Christine Lamenha Luna, Cláudia Miriam
Scheid, Flavia Pereira Puget, João Francisco A. Vitor, Marcel Vasconcelos Melo, Marcelo Guilherme G. Mazza,
Marcos Roberto T. Halasz e Sílvia Cristina A. França.

Rio de Janeiro, Julho de 2001

Giulio Massarani

5
 Capítulo 1
Fluidodinâmica da Partícula Sólida

1. Equação do Movimento da Partícula

A fluidodinâmica em sistemas particulados pode ser estudada tomando como ponto de

partida a fluidodinâmica da partícula isolada. A determinação das propriedades do todo pela
extrapolação do comportamento de um elemento da estrutura complexa é intuitiva e didática,
embora, na maioria das situações, esta estratégia exija um grande esforço de imaginação
combinando a um procedimento matemático complicado e duvidoso.

O capítulo 1 procura reunir o conhecimento comum que diz respeito à fluidodinâmica

da partícula, consolidado na literatura a partir do trabalho pioneiro de Stokes sobre a interação
fluido newtoniano-partícula esférica rígida no movimento relativo lento.

C.R. Stokes, "On the Effect of the Internal Friction of Fluids on the Motion of
Pendulums", Trans. Cambridge Phil. Soc., 9,8 (1850).

A fluidodinâmica da partícula pode ser descrita através de um conjunto de equações

que inclui a equação do movimento da partícula, as equações da continuidade e movimento
para o fluido, a condição de aderência na interface fluido-partícula e mais as equações
constitutivas para o fluido e as condições limites pertinentes ao problema específico. A
análise limita-se à fluidodinâmica da partícula rígida, incluindo-se nesta categoria não apenas
as partículas sólidas como também gotas e bolhas de dimensões diminutas. A partícula tem
massa mP , densidade uniforme ρ S , volume V P e a superfície em contato com o fluido é S P .
As equações que seguem são estabelecidas em base a um referencial inercial.

Equação do movimento da partícula

mP (a S ) C = ∫ TF ndS + ρ SV P b . (1)
SP

Equações da continuidade e movimento para o fluido

∂ρ F
+ div(ρ F v F ) = 0 (2)
∂t

 ∂v 
ρ F  F + (grad v F )v F  = div TF + ρ F b. (3)
 ∂t 

Condição de aderência sobre a superfície da partícula

( v F ) Q = ( vs ) C + ω × rQC . (4)

6
Nestas equações, em relação à partícula, ( vs ) C e ( a s ) C são respectivamente a velocidade e a
aceleração de seu centro de massa, ω a velocidade angular e rQC o vetor posição do ponto Q
sobre a superfície da partícula em relação ao centro de massa. Quanto ao fluido, ρF ,v F e TF
são respectivamente a densidade, o campo de velocidades e o tensor tensão que atua sobre
esta fase. b é a intensidade do campo exterior.

A força de interação fluido-partícula pode ser decomposta na força resistiva l e no

empuxo,

∫SP TFndS = l - ρFVPb (5)

sendo nula a força resistiva quando a velocidade relativa entre as fases for nula. A equação
do movimento da partícula toma a forma

m P (as )C = l + (ρS − ρF )VPb (6)

A análise limita-se, deste ponto em diante, ao movimento de translação da

partícula, para atender às necessidades do próximo capítulo sobre a separação sólido-fluido
em sistemas diluidos. Mesmo neste caso relativamente simples, as expressões analíticas
conhecidas para representar a força resistiva restringem-se a algumas configurações
caracterizadas pela forma regular da partícula e pelo movimento relativo partícula-fluido
suficientemente lento, o regime de Stokes, quando a equação do movimento para o fluido,
equação (3), pode ser linearizada.

Os resultados reunidos na tabela (1), alcançados através das equações (1) a (5), são em
maioria exatos ou encerram alguma sorte de aproximação, preservando, no entanto, a forma
analítica do resultado (Berker, 1963). Trata-se de um repertório clássico de soluções que
forma a base para o estudo da fluidodinâmica da partícula.

Os resultados mostram que:

a) A força resistiva exercida pelo fluido sobre a partícula depende das dimensões e
forma da partícula;

b) A força resistiva depende do campo de velocidades do fluido não pertubado pela

presença da partícula;

c) A força resistiva é influenciada pela presença de contornos rígidos e pela presença

de outras partículas;

d) No movimento acelerado da partícula a força resistiva depende da história da

aceleração da partícula.

7
Tabela 1 - Força resistiva fluido-partícula no movimento de translação da partícula no regime de Stokes. O fluido é newtoniano e tem
viscosidade µ. uF é o campo de velocidades do fluido não perturbado pela presença da partícula e vs é a velocidade de translação da partícula
(Berker, 1963).

Descrição uF vS l
Esfera fixa com diâmetro D, (uF ) x = U ∞
escoamento permanente. ( uF ) y = ( uF ) z = 0 vS = 0 l x = 3πµDU∞

Translação retilínea e uniforme de ( vS ) x = v

esfera com diâmetro D, fluido uF = 0 ( vS ) y = ( vS ) z = 0 l x = −3πµDv
inicialmente em repouso
Elipsóide fixo, semi-eixos a, b, c, lx = 3πµD' U ∞
escoamento permanente. 32 πabc
(uF ) x = U ∞ D' =
2 2 2 vS = 0 3 ψ o + a 2α o
x y z ( uF ) y = ( uF ) z = 0
+ + =1 ∞ du ∞ du
a2 b2 c2 ψ o = 2 πabc ∫ , α o = 2 πabc ∫
o ∆u o ( a 2 + u) ∆u

[ ]
1/ 2
∆u = (a 2 + u)(b 2 + u)(c 2 + u)
Esfera fixa com diâmetro D, πD3
escoamento permanente do fluido não uF vS = 0 l = 3πµD(uF )C + (grad P)C ,
8
pertubado pela presença da partícula
onde C denota a posição do centro
resultante do campo de pressões
de massa da partícula
piezométricas P.

8
Tabela 1 (cont.) - Força resistiva fluido-partícula no movimento de translação da partícula no regime de Stokes. O fluido é newtoniano e tem
viscosidade µ. uF é o campo de velocidades do fluido não perturbado pela presença da partícula e vs é a velocidade de translação da
partícula (Berker, 1963).

Descrição uF vs l
Translação retilínea e uniforme da
esfera com diâmetro D em presença de
duas paredes planas paralelas. O fluido
está inicialmente em repouso.
( vS ) x = v  9 1 1 
uF = 0 ( vS ) y = ( vS ) z = 0
lx = −3πµDv 1 + D + 
h1 x h2  32  h1 h 2 

fluido

Translação retilínea e uniforme da

esfera com diâmetro D ao longo do eixo
do tubo com diâmetro Dt. O fluido está
inicialmente em repouso.

( vS ) x = v  D
uF = 0 lx = −3πµDv1 + 2,1 
( vS ) y = ( vS ) z = 0  Dt 

Dt

9
Tabela 1 (cont.) - Força resitiva fluido-partícula no movimento de translação da partícula no regime de Stokes. O fluido é newtoniano e tem
viscosidade µ. uF é o campo de velocidades do fluido não perturbado pela presença da partícula e vS é a velocidade de translação da
partícula (Berken, 1963).

Descrição uF vs l
Translação retilínea e uniforme das
esferas 1 e 2 com diâmetro D1 e D2. O
fluido está inicialmente em repouso.
f2   3 D2 
(vS ) x1 = (vS ) x2 = v  f1 = 3πµD1v  1 − 8 h 
uF = 0 
(vS ) y1 = (vS ) y2 = 0  f 2 = 3πµD2 v  1 − 3 D1 
f1   8 h
(vS ) z1 = (vS ) z2 = 0
h φ
q2 9 cos φ
q1 = q2 = πµD1D2 v
v 8 h

q1
v

Esfera em translação retilínea não 1 dv

uF = 0
-lx == πD3ρF + 3πµDv
uniforme e com velocidade inicial nula. (vS ) x = v (t ), v (0) = 0 12 dt
O fluido está inicialmente em repouso. (vS ) y = ( vS ) z = 0 dv
3 t
+ D 2 (πµρ)1 / 2 ∫ dτ dτ
2 o t−τ

10
Na situação em que a partícula apresenta forma irregular e fora do regime de Stokes, não
parece haver outra alternativa senão a de tratar a força resistiva de modo empiríco,
procurando generalizar os resultados clássicos (Bird et al., 1960, p.193):

u−v
l = A ⋅ 1 ρF u − v 2 ⋅ c D ⋅ , (7)
2 u−v

onde A é uma área característica, cD o coeficiente de arraste cujo valor numérico depende da
definição de A, u é a velocidade do fluido não perturbado pela presença da partícula na
posição do centro de massa desta partícula, e v a velocidade de translação da partícula.
Considera-se na equação (7) que a força resistiva e a velocidade relativa

U = u−v (8)

tenham a mesma direção, o que implica em admitir que a forma da partícula apresenta um
certo grau de regularidade. Nestas condições, a equação do movimento da partícula toma a
forma

1
mP a S = Aρ F cD U U + + ( ρ S − ρ F )V P b . (9)
2

2. A Força Resistiva Fluido-Partícula

O estudo da fluidodinâmica da partícula requer o conhecimento da reologia do fluido e

das propriedades físicas da partícula expressas pela densidade, dimensão e forma. Entre as
múltiplas possibilidades conhecidas na caracterização da partícula e para melhor usufruir um
grande número de dados experimentais disponíveis na literatura, adotam-se neste texto o
diâmetro volumétrico DP como dimensão característica e a esfericidade φ na
caracterização da forma da partícula (Allen, 1981).

O diâmetro volumétrico é definido como sendo o diâmetro da esfera com o mesmo

volume que a partícula,

1/ 3
6 
DP =  V P  . (10)
π 

O valor desta propriedade para partículas de forma irregular pode ser determinado com o
auxílio da picnometria clássica ou, na situação em que as partículas são diminutas, através da
análise granulométrica realizada no Coulter Counter (Allen, 1981).

A esfericidade é definida como sendo o cociente entre a superfície da esfera com o

mesmo volume que a partícula e a superfície S P ,

φ = πDP2 / S P . (11)

11
 A esfericidade é um fator de forma empírico que pode ser determinado por
permeametria, técnica que será apresentada em detalhes no capítulo 3. É a partícula esférica
que apresenta o maior valor da esfericidade, φ=1; as partículas que ocorrem usualmente,
como aquelas resultantes dos processos de moagem, apresentam a esfericidade na faixa de 0,5
a 0,7.

O coeficiente de arraste cD , presente na equação que define a força resistiva fluido-

partícula, equação (7), pode ser calculado através da medida da velocidade terminal da
partícula vt , isto é, a velocidade constante atingida pela partícula quando lançada no fluido
inicialmente em repouso. Definindo a área característica desta equação como sendo a área da
seção transversal da esfera de diâmetro DP ,

A = πDP2 / 4 (12)

resulta no campo gravitacional, a partir das equações (8) e (9).

U z = 0 − vt = − v t (13)

4 (ρ S − ρ F ) DP g
cD = .
3 ρ F vt2 (14)
vt
Um grande número de experiências conduzidas com partículas isométricas, isto é, partículas
esféricas ou na forma de poliedros regulares (tetraedro, cubo, octaedro, icosaedro e
dodecaedro), parecem indicar que o valor do coeficiente de arraste depende apenas do número
de Reynolds,

D P vt ρ F
Re = (15)
µ

e da esfericidade (Pettyjohn e Christiansen, 1948). Generalizande este resultado,

4 (ρ S − ρ F ) DP b
cD = = f1 ( Re, φ) (16)
3 ρ FU 2

DPUρ F
Re = (17)
µ

b = b , U = U = u−v . (18)

A partir da equação (16):

Re = f 2 (cD Re 2 , φ) (19)

Re = f 3 (cD / Re, φ) (20)

12
onde os grupos adimensionais cD Re 2 e cD / Re são assim calculados

4 ρ F (ρS − ρ F )bDP3
cD Re2 = (21)
3 µ2

4 (ρ S − ρ F )µb
cD / Re = . (22)
3 ρ2F U 3

Cabe ressaltar que a correlação expressa pela equação (16) é o ponto de partida para o
estabelecimento das equações (19) e (20) e que pode ser utilizada com vantagem no estudo da
dinâmica da partícula em fluido não newtoniano pelo fato da viscosidade estar presente
apenas no número de Reynolds. A equação (19) presta-se para o cálculo de U , pois cD Re 2
não inclui esta variável; analogamente, a equação (20) deve ser utilizada no cálculo de DP já
que c D /Re não inclui esta variável. Nestas duas últimas situações, U e DP são calculados a
partir do número de Reynolds.

As correlações apresentadas nas tabelas (2) a (4) referem-se à fluidodinâmica da

partícula isométrica isolada em fluido newtoniano. Embora a tabela (3) inclua a partícula
esférica, recomenda-se neste caso, para maior precisão, a utilização da tabela (2). A tabela (4)
fornece diretamente as expressões para a velocidade relativa fluido-partícula e para o
diâmetro da partícula quando prevalece o regime de Stokes ou o de Newton , isto é, quando
Re < 0,5 ou 103 < Re < 2 × 105 . As correlações das tabelas (2) e (3) foram estabelecidas
através do Método das Duas Assíntotas de Churchill (1983).

y ( x ) = [ yon ( x ) + y∞n ( x )]1/ n , (23)

onde yo ( x ) e y∞ ( x ) referem-se, respectivamente, aos regimes de Stokes e Newton, e o

“valor ótimo” de n é determinado a partir de dados experimentais, dentro de algum critério
estatístico.

Entre outras correlações apresentadas na literatura para a fluidodinâmica da partícula

isométrica, cabe mencionar as de Concha e Barrientos (1986) e Haider e Levenspiel (1989).
Estas correlações, baseadas essencialmente nos dados experimentais de Pettyjohn e
Christiansen (1948), são de complexidade e precisão equivalentes àquelas apresentadas na
tabela (3).

Em algumas situações foram levantadas correlações específicas para descrever a

fluidodinâmica da partícula não-isométrica (Concha e Christiansen, 1986), porém, na falta
destas, utilizam-se os resultados relativos à partícula isométrica, caracterizando a forma da
partícula não-isométrica através da esfericidade.

13
Tabela 2 - Fluidodinâmica da partícula esférica isolada:
Correlações de Coelho & Massarani (1996) com base nos dados de Lapple & Shepherd
(1940) e Pettyjohn & Christiansen (1948).

Re < 5 × 104

Correlação n Valor Médio e Desvio Padrão

1/ n ( cD ) exp
 24  n  = 1, 00 ± 0, 09
c D =   + 0,43n  0,63 ( cD ) cor
 Re  

−1/ n ( Re) exp

 2  −n  2  −n/2  = 1,00 ± 0,06
c Re c Re 0,95
Re =  D  + D   ( Re) cor
 24   0,43  
 

 24  n / 2  0,43  n 
1/ n ( Re) exp
0,88 = 1,00 ± 0,09
Re =   +   ( Re) cor
 cD / Re   cD / Re  
 

DPUρF 4 ρF (ρS − ρF )bDP3 4 (ρS − ρF )µb

Re = , cD Re 2 = 2
, cD / Re =
µ 3 µ 3 ρF2 U 3

14
Tabela 3 - Fluidodinâmica da partícula isométrica isolada:
Correlações de Coelho & Massarani (1996) com base nos dados de Pettyjohn & Christiansen
(1948).

0,65 < φ ≤ 1 e Re < 5 × 104

Correlação n Valor Médio e Desvio Padrão

 24  n 
1/ n ( cD ) exp
n 0,85 = 1, 00 ± 0, 13
cD =   + K2 ( cD ) cor
 K1Re  
 

 −n −n/2  −1/ n ( Re ) exp

K c Re2   c Re 2 
1,2 = 1, 00 ± 0, 10
Re =  1 D  + D   ( Re ) cor
 24   K2  
 

1/ n (Re)exp
 
n/2
 K2  
n
Re = 
24 = 1,00 ± 0,14
 +   1,3 (Re)cor
 K1 (cD / Re)   cD / Re  
 

DPUρF 4 ρF ( ρS − ρF )bDP3 4 ( ρS − ρF )µb

Re = , cD Re 2 = , c D / Re =
µ 3 µ2 3 ρF2 U 3

K1 = 0, 843 log10 ( φ / 0, 065), K2 = 5, 31 − 4 , 88φ

15
Tabela 4 - Fluidodinâmica da partícula isométrica isolada:
Cálculo da velocidade e do diâmetro da partícula (Pettyjohn & Christiansen, (1948).

0,65 < φ ≤ 1

Variável a Ser Estimada Regime de Stokes Regime de Newton

Re < 0, 5 103 < Re < 5 × 104

cD 24 K2
K1 Re

1/ 2
( ρS − ρF )bK1 D 2p  4( ρS − ρF )bD p 
U
18µ  
 3ρF K2 

1/ 2
 18µU  3ρF K2U 2
Dp   4( ρS − ρF )b
 ( ρS − ρF )bK1 

K1 = 0, 843 log10 ( φ / 0, 065), K2 = 5, 31 − 4 , 88φ

16
 Exemplo

Deseja-se estudar a possibilidade de separar o minério A do minério B através da

elutriação com corrente ascendente de água.

A+B
A
Propriedades do minério A: ρSA = 2, 2 g / cm3 e φA = 0, 70

Propriedades do minério B: ρSB = 3,2 g / cm3 e φB = 0,85

Água Faixa granulométrica da mistura A+B: 0,149 < DP < 0,595mm ,

correspondendo às peneiras 28/100 # Tyler.
A+B

A velocidade de elutriação de água (20ºC) que permite recuperar a maior quantidade

possível do produto A puro é igual à velocidade terminal da menor partícula de B, isto é,
DPB = 0,149mm . Resulta da tabela 3, utilizando as propriedades de B: cD Re 2 = 95,2 ,
Re = 2 ,93 e, deste último, uF = vtB = 1,97cm / s .

Conhecida a velocidade de elutriação, é possível calcular o diâmetro da maior

partícula de A presente no produto arrastado. A tabela (3) leva aos seguintes resultados
utilizando as propriedades de A: cD / Re = 2,05 , Re = 3,97 e, deste último, DPA = 0,202mm .

Em conclusão, a velocidade de elutriação u = 1,97cm / s leva a um produto de topo

constituido de A puro na faixa granulométrica 0,149 < DP < 0,202mm ; o produto de fundo é
constituido de uma mistura de A (0,202 < DP < 0,595mm) e de B (0,149 < DP < 0,595mm) .
Cabe salientar que esta análise trata apenas das condições de separabilidade dos componentes
A e B na elutriação e nada informa sobre a cinética de separação. Pode-se esperar que as
partículas maiores de A sejam arrastadas muito lentamente e que as partículas menores de B
sedimentem também muito lentamente.

Efeito da presença de fronteiras rígidas

Resultados analíticos reunidos na tabela (1) evidenciam que a fluidodinâmica da

partícula é influenciado pela presença de fronteiras rígidas, resultando uma redução na
velocidade terminal em relação à velocidade terminal da partícula isolada, v∞ .

17
 Almeida (1995) estudou experimentalmente o
movimento da partícula isométrica ao longo do eixo
principal de um tubo cilíndrico com diâmetro Dt ,
resultando a figura (1) e as correlações empíricas
apresentadas na tabela (5). Cabe ressaltar que as
correlações clássicas de Francis (1933), regime de Stokes,
e de Munroe (1888), regime de Newton, válidas para
esferas, podem ser utilizadas também para partículas
vt isométricas.

Dt

0
vt 1
0,9 0,05
v∞ 0,1
0,8
0,7
0,2
0,6
0,5 0,3
0,4
0,4
0,3
Dp Francis (1933)
β= = 0,5
0,2 Dt
Almeida (1995)
Munroe (1888)

-2
10
-2 -1 2 3
10 10 1 10 10 10
Dpv∞ρF
Re∞ =
µ

Figura 1 - Efeito de parede na velocidade terminal da partícula isométrica (Almeida, 1995).

18
Tabela 5 - Efeito de parede na fluidodinâmica da partícula isométrica em fluido newtoniano
(Almeida, 1995): 0,65 < φ ≤ 1 e 0 < D P / D t ≤ 0,5 .

D P v∞ ρ F vt
Re∞ = kP = , β = DP / Dt
µ v∞

4
< 0,1  1− β 
kP =  
(Francis, 1933) 1 − 0,475β 

10
kP =
1 + A Re∞B
3
0, 1 − 10
A = 8,91e 2 ,79β , B = 117
, × 10−3 − 0,281β

> 103 k P = 1 − β 3/ 2
(Munroe, 1888)

24 e 3,54β
Re = n
, n = 0,85 para Re < 35
K1 (cD − K2n )1/ n

4 (ρS − ρF )gD P
cD =
3 ρF v 2t

φ
K1 = 0,843 log10 , K2 = 5,31 − 4 ,88φ
0,065

19
 Exemplo

Deseja-se planejar uma experiência que consiste na medida da velocidade terminal

limitando, com a escolha adequada do diâmetro do cilindro de testes, o efeito de parede a 5%,
v
isto é, k P = t > 0,95 . A partícula tem diâmetro DP = 5mm . Utilizando as correlações de
v∞
Francis e Munroe, tabela (5),

Regime de Stokes : Dt / DP > 41, Dt > 205mm ;

Regime de Newton: Dt / DP > 8, Dt > 40mm .

Os resultados evidenciam que o efeito da parede e bem mais agudo no regime de

Stokes que no regime de Newton.

Influência da concentração de partículas

Um grande número de dados experimentais apresentados na literatura evidencia que a

velocidade terminal de uma partícula tem seu valor substancialmente reduzido pela presença
de outras partículas. Esta redução, tanto mais sensível quanto maior a concentração de
sólidos, é da ordem de 5% para concentrações de apenas 2%, como mostra a equação de
Einstein (Govier e Aziz, 1972, p. 98).

vt / v∞ = 1 / (1 + 2 ,5cV ) , (24)

onde v∞ é a velocidade terminal da partícula isolada e cV a fração volumétrica da fase sólida

na suspensão.

O efeito da presença da fase particulada na fluidodinâmica de suspensões é

comumente expresso através de correlação do tipo (Richardson e Zaki, 1954).

U / v∞ = f ( Re∞ , ε ) , (25)

onde U é o módulo da velocidade relativa fluido-partícula,

U = v−u ,

Re∞ o número de Reynolds referente à velocidade terminal da partícula isolada,

D P v∞ ρ F
Re∞ = ,
µ

ε, a porosidade, é a fração volumétrica de fluido na suspensão,

ε = 1 − cV .

20
 As correlações referentes à equação (25) podem ser determinadas através da
experimentação conduzida na sedimentação em batelada e na fluidização homogênea: no
primeiro caso U = v / ε , onde v é a velocidade da frente de sedimentação; no segundo caso
U = QF / (εA) , sendo QF a vazão de fluido e a A área da seção transversal de fluidização
(Barnea e Mizrahi, 1973). A experimentação torna-se imprecisa quando a faixa
granulométrica das partículas sólidas é extensa e quando a concentração de sólidos é
reduzida, inferior a 5% em volume, resultando nas duas situações uma interface fluido-
suspensão pouco nítida por problemas de segregação de partículas.

A maioria das correlações apresentadas na literatura referem-se a amostras com

partículas "arredondadas", em faixa granulométrica "estreita" representada por um diâmetro
médio que possivelmente não caracteriza a fluidodinâmica da suspensão. Como conseqüência
da caracterização incompleta do sistema particulado, as correlações da literatura podem
diferir substancialmente entre si. São apresentadas na tabela (6) as correlações de Richardson
e Zaki (1954) para partículas arredondadas, a de Politis e Massarani (1989) para partículas
irregulares e outras resultantes dos dados experimentais reunidos por Concha e Almendra
(1978). Na figura (2) é feita a comparação entre os resultados de Richardson e Zaki (1954) e
Almendra (1979) para partículas arredondadas: as maiores discrepâncias ocorrem quando a
porosidade é elevada e na região intermediária entre os regimes de Stokes e Newton.

21
Tabela 6 - Influência da concentração de partículas na fluidodinâmica de suspensões.

A. Correlação de Richardson e Zaki (1954) para partículas arredondadas:

U / v∞ = ε n , n = n( Re∞ )

Re∞ < 0,2 0,2-1 1-500 > 500

n 3,65 4 ,35 Re∞−0 ,03 − 1 4 ,45 Re∞−0 ,1 − 1 1,39

B. Correlação de Politis e Massarani (1989) para partículas irregulares (areia, hematita,

itabirito, dolomita e quartzo, 0,47<φ<0,80).

−0 ,14
U / v∞ = ε 5,93 Re∞ , 9,5 < Re∞ < 700 .

O diâmetro médio é a média aritmética da abertura das peneiras de corte.

C. Correlações empíricas estabelecidas com base nos dados experimentais reunidos por
Concha e Almendra (1978) (Massarani e Santana, 1994)

U  0,83ε 3,94 , 0,5 < ε ≤ 0,9

Re ∞ < 0,2 , =
v∞ 4,8ε − 3,8, 0,9 < ε < 1

U 1
1 < Re ∞ < 500, = −B
, 0,5 < ε < 0,95
v∞ 1 + A Re ∞
A = 0,28ε −5,96 , B = 0,35 − 0,33ε
U
Re ∞ > 2 × 103 , = 0,095exp(2,29ε ), 0,5 < ε < 0,95.
v∞

22
 1
U/v∞
ε = 0,95
0,8
0,95

0,90
0,6 0,90

0,80
0,4
0,80
0,70

0,2 0,70 Richardson e Zaki

0,60 Almendra
0,60
0
-2 -1 2 3 4
10 10 1 10 10 10 10
Re∞

Figura 2 - Influência da concentração de partículas na fluidodinâmica de suspensões:

comparação entre os resultados de Richardson e Zaki (1954) e Almendra (1979).

23
 Outra estratégia que pode ser adotada na análise de fluidodinâmica de suspensões
consiste em considerar o comportamento isolado de uma partícula no seio da mistura sólido-
fluido, mistura esta caracterizada pela densidade e viscosidade ρSusp e µSusp(Govier e Aziz,
1972, p.98; Massarani e Santana, 1994). Assim, no regime de Stokes, tabela (4),

(ρS − ρSusp )gK1D 2P

v= = εU . (26)
18µSusp

Sendo

( ρ S − ρ F )gK1 DP2
v∞ = (27)
18µ

ρ S − ρ Susp = ε (ρ S − ρ F ) , (28)

resulta, combinando as equações (25) a (28),

µ
µ Susp = = µ / f ( Re∞ ,ε) . (29)
U / v∞

Finalmente, cabe indagar em que medida podem estar relacionados entre si os

resultados clássicos da fluidodinâmica nos meios de densos, estabelecidos no contexto da
Teoria de Misturas, e os da fluidodinâmica de suspensões estabelecidos a partir do
comportamento da partícula isolada. O assunto será abordado nos capítulos 3 e 4. Demonstra-
se, por exemplo, que no regime de Stokes

1 ( DP φ ) 2 ε 2
U= ⋅ ⋅ ⋅ (ρ S − ρ F ) g (30)
µ 36β 1 − ε

ou, de modo equivalente,

4 36β 1 − ε 1
cD = ⋅ ⋅ ⋅ , (31)
3 φ 2 ε 2 Re

sendo

β = 3,82 / (1 − ε 7,21 ) e ε < 0,97 .

24
 Exemplo

Deseja-se calcular a porosidade no transporte vertical ascendente, em duto com

diâmetro Dt = 5,1cm , de partículas sólidas com as seguintes propriedades: diâmetro
DP = 1mm , densidade ρ S = 3g / cm3 e esfericidade φ = 0,75 .

a) O fluido é água (ρ F = 1g / cm3 e µ = 0,9cP) e as vazões de fluido e sólido são

respectivamente QF = 15m3 / h e QS = 3m3 / h .

b) O fluido é ar a 20ºC e 1 atm (ρ F = 1,2 × 10−3 g / cm3 e µ = 1,8 × 10-2 cP) e as

vazões de fluido e sólido são respectivamente QF = 39,9m3 / h e QS = 1,32m3 / h .

A porosidade no transporte vertical pode ser calculada resolvendo a equação (25),

QF QS
U= − = v∞ f ( Re∞ , ε ) ,
εA (1 − ε ) A

onde A = 20,4cm2 é a área da seção transversal do duto. Uma estimativa do valor da

porosidade pode ser alcançada a partir do conhecimento das vazões de cada fase,

QF uεA ε
α= = = . (32)
QF + QS uεA + v (1 − ε ) A ε + (1 − ε ) v
u

v
Quando tende a 1, ε tende a α.
u

Na solução deste exemplo admite-se que as correlações de Richardson e Zaki (1954), tabela
(6), sejam válidas apesar das partículas não serem arredondadas.

2
c D Re∞ Re∞ n v∞ α ε u v u
(eq. 21) (tab. 3) (tab. 6) (cm/s) (eq. 32) (eq. 25) (cm/s) (cm/s) v
Trans.
Hidráulico 3,23x104 134 1,73 12,1 0,833 0,829 246 239 1,03

Trans.
Pneumático 1,45x105 295 1,52 443 0,968 0,921 592 228 2,60

25
 O fato da densidade e viscosidade da água serem muito maiores do que estas
propriedades para o ar explica os resultados esperados de que a velocidade de deslizamento
u − v é muito menor no primeiro caso do que no segundo.

3. O Movimento Acelerado da Partícula

O movimento retilínio acelerado de uma esfera no seio de um fluido newtoniano,

regime de Stokes, foi estudado no final do século passado por Basset (Berker, 1963, p.241).
No caso da queda livre da partícula partindo do repouso em fluido inicialmente estagnado a
força resistiva toma a forma indicada na tabela (1):

dv
ρ dv 3 t
l = F VP + 3πDµv( t ) + D 2 (πµρF )1 / 2 ∫ dτ dτ . (33)
2 dt 2 0 t−τ

O primeiro termo do segundo membro da equação fornece o valor da força resistiva que o
fluido ideal em escoamento potencial exerce sobre a partícula; o segundo termo exprime o
resultado clássico de Stokes para o movimento retilíneo e uniforme de uma esfera em fluido
viscoso; o terceiro termo evidencia a ação “hereditária” do fluido sobre a partícula, pois
explicita o fato de que a força resistiva depende da história da aceleração da partícula.

Substituindo a equação (33) na equação do movimento da partícula, equação (6),

resulta a equação integro-diferencial

dv
 ρ  dv 3 t dτ dτ ,
 ρ S + F V P = (ρ S − ρ F )V P g − 3πDµv (t ) − D 2 ( πµρ F )1/ 2 ∫ (34)
 2  dt 2 0 t−τ

que pode ser resolvida analiticamente por diferentes técnicas (Clift et al., 1979, p. 285;
Hackenberg, 1991).

O resultado expresso pela equação (34) mostra que a aceleração inicial da partícula

2 (ρ S − ρ F )
a (0) = g, (35)
2ρ S + ρ F

tende ao valor da intensidade do campo gravitacional g quando ρS >> ρF e que esta

aceleração é nula no caso limite em que as densidades do fluido e da partícula forem iguais
entre si. Desprezando o efeito da história da aceleração da partícula, a integração da
equação (34) fornece

vt  36µ 
= exp  2
t , (36)
vt − v  (2ρ S + ρ F ) D 

onde vt é a velocidade terminal da partícula

26
 (ρ S − ρ F ) gD2
vt = .
18µ

Exemplo

O diâmetro da esfera sólida de densidade ρ S = 3g / cm3 em queda livre no ar

(ρ S = 1,2 × 10−3 g / cm3 e µ = 1,8 × 10−4 P) e na água (ρ F = 1g / cm3 , µ = 10−2 P ) , no limite
de validade do regime de Stokes,

Dvt ρ F
Re∞ = = 0,5 ,
µ

é de respectivamente 43µm e 77µm. Para os sistemas assim definidos, a integração da

equação (34) e a equação (36) conduzem à figura (3) (Hackenberg, 1991). Pode-se observar
que o regime permanente é atingido em fração de segundo, sendo que a água leva a uma
resposta mais rápida inicialmente e mais retardada ao final. O efeito da história da aceleração
da partícula é importante no caso da água e desprezível no caso do ar.

1
v/vt
água
0,8

ar
0,6

Eq. (34)
0,4
ar
Eq. (36)
água

0,2

0
0 0,02 0,04 0,06 0,08 0,10 0,12 0,14 t (s)
3
Figura 3 - Movimento acelerado das esferas com densidade 3g / cm e diâmetro D = 43µm
em ar e D = 77µm em água (Re ∞ = 0,5) .

Face às evidentes dificuldades tanto na abordagem teórica quanto experimental, a

literatura evidencia uma grande carência de informações relativas ao movimento acelerado da

27
partícula fora do regime de Stokes e quando estas não são esféricas (Marchildon e Gauvin,
1979; Renganathan et al., 1989). Na situação em que ρS >> ρF, Renganathan et al. (1989), em
abordagem empírica, consideram a queda livre da partícula como descrita pela equação do
movimento

dv πD 2 ρ v2
mP = VPρS g − cD F (37)
dt 4 2

em que o coeficiente de arraste cD = f (Re) preserva a forma funcional das correlações

alcançadas no movimento estabelecido, cD∞ = f ( Re∞ ) .

4. Dinâmica da Partícula em Fluido Não Newtoniano

Os estudos teóricos relativos ao escoamento de fluidos não newtonianos nas

vizinhanças de esferas rígidas restringem-se aos casos em que prevalece o regime de Stokes.
Neste sentido, cabe mencionar os trabalhos de Caswell (1962, 1970).

A estratégia usada neste capítulo e nos seguintes consiste em estender a formulação

clássica sobre a dinâmica da partícula sólida em fluidos newtonianos para contemplar também
uma classe ampla de fluidos não newtonianos: o elo de ligação é a viscosidade efetiva µef que
pode ser calculada através da tensão cisalhante S, uma propriedade material do fluido, e da
taxa de deformação característica λ∗ , uma propriedade cinemática de escoamento (Massarani
e Silva Telles, 1978),

µ ef = S ( λ* ) / λ* (38)

A taxa de distensão característica λ∗ pode ser determinada empiricamente através da

medida experimental da velocidade da partícula com o auxílio, por exemplo, das relações
cD × Re apresentadas nas tabelas (3), (5) e (6) na seguinte seqüência:

λ* a partir
4 (ρ S − ρ F ) D P b D vρ
vt → c D = → Re → µ ef = P t F → de .
3 ρ F vt2 Re * *
µ ef = S ( λ ) / λ

A propriedade cinemática λ* pode ser representada do modo,

v
λ* = α (39)
D

onde v* e D* são respectivamente uma velocidade e dimensão características e α um fator de

configuração adimensional. Estão reunidos na tabela (7) os resultados obtidos para a dinâmica
da partícula isolada e nos casos em que são levados em conta os efeitos de parede e de
concentração.

28
 Tabela 7 - Dinâmica da partícula sólida em fluido não-newtoniano

Descrição Taxa de Distensão λ∗ Referência

Partículas esféricas e não v Laruccia

0,39(−8,85φ2 + 13,46φ − 3,62)
esféricas isoladas DP (1990)

Deslocamento da partícula v Almeida

esférica ao longo do eixo 0,39e6,81β ⋅ (1995)
D
principal do tubo (β = D / D t < 0,5)

Efeito de concentração na 1− ε U Silva Telles e

fluidodinâmica de 9 ⋅ Massarani (1979)
φε D P
partículas
(ε < 0,9)

Exemplo

Reômetro de Stokes para fluidos não newtonianos

Deseja-se determinar a relação S = S ( λ ) para um fluido

não newtoniano através da medida da velocidade de
deslocamento de esferas neste fluido. Diâmetro do tubo,
Dt = 20mm . Densidade do fluido, ρ F = 1,15g / cm3 .
Dados:

Exp. ρS D v
β
nº (g/cm³) (cm) (cm/s)
1 2,55 0,20 0,10 0,72
2 2,55 0,50 0,25 3,61
3 3,98 0,30 0,15 3,78
v
4 3,98 0,50 0,25 8,87
Dt
5 7,60 0,30 0,15 9,85
6 7,60 0,50 0,25 22,3

29
Resulta:

cD Re D P v ρF
λ* ( s −1 ) µef = S(λ* ) = µ ef λ*
Exp. nº (eq. 14) (tab.5) Re
(tab. 7) (dyn / cm 2 )
(P)

1 2,77 614 0,056 2,96 8,20

2 15,1 61,1 0,957 2,17 32,8
3 13,6 67,6 0,606 2,15 29,2
4 38,0 20,5 2,94 1,73 65,7
5 35,6 22,7 1,86 1,83 65,2
6 95,5 7,38 8,75 1,47 140

O resultado pode ser expresso de modo conveniente através de

0,80 2 −1
S = 3,65λ dyn / cm , válido para 2 < λ < 100s .

30
 Problemas: Fluidodinâmica da Partícula Sólida

1. Foram os seguintes os resultados obtidos na elutriação de 25 g de um pó industrial com

água a 30°C, numa vazão de 37 cm3/min:

Elutriador Diâmetro do Massa

tubo (cm) recolhida (g)
1 3,0 4,62
2 4,0 6,75
3 6,0 7,75
4 12,0 4,42

Determinar a distribuição granulométrica da amostra em termos do diâmetro de Stokes,

sabendo-se que a densidade do sólido é 1,8 g/cm3.

Resposta:

Elutriador Diâmetro (cm) Velocidade do fluido DP X

(cm/s) (µm)

1 3,0 8,72×10-2 44,9 0,815

2 4,0 4,91×10-2 33,7 0,545

3 6,0 2,18×10-2 22,4 0,235

4 12,0 5,45×10-3 11,2 0,058

31
2. Calcular a velocidade de sedimentação de uma suspensão de partículas em querosene.

Propriedades do fluido: densidade 0,9 g/cm3 e viscosidade 2,3 cP.

Propriedades das partículas: densidade 2,3 g/cm3, diâmetro médio 0,8 mm, esfericidade 0,8.
Concentração de sólidos na suspensão: 260 g/l de suspensão.

Resposta:

Porosidade da suspensão: 0,887.

Velocidade terminal da partícula isolada: 7,71 cm/s
Velocidade de sedimentação da suspensão: 5,25 cm/s.

3. Os seguintes dados foram obtidos em ensaios de sedimentação de partículas de Al2O3 em

água, a 25°C:

c( gAl2 O3 / cm3 de suspensão) 0,041 0,088 0,143 0,275 0,435

v (cm / min) 40,5 38,2 33,3 24,4 14,7

A densidade das partículas é 4,0 g/cm3 e a esfericidade é estimada em 0,7.

a) Determinar, pela extrapolação dos dados, a velocidade terminal das partículas à diluição
infinita e, a partir deste valor, calcular Dp (diâmetro da esfera de igual volume que a
partícula);
b) Comparar os resultados experimentais com as estimativas segundo a correlação empírica
de Richardson & Zaki.

Resposta:

Dados experimentais:
v = −223 + 267 ε cm / min( R 2 = 0,996) .
Velocidade terminal calculada por extrapolação dos dados experimentais: 43,3 cm/min.
Diâmetro volumétrico das partículas: 72 µm.

ε 0,990 0,978 0,964 0,931 0,891

v (cm / min) 40,5 38,2 33,3 24,4 14,7
v = v∞ ε 4,43 (cm / min) 41,4 39,2 36,8 31,5 25,9

4. Michael e Bolger (IEC Fundam., 1, 24, 1962) desenvolveram um método que permite a
caracterização de partículas floculadas (diâmetro e densidade médios, grau de floculação e
velocidade de sedimentação dos flocos). Uma vez determinada experimentalmente a
velocidade de sedimentação da suspensão v a diferentes concentrações co, os parâmetros
desejados podem ser estimados através do seguinte sistema de equações:

v = v∞ (1 − kco ) 4 ,65 (Correlação de Richardson e Zaki)

32
 D 2fl (ρ fl − ρ F ) g
v∞ = (Equação de Stokes)
18µ

ρS − ρ F
ρ fl − ρ F = (Balanço de massa),
kρ S

onde

v - velocidade de sedimentação da interface lodo-líquido clarificado no ensaio em

batelada;
v∞ - velocidade terminal do floco à diluição infinita;
k - volume de flocos por unidade de massa de sólido seco (fornece o grau de floculação);
co - concentração em massa de sólido seco por unidade de volume de supensão;
D fl - diâmetro médio dos flocos;

ρ fl - densidade média dos flocos;

ρ F - densidade do fluido;
ρ S - densidade do sólido seco;
g - aceleração da gravidade;
µ - viscosidade do fluido.

Calcular as propriedades caracaterísticas (v∞ , D fl , k e ρ fl ) dos flocos de hidróxido de cálcio

de uma suspensão aquosa (agente de floculação: alúmen) sabendo-se que a 25ºC:

co ( g / cm3 ) 6×10-3 8×10-3 10×10-3 12,5×10-3 15×10-3 20×10-3 25×10-3 30×10-3

v (cm / min) 4,77 4,32 3,65 3,04 2,33 2,08 1,37 0,30.

A densidade do sólido seco é 2,20 g/cm3.

Resposta:

v∞ = 7,89 cm/min e k = 14,7 cm3/g (1ª equação).

ρ fl = 1,037 g/cm3 (3ª equação).
D fl = 256 µm (2ª equação).

5. Determinar as respectivas velocidades de elutriação para separar pó de diamante nas faixas

0-1 µm, 1-2 µm, 2-3 µm (diâmetro da esfera de igual volume que a partícula). A densidade

33
do diamante é 3,5 g/cm3 e a esfericidade das partículas 0,7. O fluido de arraste é água a 20ºC.
(P.Grodzinski, “Diamond Technology”, NAG Press Ltd., Londres, 2ª edição, p. 349, 1953).

Resposta:

Faixa granulométrica (µm) 0-1 1-2 2-3

Velocidade de elutriação (cm/h) 0,427 1,71 3,84.

6. Uma mistura finamente dividida de galena e calcário na proporção 1:4 em massa é sujeita à
elutriação com corrente ascendente de água com velocidade de 0,5 cm/s. A distribuição
granulométrica dos dois materiais é a mesma:

Dp (µm) 20 30 40 50 60 70 80 100
100X 15 28 43 54 64 72 78 88.

Calcular a percentagem de galena no material arrastado e no produto de fundo.

Galena: densidade 7,5 g/cm3 e esfericidade das partículas 0,8.
Calcário: densidade 2,7 g/cm3 e esfericidade das partículas 0,7.
Temperatura da água: 20ºC.

Resposta:

1
Análise granulométrica da alimentação X = , Dp em µm.
2, 27
 44,6 
1+  
 Dp 
 

4 (ρ S − ρ F )µg Diâmetro % Massa

cD / Re =
Material 3 ρ2F u 3 Re Crítico, µm Arrastada
Calcário 178 0,395 73,8 0,76
Galena 680 0,196 39,2 0,43

% galena na alimentação : 20,0

% galena no produto de fundo: 37,3
% galena no produto de topo : 12,4.

7. O separador de poeira opera em 3 compartimentos, como mostra o esquema abaixo

representado. Estimar a faixa granulométrica das partículas retidas em cada compartimento
sabendo-se que a vazão de gás (ar a 20ºC e 1 atm) é 140 m3/min, a densidade das partículas é
3 g/cm3 e sua esfericidade 0,75.

34
Resposta:

4 (ρ S − ρ F )µg Faixa
Compartimento L(m) vt = H u / L cD / Re = 3 Re DP
ρ2F vt3 Granulométrica
( m / s) (µm)
(µm)
1 1,5 0,390 8,27 1,84 72,2 >72,2
2 3 0,195 66,2 0,640 49,2 49,2-72,2
3 4,5 0,130 223 0,347 40,0 40,0-49,2

8. Dimensionar um rotâmetro tronco de cone-esfera para medir a vazão de água (20ºC) na

faixa de 1 a 3 m3/h. O flutuador é uma esfera de aço com 1 cm de diâmetro e densidade 7,7
g/cm3. Que faixa de vazões este mesmo rotâmetro mediria se o fluido fosse ar a 20ºC e 1
atm?

Resposta:

A altura h não influencia o desempenho do rotâmetro. Pode ser da ordem de 20 cm por

questão de comodidade e precisão na leitura da escala do aparelho.
Faixa de vazão de ar: 30,4 a 97,8 m3/h.

35
 Bibliografia

Allen, T., “Particle Size Measurement”, Chapman & Hall, Londres, 3¦ edição, 454 p. (1981).

Almendra, E.R., “Velocidade de Sedimentação de Sistemas Particulados”, Tese de M.Sc.,

Programa de Eng. Metalúrgica e de Materiais, COPPE/UFRJ, 88 p. (1979).

Almeida, O.P., “Estudo do Efeito de Fronteiras Rígidas Sobre a Velocidade Terminal de

Partículas Isométricas”, Tese de M.Sc., Programa de Eng.Química, COPPE/UFRJ,
86p.(1995).

Barnea, E. e Mizrahi, J., “A Generalized Approach to the Fluid Dynamics of Particulate

Systems. Part I. General Correlation for Fluidization and Sedimentation in Solid Multiparticle
Systems”, Chem. Engineering J., vol. 5, 171-189 (1973).

Berker, R., “Intégration des Équations du Mouvement d'un Fluid Visqueux Incompressible”,
Handbuch der Physik (Flügge, S., Ed.), vol. VIII/2, Spring-Verlag, Berlin, 384 p. (1963).

Bird, R.B., Stewart, W.E. e Lightfoot, E.N., “Transport Phenomena”, J. Wiley, N.Iorque,
780 p.(1960).

Caswell, B. e Schwarz, W.H., “The Creeping Motion of a Non-Newtonian Fluid Past a

Sphere”, J. Fluid Mech., vol. 13, 417-426 (1962).

Caswell, B., “The Effect of Finite Boundaries on the Motion of Particles in Non-Newtonian
Fluids”, Chem. Engr. Science, vol. 25, 1167-1176 (1970).

Churchill, S.W., “The Development of Theoretically Based Correlations for Heat and Mass
Transfer”, Rev. Latino Am. Transf. Cal. Mat., vol. 7, 207-229 (1983).

Clift, R., Grace, J.R. e Weber, M.E., “Bubbles, Drops, and Particles”, Academic Press, N.
Iorque, 380 p. (1979).

Coelho, R.M.L. e Massarani, G., “Fluidodinâmica de Partículas: Ainda sobre Correlações em

Base aos Dados Experimentais de Pettyjohn e Christiansen”, Relatório LSP/COPPE 1/96
(1996).

Concha, F. e Barrientos, A., “Settling Velocities of Particulate Systems. 4. Settling of

Nonspherical Isometric Particles”, Int. J. Mineral Processing, vol. 18, 297-308 (1986).

Concha, F. e Christiansen, A., “Settling Velocities of Particulate Systems. 5. Settling

Velocities of Suspensions of Particles of Arbitrary Shape”, Int. J. Mineral Processing, vol 18,
309-322 (1986).

Francis, A.W., “Wall Effect in Falling Ball Method for Viscosity”, Physics, vol. 4, 403-406
(1933).

36
Govier, G.W. e Aziz, K., “The Flow of Camplex Mixtures in Pipes”, Van Nostrand, N.Iorque,
792 p. (1972).

Hackenberg, C.M., “Movimento Gravitacional Acelerado de Esfera Sólida em Fluido

Newtoniano no Regime de Stokes”, Caderno de Eng. Química da RBE, vol. 7 nº 2, 27-41
(1991).

Haider, A. e Levenspiel, O., “Drag Coefficient and Terminal Velocity of Spherical and
Nonspherical Particles”, Powder Technology, vol. 58, 63-70 (1989).

Lapple, C.E. e Shepherd, C.B., “Calculation of Particle Trajectories”, IEC, vol. 32, nº 5,
605-617 (1940).

Laruccia, M.B., “Velocidade de Sedimentação em Fluidos Não Newtonianos: Efeito da

Forma e da Concentração de Partículas”, Tese de M. Sc., Faculdade de Engenharia Química,
UNICAMP, 143 p. (1990).

Marchildon, E.K. e Gauvin, W.H., “Effects of Acceleration, Deceleration and Particle Shape
on Single-Particle Drag Coefficients in Still Air”, AIChE J., vol. 25, nº 6, 938-948 (1979).

Massarani, G. e Silva Telles, A., “Escoamento de Fluidos Não Newtonianos na Vizinhaça de

Partículas Sólidas”, Revista Brasileira de Física, vol. 8, nº 3, 150-166 (1979).

Massarani, G. e Santana, C.C., "Fluidização Homogênea: Caracterização Fluidodinâmica de

Sistemas Particulados com Porosidade Elevada", Congresso Europeu de Fluidização, Las
Palmas de Gran Canária, Espanha, 83-91 (1994).

Munroe (1988), citado em Fidleris, V. e Whitmore, R.L., “Experimental Determination of the

Wall Effect for Spheres Falling Axially in Cylindrical Vessels”, British J. of Applied Physics,
vol. 12, 490-494 (1961).

Pettyjohn, E.S. e Christiansen, E.B., “Effect of Particle Shape on Free-Settling Rates of

Isometric Particles”, Chem. Eng. Progress, vol. 44, nº 2, 157-172 (1948).

Politis, T. e Massarani, G., “Relative Solid-Fluid Velocity in Irregular Particle Systems”,

Fluid/Particle Separation J., vol 2, nº 4, 214-216 (1989).

Renganathan K., Turton, R. e Clark, N.N., “Accelerating Motion of Geometric and Spherical
Particles in a Fluid”, Powder Technology, vol. 58, 279-284 (1989).

Richardson, J.F. e Zaki, W.N., “Sedimentation and Fluidization: Part I”, Trans. Instn. Chem.
Engrs., vol. 32, 35-42 (1954).

Silva Telles, A. e Massarani, G., “Escoamento de Fuidos Não Newtonianos em Sistemas

Particulados”, Revista Brasileira de Física, vol. 9, nº 2, 535-547 (1979).

37
 Capítulo 2
A Decantação

1. A Trajetória da Partícula

O processo de separação sólido-fluido conduzido a partir de suspensões diluidas - a

decantação - pode ser analisado através do estudo da trajetória das partículas no interior do
equipamento de separação (Brauer, 1982).

Considera-se nesta análise que:

a) As partículas sejam caracterizadas individualmente através do diâmetro volumétrico

DP e da esferacidade φ ;

b) A distribuição de tamanhos das partículas, isto é, a análise granulométrica, seja

expressa por X = X ( DP ) , sendo X a fração em massa das partículas com diâmetro menor que
DP ;

c) O campo de velocidades do fluido não pertubado pela presença das partículas seja
u = u( x ) ;

d) Os efeitos da aceleração e concentração de partículas sejam desprezíveis no

comportamento dinâmico destas partículas.

Como visto no capítulo anterior, a equação do movimento de translação da partícula é

expressa por

0 = l + (ρ S − ρF )VP b (1)

l = A ρ Fc D U U , U = u − v (2)
2

DPUρ F
cD = f ( Re, φ), Re = , U= U . (3)
µ

Nestas equações ρ F e µ são respectivamente a densidade e a viscosidade do fluido, ρ S a

densidade das partículas, l a força resistiva que o fluido exerce sobre a partícula, A e VP a
área projetada ( πDP2 / 4 ) e VP o volume ( πDP3 / 6) da partícula, b a intensidade do campo
exterior, cD o coeficiente de arraste, u, v e U respectivamente a velocidade do fluido, a
velocidade da partícula e a velocidade relativa fluido-partícula.

Seja a situação simples em que se deseja determinar o diâmetro da partícula que

percorre a trajetória assinalada na figura (1) representando uma fenda retangular com
38
dimensões B, H e L. Na situação de maior interesse tecnológico H<<B, o que equivale a
considerar o escoamento como ocorrendo entre placas paralelas. O efeito da aceleração da
partícula não é levado em conta.

B L

H
h y u
v
x
(corte transversal) (corte longitudinal)

Figura 1 - Fluidodinâmica da partícula na fenda de seção retangular

Equação do movimento da partícula:

A
Componente na direção x , 0 = ρ F cDU (ux − v x ) + 0 ;
2

A
Componente na direção y , 0 = ρ F cDU (0 + v y ) − (ρ S − ρ F )VP g
2

Resulta da primeira equação que v x = ux e , portanto,

[ ]
1/ 2
U = (ux − v x ) 2 + (0 + v y ) 2 = vy .

Substituindo este resultado na segunda equação, vem

1/ 2
 
 (ρ S − ρ F )VP g 
vy =  
A
 ρ F cD 
 2 

que representa, segundo a equação (14) do primeiro capítulo, a velocidade terminal da

partícula isolada, vt . Portanto, desprezando o efeito da aceleração da partícula:

Na direção do escoamento do fluido, v x = ux ;

Na direção normal ao escoamento do fluido, v y = vt .

Voltando à figura da fenda de seção retangular, pela composição do movimento da

partícula

39
 h L
= , (4)
vt u h

onde u h
é a velocidade média do fluido em 0 ≤ y ≤ h . Portanto,

hu h 4 (ρ S − ρ F ) µg
vt = e cD / Re =
L 3 ρ2F vt3

que permite calcular Re (tabela 2 do capítulo 1) e dele o valor do diâmetro da partícula

desejado.

A situação mais desfavorável para a captura da partícula corresponde à posição h=H, a

espessura de separação da câmara. O diâmetro crítico D pc especifica as condições limites de
separabilidade no equipamento em análise: partículas com diâmetros maior que D pc são
coletadas com eficiência de 100% independentemente da posição em que ingressam na
câmara de separação. A equação (4) toma a forma

H L
= , (5)
vt u

a equação de projeto para a separação de partículas na fenda de seção retangular.

Exemplo

O sedimentador lamelado é constituido por um conjunto de fendas com seção

retangular em que a espessura de separação em cada fenda H/n é muito menor que o
comprimento L. A inclinação das lamelas, da ordem de 40º , permite a retirada contínua por
gravidade das partículas depositadas nestas lamelas.

L
H/n
H

θ
Suspensão

40
 A velocidade terminal da partícula com diâmetro crítico D pc é, analogamente à
equação (5),

H u H u
vt = ≅ .
L cos θ + H sen θ L cos θ

Sendo Q = BH u a vazão de suspensão que alimenta o sedimentador lamelado, resulta a

equação de projeto

vt = Q / A* ,

onde A* = nBL cos θ é a área projetada das n lamelas ativas no plano horizontal.

Exemplo

Deseja-se determinar o tempo consumido para que uma partícula se desloque, num
campo centrífugo, da posição radial r até a parede do equipamento de separação.

Ω
v
vθ
vr
r

Fluido

Na situação representada na figura, os componentes da velocidade do fluido são

ur = 0 e uθ = rΩ (Bird et al., 1960, p.96) e do campo centrífugo br = vθ2 / r e bθ = 0 , sendo ê
a velocidade angular da carcaça cilíndrica. Desprezando a aceleração da partícula, resulta da
equação do movimento:

1/ 2
 
dr  (ρ S − ρ F )VPbr 
vr = = vt =   (6)
dt A
 ρ F cD 
 2 

vθ = uθ = rΩ , (7)

41
onde vt é a velocidade terminal da partícula no campo centrífugo. A integração da equação
(6) para a partícula esférica e regime de Stokes leva ao valor do tempo desejado,

18µ
t= .ln(R/r) (8)
(ρS − ρF )ΩD 2

2. Separação Sólido-Fluido na Fenda de Seção Retangular

O projeto e a análise do desempenho do equipamento de separação sólido-fluido

podem ser realizados em base aos seguintes resultados:

a) Equação que relaciona o diâmetro de corte D∗ às propriedades físicas do sistema

particulado, às dimensões do equipamento e às condições operacionais;

b) Função eficiência individual de coleta relativa à partícula com diâmetro D,

η = η( D / D* ) (9)

que depende da configuração do equipamento, do regime de escoamento do fluido e da

dinâmica da partícula;

c) Função eficiência global de coleta que depende da distribuição granulométrica do

conjunto de partículas, X = X ( D) ,

1
η = ∫ η( D / D* )dX ; (10)
0

d) Equação que relaciona queda de pressão e vazão de fluido no equipamento de

separação.

O diâmetro de corte pode ser especificado de diferentes formas; neste texto é definido
como sendo o diâmetro das partículas que são coletadas com eficiência de 50% no
equipamento de separação.

Na análise da separação sólido-fluido em camada delgada ( H << B ) conduzida no

equipamento representado na figura (1) serão consideradas as seguintes hipóteses:

a) As partículas estão igualmente distribuidas na alimentação, x = 0 ,

independentemente do valor do diâmetro. Portanto, a eficiência de coleta da partícula com
diâmetro D que percorre a trajetória assinalada na figura é

η( D) = h / H , (11)

estando o diâmetro de corte associado a h = H / 2 .

42
 b) O escoamento de fluido na fenda é laminar, resultando (Bird et al., 1960, p.62)

 y  y 2
u = 6 u  −   (12)
 H  H  

1 h 1 1 h  h
u h = ∫ udy = 6 u  − 
h 0 2 3 HH (13)
u = u H/2 = u H

1 BH 3  ∆p 
Q = HB u = −  , (14)
12 µ  L 

onde Q é a vazão de fluido e ∆p a queda de pressão no equipamento.

c) Prevalece o regime de Stokes para as partículas sólidas (capítulo 1, tabela 4)

(ρ S − ρ F ) gD 2
vt = K1 (15)
18µ

K1 = 0, 843 log10 ( φ / 0, 065) .

Combinando as equações (4), (11), (13) e (15) resulta

L ( vt ) D
−1 2
h 1 h 1 uh 1  h  1 1 h   D 
η( D) = = = =  −    . (16)
H 2 ( H / 2) 2 L(vt ) D* 12  H  2 3 H    D* 
u

Portanto, a função eficiência individual de coleta η = η( D / D* ) para o equipamento em

questão, dentro das hipóteses consideradas, é

 1 D 
2
(3 − 2 η) η2 =  *  , D / D* ≤ 2
 2D  (17)
η = 1 *
, D/ D ≥ 2 .


A relação entre o diâmetro de corte D*, as propriedades físicas do sistema particulado,

as dimensões do equipamento e as condições operacionais pode ser estabelecida combinando
as equações (4) e (15)

1/ 2
* 9µQ 
D =  . (18)
 BLK1 g (ρ S − ρ F ) 

43
 Cabe ainda mencionar que quando o escoamento de fluido é turbulento,

u h
≅ Q / BH ,

resultando da equação (16) a denominada “eficiência teórica” do equipamento de separação

(Perry e Green, 1984, p.20-86):

 1 D 
2
η =  *  , D / D* ≤ 2
 2D  (19)
η = 1 , D / D* ≥ 2 .


3. O Conceito Sigma e a Especificação de Centrífugas

A trajetória da partícula assinalada no esquema da centrífuga tubular, figura (2),

permite especificar o valor D do diâmetro das partículas que são coletadas com eficiência de
100%. Para facilitar a análise, considera-se que as partículas sejam esféricas e que prevaleça o
regime de Stokes.

fluido fluido
u u

v L

Q
R0
R

Figura 2 - Esquema da centrífuga tubular

44
 Resulta da composição do movimento da partícula com diâmetro crítico, utilizando a
equação (8),

L L 18µ
t= =
Q
= ln (R / R 0 )
u (ρS − ρF )ΩD
π(R 2 − R 02 )

ou, explicitando a vazão de líquido,

(ρS − ρF )gD 2 π(R 2 − R 02 )LΩ 2

Q= ⋅ . (20)
18µ g ⋅ ln (R / R 0 )

Este último resultado mostra que a capacidade da centrífuga pode ser expressa pelo produto
de dois termos, um que caracteriza o sistema particulado (a velocidade terminal da partícula
no campo gravitacional) e o outro que caracteriza a configuração, as dimensões e rotação da
centrífuga, o fator sigma:

Q = vt Σ . (21)

A equação (21) constitui a base para a especificação da centrífuga para uma dada tarefa,
conhecendo o desempenho de uma centrífuga de laboratório, ambas do mesmo tipo, operando
com a mesma suspensão (Svarovsky, 1981):

 Q  Q
  =  . (22)
 Σ 1  Σ  2

4. Ciclones a Gás e Hidrociclones

A separação de particulas no interior do ciclone é efetuada pela ação do campo

centrífugo resultante da configuração do equipamento e do modo com que a suspensão o
alimenta.

O estudo da fluidodinâmica da partícula no ciclone vem recebendo contribuições

teóricas significativas, o que faz prever que em futuro próximo o projeto e a análise do
desempenho deste equipamento deixem de ser fundamentalmente empíricos: Leith e Licht
(1972), Bloor et al. (1980), Mothes e Löffler (1985), Barrientos e Concha (1992).

Tal como foi abordado no item 2 deste capítulo, procura-se estabelecer para ciclones
com diferentes configurações as equações que fornecem a relação entre diâmetro de corte,
propriedades físicas do sistema, dimensões do equipamento e condições operacionais, a
função eficiência de coleta relativa à partícula de diâmetro D, a expressão para a eficiência
global de coleta e a equação que relaciona vazão e queda de pressão no ciclone. Cabe
ressaltar que a configuração do ciclone caracteriza-se por uma relação específica entre suas
dimensões, expressa usualmente em termos do diâmetro da parte cilíndrica do equipamento,
Dc .

45
 Serão estudados neste item os ciclones a gás nas configurações Lapple e Stairmand e
os hidrociclones nas configurações Rietema e Bradley. Enquanto que os ciclones Lapple e
Stairmand são amplamente utilizados na indústria, os hidrociclones Rietema e Bradley
recebem o rótulo de equipamento de pesquisa e são distintos daqueles disponíveis
comercialmente (Pereira e Massarani, 1995).

As configurações

Estão especificadas na figura (3) as configurações dos ciclones a gás Lapple e

Stairmand, e na figura (4) as configurações dos hidrociclones Rietema e Bradley.

Diâmetro de corte na separação centrífuga

1/ 2
D*  µDc 
= K  ⋅ f ( RL ) ⋅ g (cv ) , (23)
Dc  Q(ρ S − ρ F ) 

onde Dc é o diâmetro da parte cilíndrica do ciclone, K um parâmetro que depende da

configuração, µ e Q são a viscosidade e a vazão de fluido que alimenta o ciclone, f é um fator
de correção que leva em conta o fato de que uma fração das partículas sólidas é coletada no
"underflow" sem a ação do campo centrífugo (efeito "T") e g um fator que leva em conta a
concentração volumétrica de sólidos na alimentação, cv (Massarani, 1991).

O fator f está relacionado ao quociente entre as vazões de fluido no "underflow"e na

alimentação, RL ,

f ( RL ) = 1 + AR L (24)

RL = B( Du / Dc ) C , (25)

e os parâmetros A, B, e C relacionados à configuração do ciclone, Du e Dc respectivamente os

diâmetros do "underflow" e da parte cilíndrica do equipamento.

Para partículas arredondadas o fator g pode ser expresso através da seguinte equação
empírica:

g (cv ) = 1 / [4,8(1 − cv ) 2 − 3,8(1 − cv )]0,5 . (26)

Os ciclones a gás operam com suspensões mais diluidas do que os hidrociclones e

freqüentemente a descarga de sólido é feita de modo intermitente a partir do barril acoplado
ao "underflow" do equipamento. Por estas razões, considera-se que para os ciclones a gás f e
g não influenciam o valor do diâmetro de corte, equação (23), ou seja, f = g = 1 .

Os valores dos parâmetros de configuração A, B, C e K estão reunidos na tabela (1),

cuja validade está restrita às condições operacionais assinaladas na própria tabela.

46
 Dc

Bc Do

Sc H c

Lc

Ciclone

Zc Lapple Stairmand
Bc / Dc 0,25 0,20
Do / Dc 0,50 0,50
Hc / Dc 0,50 0,50
Lc / Dc 2 1,50
Sc / Dc 0,62 0,50
Zc / Dc 2 2,50
Du / Dc 0,25 0,37

Du

47
 Figura 3 - Configuração dos ciclones a gás Lapple e Stairmand
Dc

Do
Di

l L1

Hidrociclone

Rietema Bradley
L Di / Dc 0,28 1/7
Do / Dc 0,34 1/5
L / Dc 5 -
L1 / Dc - 1/2
l/ Dc 0,40 1/3
θ 10º-20º 9º

Du

Figura 4 - Configuração dos hidrociclones Rietema e Bradley

48
Tabela 1 - Parâmetros de configuração do ciclone e condições operacionais recomendadas.

Configuração K A B C β u* ou Re** Du / Dc
(eq. 23) (eq. 24) (eq. 25) (eq. 25) (eq. 32)
Lapple 0,095 - - - 315 5 < u < 20m / s 0,25
Stairmand 0,041 - - - 400 10 < u < 30m / s 0,37
Rietema 0,039 1,73 145 4,75 1200 5 × 103 < Re < 5 × 104 0,10-0,30
Bradley 0,016 1,73 55,3 2,63 7500 3 × 103 < Re < 2 × 104 0,07-0,15

*u é a velocidade média do fluido na seção de entrada do ciclone, u = Q

Bc Hc
** Re = Dc ucρ F , onde u é a velocidade média do fluido na seção cilíndrica do
c
µ
Q
ciclone, uc = .
πDc / 4

Função eficiência individual de coleta no campo centrífugo

A eficiência individual de coleta relativa à partícula com diâmetro D pode ser

expressa pelas correlações empíricas:

Ciclones Lapple e Stairmand

( D / D* ) 2
η( D / D* ) = ; (27)
1 + ( D / D* ) 2

Hidrociclones Rietema e Bradley

exp(5D / D* ) − 1
η( D / D* ) = . (28)
exp(5D / D* ) + 146

Conhecida a distribuição granulométrica das partículas, X = X ( D) , é possível

estabelecer o valor da eficiência global de coleta no campo centrífugo,

1
I = ∫ ηdX (29)
0

e a eficiência global alcançada no ciclone, incluindo o efeito "T",

η = (1 − R L ) I + RL , (30)

sendo RL o quociente entre as vazões de fluido no "underflow" e na alimentação.

A integração da equação (29) para a situação bastante comum em que a distribuição

granulométrica pode ser representada pelo modelo de Rosin-Rammler-Bennet,

49
 'n
X ( D) = 1 − e − ( D / D ) , (31)

toma a forma (Massarani, 1991):

Ciclones Lapple e Stairmand

111
, n
0,118 + n D'
I= ⋅ ; (32)
1,81 − 0,322n + ( D' / D* ) D
*

Hidrociclones Rietema e Bradley

113
, n
0,138 + n D'
I= ⋅ . (33)
1,44 − 0,279n + ( D' / D* ) D
*

Cabe ressaltar que na equação (31) X é a fração em massa das partículas com diâmetro
menor que D e que D' e n são os parâmetros do modelo, respectivamente o diâmetro da
partícula que corresponde a X = 0, 632 e a dispersão.

A relação vazão - queda de pressão

A expressão clássica que relaciona vazão e queda de pressão na Mecânica dos Fluidos,
regime turbulento estabelecido, é utilizada também para os ciclones,

−∆p
β= (34)
ρ F uc2 / 2

Q
uc = , (35)
πDc2 / 4

sendo a queda de pressão medida entre o "overflow" e a alimentação. O valor de β depende da

configuração do ciclone, como mostra a tabela (1).

Exemplo

O diâmetro de corte na operação do ciclone

A separação sólido-fluido no ciclone pode ser considerada, numa análise grosseira,

como acorrendo em camada delgada, num campo centrífugo com intensidade constante. O
diâmetro de corte D* está associado à metade da espessura de separação, isto é, a Bc / 2 .

50
 trajetória da

Bc
partícula

Vista superior do
ciclone
Do
Dc

Admitindo que as partículas sejam esféricas e que prevaleça o regime de Stokes, resulta que o
tempo de residência do fluido e da partícula com diâmetro D* é dado por

Va Bc / 2
= (36)
Q (ρ S − ρ F ) D*2br
18µ

br = rΩ 2 = uΩ (37)

2πN e
Ω= . (38)
Va / Q

Nestes resultados, Va é o volume que o fluido ocupa no ciclone (pode se formar no

hidrociclone um nucleo de ar no interior do equipamento), u e Q respectivamente a
velocidade média na seção de entrada e a vazão de fluido, e Ne o número de espiras de fluido
que se formam no interior do ciclone. Combinando as equações (36) a (38):

1/ 2
* 9µBc 
D =  . (39)
 2 πN e u ( ρ S − ρ F ) 

No caso particular do ciclone Lapple, verifica-se por simples visualização que N e ≅ 5 .

Lembrando que para esta configuração

Q Q
u= =8 2,
Bc Hc Dc

vem para a equação (37),

1/ 2
D*  µDc 
= K 
Dc  Q (ρ S − ρ F ) 

51
 K = 0, 095 .

resultados que confirmam a equação (23) e o valor do parâmetro de configuração K, tabela

(1).

Exemplo

Deseja-se especificar uma bateria de ciclones Lapple para operar com 100m3 / min de
gás carregado com cinzas de carvão. Densidade das partículas de carvão, ρ S = 2,3g / cm3 .
São as seguintes as propriedades do gás: ρ F = 4,43 × 10−3 g / cm3 e µ = 0,035 cP . A bateria
deve funcionar com descarga de sólida intermitente e deseja-se uma eficiência global de
coleta superior a 85%. Distribuição granulométrica das partículas

1,5
X = 1 − e − ( D / 37 ,7 ) , D em µm.

Cálculo do diâmetro de corte

Resulta da equação (32), fazendo I = 0,85, D ′ = 37,7µm e n = 1,5 : D* = 6,0µm .

Estimativa de Dc e do número de ciclones em paralelo

Fazendo na equação (39) u = 15,0m / s , como recomendado na tabela (1), e lembrando

que B c ⋅ H c = D c 2 / 8 , vem que o diâmetro da parte cilíndrica do ciclone é dado por
Dc = 49,6cm . Portanto, sendo Q1 = uDc2 / 8 , resulta que o número de ciclones na bateria é

N = Q / Q1 = 3,6 .

Novo cálculo de Dc considerando 4 ciclones em paralelo

A vazão em cada ciclone é Q1 = Q / 4 = 25m3 / min . Vem da equação (23) que

Dc = 48cm , o que leva a uma velocidade u = 14,5m / s , valor este dentro da faixa
recomendada para a operação do ciclone Lapple.

Cálculo da potência do soprador

Considerando apenas a perda de carga nos ciclones, a potência requerida para a

separação é dada pela equação

Q∆p1
P= (40)
75 E

52
com P em cv, a vazão total Q em m³/s e a queda de pressão num ciclone ∆p1 em mm de
coluna de água. E é a eficiência elétrica do motor, da ordem de 0,5 para motores de baixa
potência. Resulta das equações (34) e (40) e da tabela (1): P = 2 cv .

Conclusões

A unidade: bateria com 4 ciclones Lapple em paralelo, diâmetro da parte cilíndrica

Dc= 48 cm.

Capacidade da unidade: 100m3 / min de ar carregado com cinzas de carvão

( D ′ = 37,7µm, n = 1,5) .

Eficiência global de coleta de partículas: 85%.

Potência do soprador, considerando apenas as perdas nos ciclones: 2 cv.

Problemas: Decantação

1. Calcular o diâmetro da menor partícula que é coletada com eficiência de 100% na câmara
de poeira abaixo esquematizada.

Propriedades físicas do fluido: densidade 1,2 × 10−3 g / cm3 e viscosidade 1,8 × 10−2 cP .
Propriedades físicas das partículas: densidade 2,5g / cm3 e esfericidade 0,7.
Dimensões da câmara: 2 × 2 ×16m , sendo a distância entre as lamelas de 10cm (a espessura
das lamelas é desprezível).

Vazão de suspensão na alimentação: 4 m3/s.

Considerar as seguintes situações diferentes:

a) A suspensão tem concentração volumétrica em sólido inferior a 0,2%;

b) Esta concentração é de 5%.

---- "Trajetória crítica"da menor partícula coletada com eficiência de 100%.

53
Resposta:

Velocidade terminal da menor partícula coletada com eficiência de 100%: 0,625cm / s .

Diâmetro da menor partícula coletada com eficiência de 100%, diluição infinita: 9,74µm .
Idem, 5% em volume de sólido: 11,0 µm.

2. Uma suspensão diluída de cal em água contém areia como produto indesejável.
Determinar, na operação a 25ºC:

a) A vazão de alimentação para a separação completa da areia no tanque com

dimensões 0,3 × 3 × 4m ;
b) O percentual de cal perdida na separação da areia.

Faixa granulométrica de areia: 70 < D p < 250µm .

Distribuição granulométrica das partículas de cal:
D p (µm) 20 30 40 50 60 70 80 100
100X 15 28 48 54 64 72 78 88.

Densidade da cal e da areia, respectivamente, 2,2 e 2,6 g / cm3 .

Esfericidade das partículas de cal e de areia, respectivamente, 0,6 e 0,8.

Resposta:

Capacidade do sistema para a separação completa de areia: 175 m3/h.

Diâmetro da maior partícula de cal no produto: 83,5 µm.
% de cal perdida na separação da areia: 18,8.

3. Foi conduzido no laboratório um ensaio de separação de argila (densidade 2,64g/cm3) de

uma suspensão aquosa em centrífuga tubular.

Propriedades do fluido: densidade 1g/cm3, viscosidade 1cP.

Dimensões da centrífuga de laboratório: Ro = 1,1cm; R = 2,2cm e L = 20 cm.

54
Número de rotações da centrífuga de laboratório: 20000 rpm.
Vazão de suspensão da centrífuga de laboratório que permite obter um classificado
satisfatório: 28,8 L/h.
Determinar a capacidade de uma centrífuga industrial operando com a mesma suspensão a
15000 rpm. Suas dimensões são: Ro = 5,21cm; R = 8,16cm e L = 73,4cm.

(Svarovski, L., “Solid-Liquid Separation”, Butterworths, Londres, 2ª edição, p. 196, 1981).

Resposta:

Fator Σ para a unidade de laboratório: 1,42×103 cm2.

Fator Σ para a unidade industrial: 4,22×104 cm2.
Capacidade da unidade industrial: 856 L/h.

4. A Companhia Chalboud do Brasil adquiriu uma bateria de ciclones com as dimensões

especificadas na figura para coletar partículas de um fluxo de ar a 70ºC e 1 atm. A densidade
das partículas é 1,05 g/cm3.
Verificar a validade da seguinte especificação fornecida pelo fabricante do equipamento:
partículas com diâmetro maior que 20 µm são coletadas com eficiência superior a 95%
quando a velocidade do ar na seção de alimentação do ciclone é 15 m/s.

55
Resposta:

Os ciclones fornecidos estão praticamente na configuração Lapple, podendo-se esperar uma

eficiência de coleta para as partículas de 20 µm de apenas 90%.

5. O ferro-velho "Dois Irmãos" da Pavuna dispõe de um conjunto de 3 ciclones em paralelo

na configuração Lapple, estado de conservação razoável. O diâmetro dos ciclones é 20 in.
Estimar:

a) A capacidade do conjunto para u = 15 m/s;

b) O diâmetro da partícula que é coletada com eficiência de 95%;
c) A potência do soprador a ser usado na operação.

Considerar que o gás tenha as propriedades físicas do ar a 200ºC e 1 atm e que as partículas
sólidas tenham densidade 3 g/cm3.

Resposta:

Capacidade da bateria de ciclones: 87 m3/min.

Diâmetro da partícula coletada com eficiência de 95%: 20 µm.
Potência do soprador: ~3cv (eficiência 0,5).

6. Deseja-se estudar o desempenho de uma bateria constituída por 2 ciclones Lapple em série
com respectivamente 63,6 cm e 45 cm de diâmetro no tratamento de 27,7 m3/min de gás
contendo 3% em volume de sólido.

Propriedades do gás: densidade 1,1x10-3 g/cm3 e viscosidade 1,7x10-2 cP.

56
Propriedades das partículas sólidas: densidade 2,5 g/cm3 e distribuição granulométrica dada
por

  D  1,5 
X = 1 − exp −  , D em µm .
  17,3 
Pede-se:

a) A eficiência global de coleta do sistema;

b) A potência do soprador para o serviço.

Resposta:

Diâmetro de corte relativo ao 1º ciclone, concentração volumétrica em sólido 0,03: 6,42 µm.
Diâmetro de corte relativo ao 2º ciclone, concentração volumétrica em sólido considerada
como nula: 3,48 µm.
1 ∗
Eficiência global de coleta o 1º ciclone: ∫0 η1 (D / D )dX = 0,690 .

1
Eficiência global de coleta o 2º ciclone: ∫0 (1 − η1 )η 2 dX = 0,162 .
(η1 e η2 são as eficiências individuais de coleta em cada ciclone)
Eficiência global do sistema: 85,2%.
Potência do soprador: ~3cv (eficiência 0,5).

7. Uma usina em Campos, RJ, pretende secar bagaço de cana com o gás de chaminé
proveniente da caldeira (propriedades do ar a 210ºC e 1 atm).

Especificar a bateria de ciclones Lapple para a recuperação de finos secos sabendo-se que a
vazão de gás é 140 m3/min e que as partículas maiores que 40 µm devem ser coletadas com
eficiência superior a 95%. A densidade do bagaço seco é 1,55 g/cm3.

Resposta:

Bateria constituída por 2 ciclones em paralelo com diâmetro 1 m.

8. Especificar a bateria de ciclones Lapple para operar com 100 m3/min de ar (520ºC e 1 atm)
contendo cinzas de carvão. A eficiência de coleta deve ser superior a 80%. Determinar
também a potência do soprador para a operação. A densidade das partículas sólidas é
2,3g / cm3 e a distribuição granulométrica é dada por

D(µm) 5 10 15 20 30 40
100X 12 27 48 63 80 88.

57
Resposta:

  D  1,35 
Distribuição granulométrica: X = 1 − exp −   , D em µm.
  21,5 

Diâmetro de corte para uma eficiência de coleta de 80%: 4,31 µm.

Bateria de 10 ciclones Lapple em paralelo com diâmetro 30 cm.
Potência do soprador: ~2cv (eficiência 0,5).

9. Estuda-se a possibilidade de reduzir o teor de cinzas de um carvão através da separação em

hidrociclone operando em fase densa. A alimentação contém 2 partes de carvão para 1 de
cinzas, em massa. A concentração volumétrica de carvão e cinzas na alimentação é de 5%.
Carvão e cinzas apresentam a mesma distribuição granulométrica,

  D 1,35 
X = 1 − exp −    , D em µm.
  21,5  

Estimar o teor em cinzas do concentrado de carvão ("overflow") que deve ser obtido numa
bateria de hidrociclones em paralelo com 2 in de diâmetro, nas configurações (a) Bradley e
(b) Rietema operando a uma queda de pressão de 45 psi. Fornecer também a capacidade de
cada hidrociclone.

Densidade do carvão e cinzas, respectivamente, 1,25 e 2,10 g/cm3.

Propriedades do fluido: densidade 1,21 g/cm3 e viscosidade 2,7 cP.

Resposta:

Fixando a relação entre os diâmetros de descarga "underflow" e da parte cilíndrica do

hidrociclone em 0,15, na operação a 45 psi:

Configuração Bradley Rietema

Capacidade por hidrociclone (m / h)3 1,90 4,74
Relação % vazões "overflow"/ alimentação 62,3 98,2
% carvão no "underflow" 58,8 41,1
% carvão no "overflow" 75,4 75,6
% do carvão da alimentação perdido pelo "underflow" 46,4 16,0

10. Deseja-se avaliar a possibilidade da utilização de uma bateria de hidrociclones Rietema no

beneficiamento do minério M. A suspensão aquosa a ser tratada contém 120 g/l de suspensão
do minério M e 45 g/l de suspensão de argila, produto indesejável. Temperatura de operação:
30ºC.

58
Diâmetro dos ciclones: 5 cm.
Densidade do minério e da argila: 2,7 e 2,1 g/cm3.
Distribuição granulométrica de minério e de argila:

  D 1,5    D 1,2 
X M = 1 − exp −    e X A = 1 − exp −    , D µm.
  21     3,5  
 

Dentro da faixa de condições operacionais recomendadas para o hidrociclone Rietema,

fornecer um "quadro desempenho" contendo: queda de pressão, capacidade, teor de minério
no produto do "underflow" e o percentual de minério perdido na operação.

Resposta:

Fixando a relação entre os diâmetros de descarga do "underflow" e da parte cilíndrica do

hidrociclone em 0,20:

∆p Q D*M D*A % Minério da

% Minério no
(m3/h) alimentação perdido
( atm ) (µm) (µm) "underflow"
pelo "overflow"
1 2,88 14,8 18,4 88,9 37,3
2 4,08 12,5 15,5 88,6 33,4
3 5,00 11,2 14,0 88,0 31,0
4 5,77 10,5 13,0 87,8 29,6

11. Deseja-se avaliar a possibilidade da utilização de uma bateria de hidrociclones Rietema

na separação pelo "underflow" das partículas com diâmetro maior que 15 µm, para posterior
operação de moagem, de uma suspensão de minério em água. O diâmetro dos ciclones é 5
cm.

Densidade do minério: 2,7 g/cm3.

% volumétrica em sólidos da alimentação: 6,6.
Distribuição granulométrica das partículas sólidas na alimentação:

  D 1,5 
X(D) = 1 − exp −    , D em µm.
  21  

Dentro da faixa de condições operacionais recomendadas para o hidrociclone Rietema,

fornecer um "quadro desempenho" do sistema de separação contemplando: queda de pressão,
capacidade, percentual de partículas com diâmetro maior que 15 µm no "underflow", as
perdas de "grossos" pelo underflow e as perdas de finos no "overflow".

59
Resposta:

As distribuições granulométricas da alimentação e do produto do "underflow" estão

relacionadas entre sí através da seguinte equação resultante de um balanço de massa

dX u
η = RL + η(1 − R L ) ,
dX a

onde:

η - Eficiência global de separação;

Xa e Xu - Fração em massa das partículas de diâmetro menor que D na alimentação e
no "underflow";
RL - Relação entre as vazões de "underflow" e de alimentação;
η - Eficiência individual de coleta relativa às partículas com diâmetro D (inclui apenas a
separação centrífuga e não o efeito "T").

Fixando a relação entre os diâmetros de descarga do "underflow" e da parte cilíndrica do

hidrociclone em 0,20:

% Total % D < 15µm

∆p Q D*M % D > 15µm D > 15µm da alimentação
( atm ) (m3/h) (µm) "underflow" perdido pelo perdido pelo
"overflow" "underflow"
1 2,88 14,8 95,5 26,0 6,0
2 4,08 12,5 95,5 15,0 6,6
3 5,00 11,2 95,3 8,5 7,3
4 5,70 10,5 95,0 5,9 7,9

Bibliografia

Barrientos, A. e Concha, F., “Phenomenological Model of Classification in Conventional

Hydrocyclones”. Cap. 21 em “Comminution - Theory and Practice” (Kawatra, K., Ed.)
American Institute of Mining, Metallurgical and Petroleum Engineers, N. Iorque, 287-305
(1992).

Bird, R.B., Stewart, W.E. e Lightfoot, E.N., “Transport Phenomena”, J. Wiley, N. Iorque, 780
p. (1960).

Bloor, M.I.G., Ingham, D.B. e Larerack, S.D., “An Analysis of Boundary Layer Effects in
Hydrocyclone”, Int. Conf. on Hydrocyclones, Cambridge, Inglaterra, 49-62 (1980).

60
Brauer, H., “Movement of Single Particles in Various Flow Fields”, em “Advances in
Transport Processes, vol. II” (Mujumdar, A.S. e Mashelkar, R.A., Eds.), Wiley Eastern
Limited, N. Deli, 352-432 (1982).

Leith, D. e Licht, W., “The Collection Efficiency of Cyclone Type Particle Collectors. A
New Theoretical Approach”, AICHE Symposium Serie, vol. 68, nº 126, 196-106 (1972).

Massarani, G., “Projeto e Análise do Desempenho de Ciclones e Hidrociclones 2”. Caderno

de Engª Química, RBE, vol. 8, nº2, 83-93 (1991).

Mothes, H. e Löffler, F., “Motion and Deposition of Particles in Cyclones”, Ger. Chem. Eng.,
nº 8, 223-233 (1985).

Pereira, C.M.S. e Massarani, G., “Análise do Desempenho de Hidrociclones”, XXIII

ENEMP,Maringá, PR, 690-694 (1995).

Perry, R.H. e Green, D.W., “Perry's Chemical Engineering Handbook”, McGraw-Hill, N.

Iorque, 6ª edição (1984).

Svarovsky, L., “Solid-Liquid Separation”, Butterworths, Londres, 2ª edição, 556 p. (1981).

61
 Capítulo 3
Escoamento de Fluidos em Meios Porosos

1. Equações da Continuidade e do Movimento para o Fluido

A fluidodinâmica da partícula isolada constituia nos capítulos anteriores o ponto de

partida na análise dos sistemas diluidos. A utilização desta mesma estratégia no estudo de
sistemas densos em partículas oferece, no entanto, dificuldade de tal ordem de grandeza que
tornaria, por exemplo, a simples estimativa da relação vazão e queda de pressão no
escoamento em meios porosos um problema transcendental. Como conhecer a posição de
cada partícula num sistema particulado tridimensional e como estabelecer as interações
interparticulares e com o fluido em escoamento?

Adota-se neste capítulo, bem como nos próximos, um modelo matemático com base
numa extensão da Mecânica do Contínuo para contemplar as misturas. As leis básicas de
conservação, que formam o núcleo da Teoria do Transporte em misturas, foram estabelecidas
por Truesdell (1957) e generalizadas mais tarde por Kelly (1964). Partindo destas equações e
na expressão do crescimento da entropia na forma da desigualdade de Clausius-Duhem,
diversos autores desenvolveram uma Teoria de Misturas capaz de descrever diferentes
fenômenos como a difusão molecular, as reações químicas e o escoamento de fluidos em
meios porosos (Crochet e Naghdi, 1966; Müller, 1971; Silva Telles e Fernandes, 1973).

As partículas sólidas constituem uma matriz porosa indeformável neste capítulo.

Para um referencial fixo à matriz, as equações da continuidade e do movimento para o fluido,
na forma integral, podem ser escritas do modo.

∂
∫ ∂t (ερ)dV + ∫ ρεu ⋅ ndS = 0 (1)
VR SR

∫ ερudV = ∫ TndS − ∫ m *dV + ∫ ερgdV , (2)

VR SR VR VR

onde SR e VR são respectivamente a superfície e o volume da região R que encerra matriz

porosa e fluido, ρ a densidade do fluido, ε a porosidade da matriz (fração volumétrica
ocupada pelo fluido), u a velocidade intersticial do fluido, T o tensor tensão que atua na fase
fluida, m* a força exercida pelo fluido sobre a matriz porosa (por unidade de volume do
sistema) e g a intensidade do campo exterior.

A velocidade superficial do fluido q, utilizada amplamente em substituição à

velocidade intersticial u, é medida desconsiderando a presença da matriz porosa.
Confundindo os conceitos de porosidade volumétrica e superficial,

q = εu . (3)

62
 A força de interação fluido-partícula m* pode ser decomposta na força resistiva m e
no empuxo,

m* = m − (1 − ε )ρg (4)

A parte isotrópica do tensor tensão que atua sobre o fluido no meio poroso é a
pressão p,

T = − p1 + τ , (5)

sendo τ a tensão extra.

As equações da continuidade e do movimento tomam a seguinte forma diferencial:

∂
( ερ) + div (ρq ) = 0 (6)
∂t

 ∂u 
ερ + (grad u)u = − grad p − m + div τ + ρg . (7)
 ∂t 

São formulados em seguida hipóteses constitutivas que, apesar de sua simplicidade,

parecem compatíveis com o conhecimento comum pertinente ao escoamento isotérmico de
um fluido homogêneo através de meios indeformáveis: a força resistiva m e a tensão extra τ
são, para um dado sistema matriz porosa-fluido, função da velocidade superficial q relativa a
um referencial fixo à matriz,

m = f (q ) (8)

τ = G (q ) (9)

As funções f e G não são inteiramente arbitrárias pois devem satisfazer a segunda lei da
termodinâmica e o princípio da invariança às mudanças de referencial. A primeira leva,
por exemplo, a

f ( 0) = 0 10)

e a segunda, no caso dos meios isotrópicos, que f e G são funções isotrópicas, isto é, para
todo o tensor ortogonal Q

f (Qq ) = Qf ( q ) (11)

G (Qq ) = QG (q )Q T . (12)

Teoremas de representação para estas funções são conhecidos e conduzem aos resultados
(Smith, 1971):
m = β( q ) q (13)

63
 τ = α1 ( q ) 1 + α 2 ( q ) q ⊗ q (14)

As equações (11) e (12) mostram que as funções escalares α1, α2 e β são, para um dado
sistema matriz porosa-fluido, função do módulo da velocidade superficial q e que esta e a
força resistiva m têm a mesma direção.

A força resistiva m

A experimentação conduzida nos últimos 150 anos, desde o trabalho pioneiro de

Henry Darcy (1856, apêndice D), fornece para a força resistiva m a seguinte equação
constitutiva:

µ  cρ k q 
m= 1 + q , (15)
k  µ 

onde µ é a viscosidade do fluido newtoniano, k e c são parâmetros que dependem apenas de

fatores estruturais da matriz porosa quando não ocorrem interações físico-químicas entre
matriz e fluido (Massarani, 1989, p. 37). k é a permeabilidade do meio poroso, com dimensão
L2, e c é um parâmetro adimensional.

A equação (15), a forma quadrática de Forchheimer, é válida para o escoamento

viscoso em meios isotrópicos homogêneos ou heterogêneos, isto é, meios em que k e c são,
respectivamente, constantes ou variáveis com a posição no sistema (Massarani, 1967). A
equação é válida também em condições não isotérmicas, verificando-se a variação da
viscosidade e da densidade do fluido ao longo do escoamento (Massarani, 1969).

A extensão da equação (15) para contemplar o escoamento de uma classe ampla de

fluidos não-newtonianos pode ser feita substituindo µ pela viscosidade efetiva,

µ ef = S ( λ* ) / λ* (16)

1,1 q
λ* = ⋅ , (17)
( tε )1 / 2 k

onde S é a função tensão cisalhante (uma propriedade material do fluido), λ* a taxa de

distensão característica (uma propriedade cinemática do escoamento) e t um fator estrutural
da matriz porosa com valor da ordem de 2,5 (Silva Telles e Massarani, 1979).

Na situação em que o escoamento de fluido na matriz porosa é lento,

cρ k q
<< 1, (18)
µ
a forma quadrática expressa pela equação (15) recai na forma linear

64
 µ
m= q (19)
k

amplamente conhecida como "lei" de Darcy, e o escoamento darcyano está associado à

validade desta "lei".

Cabe ainda assinalar que a ocorrência na natureza de meios anisotrópicos leva ao

interesse pelo estudo do escoamento nestes meios. A "lei" de Darcy toma a forma

m = µRq , (20)

onde R é o tensor resistividade (Ferrandon, 1948; Silva Telles e Massarani, 1975).

A tensão extra τ

O conhecimento atual sobre a tensão extra τ é extremamente escasso. Silva Telles e

Fernandes (1973) constataram através de um número limitado de experiências que na equação
(14) α 2 é praticamente nulo para fluidos newtonianos podendo, no entanto, atingir valores
significativos no caso de certas soluções poliméricas não-newtonianas. Na falta de
informação sobre α1, admite-se provisoriamente que nos casos usuais

τ = α1 ( q ) 1 . (21)

A equação de Darcy

Seja a situação comum em que o meio poroso isotrópico e homogêneo é percolado por
um fluido newtoniano. Resulta das equações (3), (7), (15) e (21):

 1 ∂q 1
ερ
 ε ∂t ε


[
+ 2 (grad q )q  = − grad p − α1 ( q ) ] − m + ρg (22)

µ  cρ k q 
m= 1 + q . (15)
k  µ 

No caso em que o escoamento é uniforme, isto é, quando o campo de velocidades q é

uniforme, a equação do movimento (22) toma a forma

0 = − grad p − m + ρg (23)

conhecida como equação de Darcy e utilizada indiscriminadamente na literatura sobre a

fluidodinâmica em meios porosos (Scheidegger, 1963; Bear, 1972). A utilização da equação
de Darcy é satisfatória mesmo na situação em que o escoamento de fluido é acelerado: o
termo de aceleração na equação (22) pode tomar valores significativos face ao termo resistivo

65
m apenas quando as partículas são grandes, com diâmetro da ordem de alguns milímetros
(Massarani, 1989, p. 45).

No escoamento incompressível a equação de Darcy toma a forma

− grad P = m . (24)

sendo P é a pressão piezométrica do fluido

P = p + hρg (25)

e h a distância (positiva na direção contrária a g) do ponto em questão medida a partir de um

plano horizontal de referência. No escoamento incompressível, sendo válida a "lei"de Darcy,

µ
− grad P = q (26)
k

que combinada à equação da continuidade, equação (6), leva ao resultado clássico da

hidráulica subterrânea (Polubarinova-Kochina, 1962)

∇2 P = 0. (27)

2. Propriedades Estruturais da Matriz Porosa

Como evidenciado no item anterior, a porosidade, a permeabilidade e o fator c são os

parâmetros que caracterizam a matriz porosa na percolação de um fluido homogêneo através deste
meio.

A determinação experimental de parâmetros estruturais

A porosidade pode ser determinada com o auxílio da picnometria simples, sendo

necessária a picnometria com vácuo nas medidas com meios consolidados que apresentam
porosidade reduzida.

A permeabilidade e o fator c são determinados experimentalmente por permeametria

através de um conjunto de medidas de vazão e queda de pressão efetuadas com a amostra

L
z

amostra
Q

P1 P2
Permeâmetro

66
A equação de Darcy, equação (23), toma a forma

dp µq cρq 2
− = + (28)
dz k k

na configuração do permeâmetro. A integração desta equação leva aos seguintes resultados

para os casos em que o escoamento é incompressível ou compressível e isotérmico de um gás
perfeito:

Escoamento incompressível,

1  ∆p  µ cρ
− = + q, (29)
q L  k k

Escoamento isotérmico de gás ideal,

ρ  ∆p  µ c
− = + G
G L  k k
G = ρq (30)
ρ + ρ2 M  ∆p 
ρ= 1 = p2 − , ∆P = p 2 − p1 .
2 RT  2 

As formas lineares (29) e (30) permitem calcular com facilidade os valores de k e c.

O modelo capilar

Apesar de sua simplicidade, o modelo capilar permite carrelacionar qualitativamente a

permeabilidade com alguns parâmetros estruturais da matriz porosa. A idéia de modelar o
meio poroso através de um feixe de dutos nasceu da analogia evidente entre a equação (26),
válida para o escoamento darcyano no meio poroso,

dP µ
− = q (26)
dz k

e a equação clássica da Mecânica dos Fluidos

dP µ
− = 2 u (31)
dz Rh / β

válida para o escoamento laminar e incompressível em dutos retilíneos. Na equação (31) u é a

velocidade média do fluido, Rh o raio hidráulico do duto, isto é, a razão entre a área da seção
de escoamento e o perímetro de molhamento, e β um fator adimensional que depende da
forma da seção transversal do duto, como mostra a tabela (1).

67
 Associando a velocidade u do fluido no duto à velocidade intersticial q/ε no meio
poroso, resulta das equações (26) e (31) a relação entre a permeabilidade e o raio hidráulico
da matriz porosa

k = εRh2 / β . (32)

O raio hidráulico depende da porosidade e da superfície específica do meio poroso:

Área da seção de escoamento

Rh = =
Perímetro de molhamento

(Área da seção de escoamento) (Comprimento do meio poroso)

Volume do meio poroso (33)
= =
(Perímetro de molhamento) (Comprimento do meio poroso)
Volume do meio poroso
ε
= .
SV

A superfície específica da matriz porosa depende da distribuição granulométrica e dos

fatores de forma B associado à superfície das partículas e C ao volume:

X = X ( D) (34)

S p = BD 2 , V p = CD3 (35)

onde X(D) é a fração em massa das partículas da amostra com diâmetro menor que D. Deste
modo, para a massa m de partículas com densidade ρ S ,

Superfície da matriz porosa

SV = =
Volume do meio saturado com o fluido

1 m / ρs
∫ 0 BD
2
⋅ dX
CD 3 B 1 dX
= = (1 − ε ) ∫
m / ρs C 0 D
1− ε (36)

68
Tabela 1 - Escoamento laminar e incompressível em dutos retilíneos: fator de forma β
(Berker, 1963)

Seção Transversal do Duto Fator β Observações

Círculo

2 —
a

Região anular entre

circunfe-
rências coaxiais 2(1 − α ) 2
[
1 + α 2 − (1 − α 2 / ln (1 / α) ] —
0≤ α <1
2≤β<3

αa
a

Elipse
π/2
αa (1 + α 2 ) π 2 E=∫ (1 + k 2 sen 2 φ)1/ 2 dφ
0
a
4E 2
k = (1 − α 2 )1/ 2
0< α ≤1
2 ≤ β < 2,46
Retângulo

16 16 1024 ∞ tg (am)
αa
a f = − 5 α∑
(1 + α ) 2 f 3 π n = 0 ( 2n + 1)
5

0 ≤ α ≤1 m = (2n + 1) π / (2αa )
1,78 ≤ β < 3
Cardióide
y

2 π 2 (1 + 2α 2 ) 3 π 4α 
1/ 2
x I = ∫ 1 − 2 cos θ dθ
θ (1 + 4α 2 )(1 + 4α 2 − 2α 4 ) I 2 0 4α + 1 
0 ≤ 2α < 1
2 ≤ β ≤ 2,23

x = cos θ + α cos 2θ
y = senθ + αsen2θ

69
Definindo o diâmetro médio de Sauter do modo

1 dX
D = 1/ ∫ , (37)
0 D

resulta para a superfície específica

B 1
SV = (1 − ε ) ⋅ ⋅ . (38)
C D

De modo coerente com os capítulos anteriores, a dimensão da partícula é especificada

pelo diâmetro volumétrico Dp e sua forma através da esfericidade φ ,

πD3p / 6 = V p (39)

Área superficial da esfera com diâmetro D p

φ= . (40)
Área superficial da partícula

Portanto,

B = π / φ, C = π / 6 (41)

resultando:

6(1 − ε )
SV = (42)
Dpφ

ε ( D p φ) ε
Rh = = (43)
SV 6(1 − ε )

e, finalmente,

εRh2 ( D p φ) 2 ε 3
k= = (44)
β 36β(1 − ε ) 2

1 dX
Dp = 1 / ∫ . (45)
0 Dp

O último resultado, conhecido na literatura como equação de Blake-Kozeny ou

Kozeny-Carman, permite correlacionar, no contexto do modelo capilar, a permeabilidade com
as propriedades das partículas e a porosidade do meio. A experimentação indica que o valor
do parâmetro estrutural β está compreendido entre 4 e 5 para meios com porosidade até 50%,

70
como mostra a figura (1). Para meios expandidos o valor de β aumenta significativamente
com a porosidade quando ε > 0,75 (Massarani e Santana, 1994).
O modelo capilar pode fornecer também informações qualitativas referentes ao fator c.
Neste caso, a analogia é estabelecida entre as equações que descrevem o escoamento a altas
vazões no meio poroso e no duto retilíneo:

dP cρ 2
− = q (46)
dz k

dP fρ 2
− = u (47)
dz 8R h

onde f é o fator de atrito no duto. Associando a velocidade u do fluido no duto à velocidade

intersticial q/ε no meio poroso, resulta das equações (32), (33) e (42)

Ω
c= , (48)
ε 3/ 2

sendo Ω um parâmetro adimensional a ser determinado experimentalmente,

Ergun (1952) : Ω = 0,14

0,35 < ε < 0,50 (49)
10 6 < k < 10 − 4 cm 2 .

0,98
 k 
0,37
k  
0,01
Massarani (1989, p. 52): Ω = 0,13 0  + 0,10 0  
  k  k 

−6 2
0,15 < ε < 0,75 k0 = 10 cm . (50)
10−9 < k < 10−3 cm2

A equação de Ergun (1952), extensamente utilizada na literatura de Engenharia

Química, é a expressão da equação de Darcy, equação (28), na qual a permeabilidade e o fator
c são representados por

( D p φ) 2 ε 3
k= 2
e c = 0,14 / ε 3/ 2 ,
150(1 − ε )
. (51)
∆P (1 − ε) 2 µq 1 − ε pq 2
− = 150 + 1,75 3
L ε3 ( D p φ) 2 ε ( D p φ)

71
 Fator β
5,5
1/8"cilíndros

5,0 1/32" placas

prismas
4,5 1/4"cilíndros

cubos
4,0

3,5 1/16"placas

3,2

0,30 0,32 0,34 0,36 0,38 0,40 0,42 0,44 0,46 0,48 0,50

Porosidade, ε
Figura 1 - Fator estrutural β (Coulson e Richardson, 1977, p. 12)

A literatura referente ao modelo capilar é riquíssima e não carece certamente de

imaginação. No início, o meio poroso era representado por um simples feixe de tubos
retilíneos de um mesmo diâmetro; mais tarde, o modelo foi se sofisticando: uma trança de
tubos com seção elítica variável. A passagem, por exemplo, da representação através do feixe
de tubos retilíneos para a configuração na forma de trança levou, no caso arbitrário do duto de
seção circular, ao "conceito" de tortuosidade, t,

t = β / β tubo = 0,5β (52)

que exprime o aumento do comprimento da trança quando retificada.

O grande mérito do modelo capilar está na possibilidade de fornecer, através de uma

analogia simples, resultados qualitativos que, complementados com a experimentação em
meios porosos, permitem correlacionar a estrutura do meio com suas propriedades
macroscópicas. A necessidade da experimentação vem do fato de que a estrutura da matriz
porosa é diferente daquela que um conjunto de tubos pode oferecer mesmo em configurações
complexas.

Exemplo

Deseja-se estimar o valor da permeabilidade e do fator c para o meio poroso

constituído por partículas que apresentam a seguinte análise de peneiras

Sistema Tyler Diâmetro médio da abertura D# Massa retida Fração em massa

(peneira nº) (mm) (g) retida ∆X
- 35 + 48 0,359 47,6 0,33
- 48 + 65 0,254 52,8 0,36
- 65 + 100 0,180 45,2 0,31
Total 145,6 1,00

72
A porosidade do meio é da ordem de 42% e a esfericidade das partículas é 0,78; confunde-se,
neste cálculo aproximado, os valores do diâmetro médio de peneiras D# e do diâmetro
volumétrico D p .

Cálculo do diâmetro médio de Sauter, equação (46)

1 1 1
Dp = ≅ = = 0,246mm .
1 dX  ∆X  0,33 0,36 0,31
∫ 0 Dp ∑  D  + +
0,359 0,254 0,180
i p i

Cálculo da permeabilidade, equação (44), com β ≅ 4,5

( D pφ) 2 ε 3
k= 2
= 5, 01 × 10−7 cm2 .
36β(1 − ε )

Cálculo do fator c, equação (50), com k0 = 10 −6 cm 2

[ ]
0,98
c = 0,13( k 0 / k ) 0,37 + 0,10( k 0 / k ) 0,01 / ε 3/ 2 = 1,01 .

3. Escoamento em Meios Porosos: Aplicações Clássicas

A maioria dos problemas relativos ao escoamento de fluido em meios porosos pode

ser resolvida a partir da forma simplificada da equação do movimento expressa pela equação
de Darcy,

0 = − grad p − m + ρg (23)

µ  cρ k q 
m= 1 + q , (15)
k  µ 

A equação de Darcy não leva em conta a aceleração do fluido na percolação através da matriz
porosa, o que parece ser aceitável na maioria dos problemas de interesse tecnológico, quando
a permeabilidade é inferior a 10−5 cm2 (Massarani, 1989, p. 45).

A perda de carga no meio poroso

No escoamento unidirecional e incompressível, as equações (15) e (23) levam à

equação da perda de carga no meio poroso,

∆P  WA  L  µ cρ q 
- =  = + q (53)
ρg  g  MP ρg  k k 

73
expressa em termos da altura de coluna de fluido que escoa no meio. Na equação (53) W A é a
energia dissipada devido ao atrito por unidade de massa de fluido.

Exemplo

Deseja-se calcular o valor do desnível H para que a vazão de água na coluna de

ionização seja 4m3/h (30ºC). A perda de carga na tubulação é 7,52m de coluna de água.
Dimensões da coluna: diâmetro Dc = 30cm e altura L = 100cm. Propriedades do meio poroso:
porosidade ε = 0,42 , permeabilidade k = 4 x 10-6 cm2 e fator c = 0,40.

O balanço de energia entre 2 pontos da instalação leva a (Perry e Green, 1984, p. 5-20)

∆p W W 
+ ∆z = − ∑ A  , (54)
ρg g i  g i

equação que encerra, respectivamente, a carga de pressão, a carga de altura, a carga da bomba
e a perda de carga na tubulação e no equipamento que compõe a montagem (inclui as colunas
de recheio). No caso em estudo, entre os pontos 1 e 2 assinalados na figura,
W
∆p = 0 e = 0 por não contar a instalação com uma bomba:
g

W 
H = ∑ A  . (55)
i  g i

Perda de carga no meio poroso

q = Q / A = 1,57cm / s ( velocidade superficial na coluna)

W L  µ cρq 
  =  + k  q = 3,70m de coluna de água
 g MP ρg  k 

74
Cálculo do desnível H

W  W 
H =  A + A = 7,52 + 3,70 = 11,22m
 g  tubo  g  MP

O escoamento compressível

A permeametria com gases foi analisada anteriormente neste capítulo, resultando a

equação (30) para o caso do escoamento isotérmico e unidimensional de gás ideal. Situações
mais complexas, como as que ocorrem nos secadores de grãos, podem ser também abordadas
partindo da equação de Darcy, com o conhecimento das equações de estado para o fluido. Os
problemas (8) e (16), incluidos no final deste capítulo, tratam, respectivamente, da expansão
adiabática e do escoamento isotérmico radial de gás ideal através da matriz porosa.

O escoamento transiente

Os problemas (9) e (11) propostos no final deste capítulo tratam do escoamento em

meios porosos sujeito a uma carga de líquido variável. A formulação nestes casos parte da
equação do movimento para o fluido à qual se associa, através das condições de contorno, o
balanço de massa transiente para o volume de fluido que alimenta o sistema.

4. O Escoamento Bifásico em Meios Porosos

As equações da continuidade o do movimento apresentadas neste capítulo podem ser

estendidas para contemplar a situação em que os poros da matriz são ocupados por 2 fluidos
imiscíveis:

∂
( εsi ρi ) + div(ρi qi ) = 0 (55)
∂t

0 = − grad pi − mi + ρi g, i = 1,2 . (56)

Nestas equações si é a saturação expressa pela fração volumétrica de poros ocupada pela fase i,

s1 + s2 = 1 (57)

e mi é a força resistiva (por unidade de volume de meio poroso) exercida pelo fluido i sobre a
matriz porosa e sobre o outro fluido.

O índice i = 1 denota o fluido que "molha" a matriz porosa, isto é, aquele que
preferencialmente recobre a superficie sólida. As pressões pi na equação (56) estão
relacionadas entre si através da pressão capilar pc (Bear, 1972, p. 441).

p2 − p1 = pc ( s1 ) (58)

75
que depende, dentro do ciclo embebição-drenagem, da natureza dos fluidos e de fatores
estruturais da matriz porosa. Como conseqüência, na situação em que s1 é constante ao longo
do escoamento

grad p1 = grad p2 .

No rasto das formas constitutivas para a força resistiva mi , seja , por hipótese, para
uma dada tríade fluido1 - fluido2 - matriz porosa

mi = mi ( s1 , q1 , q2 ), i = 1,2 . (59)

A função mi não é inteiramente arbitrária, pois deve satisfazer a segunda lei da

termodinâmica e o princípio da invariança às mudanças de referencial. A primeira conduz à
conclusão óbvia que

mi ( s1 ,0,0) = 0 ;

o segundo implica que para meios isotrópicos mi seja uma função isotrópica e que portanto
(Smith, 1971)

mi = Λ 1i (Ω) q1 + Λ 2i (Ω)q 2 + Λ 3i ( Ω)grad s1 , (60)

sendo, para uma tríade específica,

Ω = {s1 q1 , q2 , grad s1 , q1 ⋅ q2 , q1 ⋅ grad s1 , q2 ⋅ grad s1} (61)

Equação de Darcy-Buckingham

Na situação em que prevalece o escoamento lento das fases, a força resistiva, equação
(60) toma a forma linear

mi = α1i ( s1 )q1 + α1i ( s1 )q2 + α 3i ( s1 )grad s1

(62)
α12 = α 21 = 0,

que combinada à equação do movimento, equação (56), leva aos resultados:

1
q1 = −
α11 ( s1 )
[grad p1 + α 31 (s1 )grad s1 − ρ1g] (63)

1
q2 = −
α 22 ( s1 )
[grad p2 + α 32 ( s1 )grad s1 − ρ 2 g] . (64)

Seja o caso da percolação de água em solos não saturados em que o gradiente de

concentração de água é dominante face ao gradiente de pressão exercido sobre esta fase
(Tobinaga e Freire, 1980). Resulta da equação (63)

76
 g
q1 = K ( s1 ) − D( s1 )grad s1
g
, (65)
ρg α 31 ( s1)
K ( s1 ) = , D( s1 ) =
α11 ( s1 ) α11 ( s1)

conhecida na literatura como equação de Darcy-Buckingham. A condutividade hidráulica K e

o coeficiente de difusão D dependem da estrutura da matriz porosa.

Generalização da Forma Quadrática de Forchheimer

Nas condições em que o gradiente de saturação não é significativo, a generalização da

equação quadrática de Forchheimer para o escoamento bifásico toma a seguinte forma, a
partir da equação (60):

[
mi = γ 1i ( s1 ) + γ 2i ( s1 ) q1 + γ 3i ( s1 ) q 2 q1 + ] (66)
+ [δ 1i ( s1 ) + δ 2i ( s1 ) q1 + δ 3i ( s1 ) q 2 ]q 2 , i = 1,2

O resultado é compatível com as tradições da Engenharia de Reservatório, através das

equações de Muskat, e com os dados experimentais relativos ao escoamento bifásico na
coluna de recheio utilizada nas operações de absorção, destilação e extração (Tobinaga,
1979).

Sendo o escoamento dos fluidos lento, a experimentação indica que apenas os

coeficientes γ 11 e δ12 são significativos,

µ1
γ 11 ( s1 ) = (67)
kk 1 ( s1 )

µ2
δ12 ( s1 ) = (68)
kk 2 ( s1 )

onde k é a permeabilidade da matriz porosa e ki (s1) a permeabilidade relativa da fase i que

depende, no ciclo de embebição-drenagem, da estrutura desta matriz porosa (Scheidegger,
1963, p. 222).

As equações de Muskat são as equações do movimento para o escoamento bifásico

lento em meios porosos. Resulta das equações (56), (66) a (68):

kk i ( s1 )
qi = − (grad pi − ρi g) . (69)
µi

A operação das colunas de recheio, amplamente utilizadas nos processos de

transferência de massa entre fases, caracteriza-se pelo escoamento bifásico em altas vazões
através de meios com permeabilidade e porosidade elevadas. Estudos conduzidos por
Tobinaga (1979) no escoamento concorrente gás-líquido levou aos seguintes resultados:

77
 (1 − εs1) 2 µ1
γ11 = α1 ⋅
(εs1)3 k

(1 − εs 2 ) 2 µ 2
δ12 = α 2 ⋅
(εs 2 )3 k
1 − εs1 ρ1
γ 21 = β1 ⋅
3
(εs1) k1 / 2
1 − εs 2 ρ2
δ32 = β2 ⋅
(εs 2 )3 k1 / 2

sendo que os coeficientes αi e βi dependem da estrutura da matriz porosa. Os coeficientes γ31,

δ11, δ21, δ31 são aparentemente pouco importantes na formulação das equações.

O valor da queda de pressão na operação em contra-corrente da coluna de recheio é

geralmente estimado através de correlações empíricas, como aquela apresentada na figura (2)
(Catálogo DC-11, Norton Chemical Process Products, 1977). O problema (13) da coletânea
reunida no final do capítulo trata da estimativa da queda de pressão numa coluna operando
com selas intalox de cerâmica, 1½".

78
Figura 2 - Queda de pressão no escoamento bifásico contracorrente em coluna recheada
(Norton, DC-11, 1977)

C - 10,8 10

ρ G ( ρL- ρ G)
0,1
6 lin

CG Fν
F - Fator de forma ha

2
4 de
125 i nu Parâmetro: queda de
G - Velocidade mássica do nd
aç pressão em mm água/m
gás (kg/m2.s) 2 83 ão de recheio
1 40
L - Velocidade mássica do
líquido (kg/m2.s) 0,6 21
0,4 8
v - Viscosidade cinemática
do líquido (cSt). 0,2
ρG - Densidade do gás 4
(kg/m3). 0,1

0,06
ρL - Densidade do líquido 0,04
(kg/m3)
0,02

0,01

0,01 0,02 0,04 0,06 0,1 0,2 0,4 0,6 1 2 4 6 10

L ρG ) 1/2
(
G ρL

Fator de forma F
Dimensão nominal (in)
Recheio Material
1/4 3/8 1/2 5/8 3/4 1 1¼ 1½ 2 3 3½
Hy-pak Metal 43 18
Super Intalox Cerâmica 60 30
Super Intalox Plástico 33 21 16
Anel Pall Plástico 97 52 40 24 16
Anel Pall Metal 70 48 33 20 16
Intalox Cerâmica 725 330 200 145 92 52 40 22
Anel Raschig Cerâmica 1600 1000 580 380 255 155 125 95 65 37
Anel Raschig Metal 1/32” 700 390 300 170 155 115
Anel Raschig Metal 1/16” 410 290 220 137 110 83 57 32
Sela Berl Cerâmica 900 240 170 110 65 45

79
Problemas: Escoamento de Fluido em Meios PorosoS

1. Determinar os valores da permeabilidade e do fator c a partir dos dados experimentais

obtidos por permeametria.

a) Meio de areia artificialmente consolidado com 5% de araldite.

Granulometria da areia: -14+20 # Tyler.
Fluido: água (densidade 1 g/cm3 e viscosidade 1,18 cP).
Comprimento do meio: 2,1 cm.
Área da seção de escoamento: 16,8 cm2.
Porosidade do meio: 0,37.
Dados de velocidade superficial e queda de pressão:
q ( cm / s) 6,33 7,47 10,2 12,7 15,2 17,7 20,3 23,9
−∆p( cmHg ) 4,69 6,24 10,4 15,2 21,2 28,0 35,9 48,9.

b) Meio não consolidado de areia.

Granulometria da areia: −35+48 # Tyler.
Fluido: ar a 25ºC e pressão atmosférica na descarga.
Comprimento do meio: 33,4 cm.
Área da seção de escoamento: 5,57 cm2.
Porosidade do meio: 0,44.
Dados de velocidade mássica e queda de pressão:
G ( g / cm 2 s) 1,59×10-3 5,13×10-3 9,49×10-3 12,3×10-3 22,4×10-3 44,6×10-3 70,3×10-3
−∆p(cm água) 6,40 20,8 38,6 50,3 92,5 197 321.

Em relação ao segundo caso, estimar os valores da permeabilidade e do fator c pelas

correlações da literatura. Considerar que a esfericidade da areia seja 0,70.

80
Resposta:

a) Permeabilidade: 6,9×10−6 cm2.

Fator c: 1,10.
b) Permeabilidade: 1,29×10−6 cm2.
Fator c: 0,75
Valores estimados (caso b):

( Dφ ) 2 ε 3
k= 2
= 1, 02 × 10−6 cm2
170(1 − ε )

0,98
  k0 
0,37
 k0 
0,01 
c = 0,13  + 0,10   ε3 / 2 = 0,81 , k0 = 10−6 cm2 .
  k   k  

2. Especificar a bomba centrífuga para a unidade de tratamento de água constituída por um

filtro de carvão (A), coluna de troca catiônica (B) e coluna de troca aniônica (C). Capacidade
da instalação: 6 m3/h. A tubulação tem 35 m de comprimento (aço comercial, 1½"# 40) e
conta com uma válvula globo (aberta) e 12 joelhos de 90º. O desnível entre os pontos 1 e 2 é
praticamente nulo. Temperatura de operação: 25ºC.

Especificação das colunas:

Coluna Altura de recheio (cm) Diâmetro (cm)

A 50 50
B 90 55
C 90 55

81
Especificação do recheio:

Coluna Granulometria Esfericidade Porosidade

# Tyler
−35+48 (30%)
A −48+65 (40%) 0,60 0,42
−65+100 (30%)
B Dp = 0,45 mm 0,85 0,37

C Dp = 0,70 mm 0,85 0,38

Resposta:

Perda de carga na tubulação incluindo acidentes: 4,45 m.

Perda de carga nas colunas:

Velocidade Diâmetro Permeabilidade Perda

Coluna superficial médio (cm2) Fator c de carga
(cm/s) (cm) (m)
A 0,849 2,45×10−2 2,80×10−7 1,16 14,7
B 0,702 4,5×10−2 1,10×10−6 1,03 5,71
C 0,702 7×10−2 2,97×10−6 0,820 2,17

Especificação da bomba: Vazão, 6 m3/h; carga, 27 m; potência, 2cv.

3. Calcular a vazão de água que a bomba centrífuga Minerva (5cv) fornece à instalação
abaixo esquematizada constituída por uma coluna recheada, 25 m de tubulação de aço de 1½”
(#40), válvula gaveta (1/4 fechada), válvula de retenção e 7 joelhos de 90º. Desnível entre os
pontos 1 e 2: -3m. Temperatura de operação: 25ºC.

A coluna recheada: diâmetro 20 cm e altura 1 m.

O recheio: diâmetro médio e esfericidade das partículas 450 µm e 0,85; porosidade do leito
0,38.

82
Curva característica da bomba:

Vazão (m3/h) 0 2 4 6 8 10 12 14
Carga (m) 53,9 53,9 53,3 52,2 50,6 46,7 41,7 32,8.

Resposta:

Perda de carga na tubulação incluindo acidentes: 7,26×10−2 Q2 m de coluna de água, com a

vazão Q expressa em m3/h.
Desnível entre os pontos 2 e 1: 3 m.
Propriedades do meio poroso: permeabilidade 1,23×10−6 cm2 e fator c 0,97.
Perda de carga na coluna recheada: 6,62 Q+0,701 Q2 m de coluna de água, com a vazão Q
expressa em m3/h.
Carga da instalação: 3+6,62 Q+0,774 Q2 m de coluna de água, com a vazão Q expressa em
m3/h.
Vazão fornecida pela bomba: 4,7 m3/h (compatibilidade entre as cargas da bomba e da
instalação).

4. Estimar a capacidade (m3/m2h) do filtro de areia abaixo esquematizado operando com água
a 20ºC. A primeira camada, com porosidade 0,37, é constituída de areia com a seguinte
granulometria:

Sistema Tyler % em massa

−14+20 20
−20+28 60
−28+35 20

A segunda camada, com porosidade 0,43, é constituída de brita com 1,3 cm de diâmetro. A
esfericidade da areia e da brita pode ser considerada como sendo 0,7.

83
Resposta:

Propriedades do leito de areia: diâmetro médio das partículas 0,70 mm, permeabilidade
1,80×10−6 cm2 e fator c 0,94.
Propriedades do leito de brita: permeabilidade 1,19×10−3 cm2 e fator c 0,38.
Capacidade do filtro: 15,1 m3/m2h.

5. Estimar o tempo consumido na percolação de 100L de óleo através de um leito de carvão

ativo com porosidade 0,42. A pressão de ar comprimido é 7 atm manométricas.

Propriedades do óleo: densidade 0,85 g/cm3 e viscosidade 35 cP.

Dimensões do leito de carvão ativo: diâmetro 30 cm e altura 50 cm.
Propriedades das partículas de carvão: esfericidade 0,6 e granulometria.

Sistema Tyler % em massa

−35+ 48 15
−48+ 65 65
−65+100 20

Resposta:

Propriedades do leito de carvão: diâmetro médio das partículas 245 µm, permeabilidade
2,8×10−7 cm2 e fator c 1,16.
Vazão de óleo que percola o leito de carvão: 4,78 L/min.
Tempo consumido na percolação de 100L de óleo: 21 min.

6. Calcular a queda de pressão no reator catalítico em leito fixo sabendo que opera isotermicamente
a 550ºC e que a pressão na descarga do reator é 1,5 atm. A porosidade do leito é 0,44.

Vazão mássica de gás com propriedades do nitrogênio: 200 kg/h.

Dimensões do leito de catalisador: diâmetro 30 cm e altura 40 cm.
Propriedades das partículas de catalisador: esfericidade 0,65 e distribuição granulométrica dada
por

1
X= , Dp em µm.
3
 123 
1+  
 Dp 
 

84
Resposta:

Propriedades do leito de catalisador: diâmetro médio das partículas 102 µm,

permeabilidade 7,08×10−8 cm2 e fator c 1,56.
Queda de pressão no reator: 8,08 atm.

7. Um conversor secundário de ácido sulfúrico tem 2,3 m de diâmetro e opera com 3 camadas
de catalisador, perfazendo um total de 1,35 m de altura. As partículas de catalisador são
cilindros eqüiláteros com diâmetro 9,5 mm. A porosidade do leito é 35%. Calcular a queda
de pressão no conversor sabendo-se que a velocidade mássica de gás é 2,6 × 103 kg / m2 h .
A alimentação é feita a 400ºC, resultando uma temperatura na descarga de 445ºC. A pressão
na descarga do conversor é 1 atm.
Composição do gás:

SO3 SO2 O2 N2
Alimentação (% molar) 6,6 1,7 10,0 81,7
Descarga (% molar) 8,2 0,2 9,3 82,3
(Coulson, J.M. e Richardson, J.F.; "Chemical Engineering", Pergamon Press, Londres, 2ª
edição, vol. 2, p. 737, 1968).

Resposta:

O escoamento será considerado como sendo isotérmico a 423ºC. Massa molecular média
32,6. A viscosidade da mistura será considerada como sendo a do nitrogênio a 423ºC e 1 atm:
0,033 cP.
Propriedades do leito de catalisador: diâmetro e esfericidade 1,09 cm e 0,87, permeabilidade
5,42×10−4 cm2 e fator c 0,54.
Queda de pressão no conversor: 40,6 cm de coluna de água.

8. Análise da expansão adiabática de um gás perfeito através de um meio poroso aberto à

atmosfera. Estabelecer a relação entre pressão no reservatório de volume V e o tempo,
admitindo que o escoamento no meio poroso seja darcyano e que a viscosidade do fluido
possa ser considerada como sendo uma constante no processo em questão. Condições iniciais
no reservatório: pressão e temperatura po To .

85
Resposta:

Expansão adiabática de gás perfeito: p −a T = po−a To = C , constante,

onde a = ( c p − cv ) / c p .
Equação do movimento do fluido escoando no meio poroso, forma integrada:

1 M k p 2 −a − patm
2−a
G= ⋅ ⋅ ⋅ , onde L é o comprimento do meio poroso.
2 − a CR µ L

Balanço de massa de gás no reservatório:

VM −a dp
− GA = (1 − a ) ⋅ ⋅p , onde A é a área da seção de escoamento.
CR dt

Combinando os resultados, resulta por integração:

1 Akt I p /p x −a
⋅ = , onde I = ∫ o atm dx .
(1 − a )( 2 − a ) µLv patm p / p atm x 2 − a − 1

Exemplo: Expansão adiabática de ar, a = 0,286,

86
9. O filtro abaixo esquematizado recebe uma vazão constante Q de líquido newtoniano.
Estabelecer a relação entre a velocidade superficial de fluido que escoa através do meio
poroso e o tempo de percolação. Estabelecer também o tempo em que ocorrerá o
transbordamento do filtro. O escoamento no meio poroso pode ser considerado darcyano. Na
condição inicial o meio poroso está saturado com líquido e l = l o

Resposta:

Forma integrada da equação do movimento do líquido escoando no meio poroso:

µ
(l-L)ρg/L = q ( t ) , onde q é a velocidade superficial do fluido.
k

Balanço de massa na camada de líquido:

Q
− q = dl/dt.
A

Combinando os resultados e integrando com a condição inicial l(0) = l 0, obtém-se a relação

entre a espessura de líquido l e o tempo de percolação t:

l = β + (lO − β ) exp(−αt )
α α

k ρg Q
α= , β = − Lα .
µL A

Substituindo o resultado na equação que fornece o balanço de massa na camada de líquido,

resulta a relação entre a velocidade superficial q e o tempo de percolação:

Q β
q= + ( lO − )α exp (−αt )
A α

87
 Q
O transbordamento do líquido no filtro se dará quando > L α, no tempo T correspondente a
A
l= H:
β l
− O
1 α
T = ln , H > lO.
α β
−H
α

10. Comparar os valores estimados de queda de pressão com os resultados experimentais

obtidos por Silva Telles e Massarani (RBF, 9, 2, 535, 1979) para o escoamento de solução
aquosa de Natrosol através de um meio poroso de areia artificialmente consolidada.
Considerar o escoamento como sendo darcyano.
Propriedades do fluido: densidade 1 g/cm3 e relação taxa de distensão - tensão cisalhante dada
por:

λ( s −1 ) 0 150 300 500 700 900 1100 1300 1500

S(dyn/cm )2 0 44,5 70,3 98,5 123 145 166 185 203
Propriedades do meio poroso: comprimento 2,0 cm, permeabilidade 1,4×10−6 cm2,
porosidade 0,38.

Dados experimentais de velocidade superficial do fluido e de queda de pressão no escoamento

através do meio poroso:

q(cm/s) 0,275 0,524 1,07 1,63

−∆p(cmHg ) 7,63 11,6 18,4 24,8

Resposta:

Relação entre a velocidade superficial do fluido e a taxa de distensão característica:

q(cm/s) 0,275 0,524 1,07 1,64

q
λ* = (s-1 ) 232 443 904 1380
k

“Modelo reológico” para 232 < λ < 1380 s −1: S = 1,63λ0,66 dyn / cm 2 .
Cálculo da viscosidade efetiva e estimativa da queda de pressão:

q(cm/s) 0,275 0,524 1,07 1,63

*
S( λ )
µ= ( P) 0,256 0,205 0,161 0,140
λ*
−∆pexp (cmHg) 7,63 11,6 18,4 24,8

−∆psim (cmHg ) 7,91 11,7 18,7 24,8.

88
11. Análise do dreno vertical darcyano. Estabelecer a relação entre a vazão de água que
percola através do dreno e o tempo de percolação, sabendo-se que as paredes são porosas e
que no tempo inicial a altura de líquido é ho.

Resposta:

Equação do movimento do fluido escoando através da parede porosa:

dp µ 1 dQ
− = ⋅ ⋅
dr k 2πr dz
− ∆p = (h - z)ρg
µ 1 R dQ
(h - z)ρg = ⋅ ⋅ ln 2 ⋅
k 2π R1 dz
πkρgh 2
Q= .
R2
µln
R1

Balanço de massa de fluido no poço:

dh
− Q( t ) = πR 12
dt
1 1 k ρgt
∴ = + .
h ho 2 R2
µ R 1 ln
R1

89
Relação entre vazão de fluido e o tempo de percolação:

πkρg 1
Q= ⋅
R 2
µ ln 2  
R1  1 kρgt 
 + 
 h0 µR 2 ln R 2 
 1
R1 
 

12. Análise do dreno-aleta darcyano. Estabelecer a relação entre a capacidade do dreno e as

propriedades e dimensões do meio poroso.

Resposta:

Expressões para a vazão no dreno:

∫( ) ∫ r (− ∂z )o, r dr.
L R
k k
Q = 2πR − ∂P dz = 2π ∂P
µ ∂r z , R µ
o o

90
Formulação do campo de pressões piezométricas:

1 ∂ ∂P 2
   ∂ P
r  + =0
 r ∂r  ∂r  ∂z 2
 ∂P 
  = 0 , P(z, R ) = p atm − zρg
 ∂r  z, o
 ∂P 
P(o, r ) = p atm + Hρg ,  =0
 ∂z  L, r

O campo de pressões piezométricas:

  2n − 1  
∞  sen  π 
2ρg H  2   sen (λ z ) I o (λ n r ) ,
P(z, r ) = p atm + Hρg − ∑ 
L n =1 λ n
+
2 
n
I o (λ n R )
λn
 
 
π
λn = (2n - 1), n = 1,2,3,...
2L

Capacidade do filtro:

 π(2n − 1) R 
I1 
Rk H ∞
2  2n − 1   1  2 L 
Q = 8ρg  L2 ∑  + sen π  .
 L  µ n =1 L π(2n − 1)  2  2n − 1 I  π(2n − 1) R 
o
 2 L 

13. Uma coluna de absorção opera com 680 kg/h de gás e 9000 kg/h de líquido, fluxos em
contracorrente. Os fluidos têm as propriedades do ar e da água a 20ºC e 1 atm. Diâmetro da
coluna: 0,42 m. Recheio: Intalox 1½", cerâmica. Calcular a queda de pressão na coluna
sabendo que a altura de recheio é 9 m.

Resposta:

A queda de pressão na fase gasosa é 75,4 cm de coluna de água.

14. Um resfriador de partículas pode operar nas configurações abaixo assinaladas. Calcular a
vazão de ar fornecida pelo compressor radial Minuano nestas duas situações. As partículas
são esferas com diâmetro 0,5 mm e a porosidade do leito é 0,40.
Dimensões do equipamento: 30×30×60 cm.
Propriedades do ar a 20ºC e 1 atm.
ρF < u >2
Queda de pressão na tubulação de 2": ∆p = 3 .
2

91
Curva característica do compressor:

Q(m3/min) 1 2 3 4
C(mm água) 1870 1200 600 0.
S S

30
30

G G

30 30 (cotas em cm)

Configuração 1 Configuração 2

Resposta:

Propriedades do meio poroso: permeabilidade 2,61×10−6 cm2, fator c 0,78.

Queda de pressão na tubulação: − ∆p = 1,21 Q 2 cm água, com a vazão expressa em m3/min.
Queda de pressão no resfriador, 1ª configuração: − ∆p = 19,1Q + 1,49 Q 2 cm água, com a
vazão expressa em m3/min.
Queda de pressão no resfriador, 2ª configuração: − ∆p = 76,7Q + 11,9 Q 2 cm água, com a
vazão expressa em m3/min.
Vazão na 1ª configuração: 2,75 m3/min (meios em paralelo)
Vazão na 2ª configuração: 1,55 m3/min (meios em série)

15. O problema de Dupuit. Desprezando os efeitos de capilaridade, estimar a vazão de

líquido newtoniano que percola, em escoamento darcyano, a parede porosa. A largura da
parede é B.

92
Resposta:

Seja o caso particular em que ( H1 − H2 ) << L . Nesta situação, a posição da interface ar-
líquido pode ser expressa por uma reta.

k dP
q=− ⋅
µ dx

k d
Q = − Bh ⋅ ⋅ (p atm + hρg)
µ dx

Q kρg H 22 − H12
=− ⋅ .
B µL 2

16. Estabelecer a relação entre a vazão mássica W (M/θ) e a queda de pressão no escoamento
permanente e isotérmico de um gás perfeito através da configuração abaixo esquematizada.

Resposta:

Equação do movimento do fluido escoando no meio poroso:

dp µ  cρ k 
0=− − 1 + qq
dr k  µ 

W = 2πrLρq = 2πrLG

dp µW 1 c W2 1
−ρ = ⋅ + ⋅ ⋅
dr 2πLk r k (2πL) r 2
2

93
 µW R c W2  1 1 
ρ(− ∆p) = ln 2 + ⋅  − 
2πLk R 1 k (2πL) 2  R1 R 2 

onde

M  p1 + p 2  M  ∆p 
ρ=  =  p2 − , ∆P = p 2 − p1 .
RT  2  RT  2 

94
Bibliografia

Bear, J., “Dynamics of Fluids in Porous Media”, American Elsevier, N. Iorque, 764 p. (1972).

Berker, R., “Intégration des Équations du Mouvement d'un Fuide Visqueux Incompressible”,
Handbuch der Physik (Flügge, S., Ed.), vol VIII/2, Spring-Verlag, Berlin, 384 p. (1963).

Coulson, J.M. e Richardson, J.F., “Tecnologia Química II: Operações Unitárias”, Fundação
Calouste Gulbenkian, Lisboa, 889 p. (1977).

Crochet, M.J. e Naghdi, P.M., “Constitutive Equation for Flow Through an Elastic Solid” Int.
J. Engng. Sci, vol. 2, 303 (1966).

Darcy, H., “Les Fontaines Publiques de la Ville de Dijon”, Dalmont, Paris,590 p. (1856).

Ergun, S., “Fluid Flow Through Packed Columns”, Chem. Eng. Progress, vol. 48, n° 2, 89-94
(1952).

Ferrandon, J., “Les Lois de l'Écoulement de Filtration”, Génie Civil, vol. 125, 24 (1948).

Kelly, P.D., “A Reacting Continuum” Int. J. Engng. Sci., vol. 2, 129 (1964).

Massarani, G., “Critique de la Validité de la Loi de Darcy en Milieu Hétérogène”, C.R. Acad.
Sc., Paris, t. 205, A, 587-589 (1967).

Massarani, G., “Critique de la Validité de l'Équation de Darcy pour un Écoulement Non-

isotherme”, C.R. Acad. Sc., Paris, t. 269, A, 812-815 (1967).

Massarani, G., “Aspectos da Fluidodinâmica em Meios Porosos”, Revista Brasileira de

Engenharia, número especial, 96 p. (1989).

Massarani, G. e Santana, C.C., “Fluidização Homogênea: Caracterização Fluidodinâmica de

Sistemas Particulados com Porosidade Elevada”, Congresso Europeu de Fluidização, Las
Palmas de Gran Canária, Espanha, 83-91 (1994).

Müller, I., “Equation for the Flow of Fluid Through a Porous Solid”, CISM Meeting at Udine,
Itália (1971).

Perry, R.H. e Green, D.W., “Perry's Chemical Engineering Handbook”, McGraw-Hill, N.

Iorque, 6ª edição (1984).

Polubazinova-Kochina, P.Ya., “Theory of Ground Water Movement”, Princeton University

Press, Princeton, 613 p. (1962).

Scheidegger, A.E., “The Physics of Flow Through Porous Media”, University of Toronto
Press, Toronto, 313p. (1963).

95
Silva Telles, A. e Fernandes, R.C., “Normal Stresses in Flows Through Porous Media”, 43rd
Annual Meeting of the Society of Rheology, Cincinatti (1973).

Silva Telles, A. e Massarani, G., “Escoamento de Fluidos Através de Meios Porosos

Anisotrópicos”, Revista Brasileira de Física, vol. 5, nº 1, 1-15 (1975).

Silva Telles, A. e Massarani, G., “Escoamento de Fluidos Não-Newtonianos em Sistemas

Particulados”, Revista Brasileira de Física, vol. 9, nº 2, 535-346 (1979).

Smith, G.F., “On Isotropic Functions of Symmetric Tensors, Skew-Symmetric Tensors and
Vectors”, Int. J. Engng. Sci, vol. 9, 899 (1971).

Tobinaga, S., “Escoamento Bifásico em Meios Porosos”, Tese de D.Sc., Programa de

Engenharia Química, COPPE/UFRJ, 203 p. (1979).

Tobinaga, S. e Freire, J.T., “Escoamento de Água em Solos Não-Saturados: Um Descrição

pela Teoria de Misturas”, RBT, vol. 11, nº 2, 105-115 (1980).

Truesdell, C., “Sulle Basi della Termodinamica”, Rend. Accad. Licei, vol. 22, 33 e 158
(1957).

96
 Capítulo 4
Fluidodinâmica em Sistemas Particulados Expandidos

1. Equações da Continuidade e do Movimento

A formulação para descrever a fluidodinâmica em sistemas particulados expandidos,

como acorre na fluidização, sedimentação livre e transporte de partículas, pode ser
estabelecida a partir das equações da continuidade e do movimento para cada fase e mais as
equações constitutivas. Analogamente ao capítulo anterior (Massarani, 1989, p. 84):

Equações da continuidade,

∂
( ερ F ) + div( ερ F v F ) = 0 (1)
∂t

∂
∂t
[(1 − ε)ρS ] + div[(1 − ε)ρSvS ] = 0 ; (2)

Equações do movimento,

 ∂v 
ερ F  F + (grad v F )v F  = − grad p − m + ρ F g (3)
 ∂t 

 ∂v 
(1 − ε )ρ S  S + (grad v S )v S  = div TS + m + (1 − ε )(ρ S - ρ F ) g ; (4)
 ∂t 

Equações constitutivas,

µ  cρ ε k U 
m = F 1 + F  εU , U = v F − v S (5)
k  µF 

TS = − p S ( ε )1 . (6)

Nestas equações, ε é a porosidade da matriz (fração volumétrica ocupada pelo fluido),

ρ F e ρ S a densidade do fluido e do sólido, v F e v S a velocidade intersticial das fases fluida e
sólida, p e pS a pressão no fluido e no sólido, m a força resistiva que o fluido exerce sobre a
matriz sólida (por unidade de volume de sistema particulado), g a intensidade do campo
exterior e TS a tensão exercida sobre a fase sólida.

A força resistiva m foi amplamente estudada na literatura e, como conseqüência, o

resultado expresso pela equação (5) mostra-se satisfatório nas situações usuais (Massarani,
1989). A experimentação torna-se, no entanto, particularmente difícil quando a porosidade do
sistema alcança valores nas proximidades de 1, como ocorre no transporte pneumático
diluído.

97
 O panorama em relação à tensão TS exercida sobre a fase sólida é bastante complexo
face às dificuldades na realização de ensaios fundamentais em Mecânica dos Sólidos. A
forma constitutiva expressa pela equação (6), apesar da simplicidade, é utilizada com
resultados satisfatórios na fluidodinâmica dos meios expandidos e no estudo da filtração e
sedimentação unidimensional com formação de substratos deformáveis (Massarani, 1989, p.
86; Gidaspow, 1994, p. 31). Como conseqüência, na situação em que a porosidade é
constante,

div TS = − grad pS ( ε ) = 0 . (7)

Formas complexas para a tensão nos sólidos se fazem necessárias quando se trata do
estudo de um leito deslizante de partículas, como ocorre na descarga de um silo e na operação
dos leitos recirculante e de jorro (Gidaspow, 1994, p. 31).

As equações do movimento para as fases fluida e sólida podem tomar formas

diferentes segundo as tradições de cada "escola" e as peculiaridades dos temas abordados, o
que acarreta um certo desconforto na tentativa de compatibilizar os resultados da literatura.
As equações (3) e (4) apresentadas neste texto satisfazem o conhecimento comum e os
resultados experimentais conhecidos, não tendo sido assinalado qualquer contra-exemplo.

a) A equação (3) conduz na estática ao resultado bem conhecido

0 = − grad p + ρ F g .

b) A equação (3) leva à definição da força resistiva m quando prevalece o escoamento

uniforme,

m = − grad p + ρ F g

e dela resultam, a partir das medidas de queda de pressão e vazão, os parâmetros estruturais k
e c da equação (5).

c) A integração da equação (3) completa, incluindo o termo da aceleração do fluido na

matriz porosa, conduz para os escoamentos incompressíveis convergente ou divergente
efetuados em calota esférica porosa aos resultados (Massarani, 1974):

Escoamento convergente,

ρ F QF2 4 µ F QF cρ F QF2
2 4
(1 − α ) = ( p1 − p0 ) − (1 − α ) − 2 3
(1 − α 3 ); (8)
8π εR0 2 πkR0 12 k π R0

Escoamento divergente,

ρ F QF2 4 µ F QF cρ F QF2
− 2 4 (1 − α ) = ( p0 − p1 ) − (1 − α ) − (1 − α 3 ) . (9)
8π εR0 2 πkR0 2
12 k π R0 3

98
Nestas equações QF é a vazão volumétrica de fluido, R0 e R1 os raios menor e maior da
calota e α = R0 / R1 . A Figura (1) mostra que os dados experimentais confirmam a validade
das equações (8) e (9) e, portanto, a validade da equação do movimento para o fluido,
equação (3).

d) A fluidização de um sistema particulado tem inicio quando no escoamento de fluido

a força resistiva iguala o peso aparente de sólido por unidade de volume,

m = (1 − ε )(ρ S − ρ F ) g . (10)

Este resultado bem conhecido pode ser resgatado

ε também a partir da equação do movimento para o sólido,
L equação (4), e da forma constitutiva para TS , equação
z
(6), considerando que na fluidização incipiente a
porosidade seja uniforme. De modo equivalente, a
distribuidor combinação das equações do movimento para as fases
leva, no caso da fluidização, ao resultado

QF

∆P
− = (1 − ε )(ρ S − ρ F ) g (11)
L

onde P é a pressão piezométrica no fluido.

99
Queda de pressão, mm de carga de água

110

Eq. (8)

90

70
p Eq. (9)

50
Q

30

Resultados experimentais
10

200 400 600 800 1000 1200

3
Vazão de fluido (cm /s)

Figura 1 - Escoamento acelerado de fluido (ar a 27,5ºC) em matriz porosa

( R0 = 4,5mm, R1 = 35,5mm, ε = 0,37, k = 1,1 × 10 -4 cm 2 , c = 0,55) (Massarani, 1974).

100
2. Caracterização dos Meios Expandidos

Seja abordada a situação em que os campos de velocidades da fases fluida e sólida e a

distribuição de porosidades são uniformes, como deve acorrer na fluidização homogênea, na
sedimentação livre e no transporte hidráulico vertical de partículas. A equação do movimento
para a fase sólida, equação (4), toma a forma

µFε cρ F ε 2 2
U+ U = (1 − ε )(ρ S − ρ F )g , (12)
k k

sendo a velocidade relativa U, respectivamente,

QF
U= (fluidização homogênea) (13)
εA

v
U = S (sedimentação livre) (14)
ε

QF QS
U= − (transporte vertical homogêneo). (15)
εA (1 − ε ) A

Nestes resultados, QF e QS são as vazões volumétricas de fluido e de sólido, A a área da

seção transversal de escoamento e vS a velocidade de deslocamento da interface líquido-
suspensão na sedimentação livre.

As relações entre porosidade, permeabilidade e fator c, parâmetros estruturais

presentes na equação (12), podem ser determinadas com facilidade através da fluidização com
líquido conduzida em duto de seção constante: o aumento da vazão de líquido acima das
condições de fluidização mínima acarreta uma expansão uniforme do leito, como mostram os
resultados obtidos por Gubulin e Massarani (1977) na fluidização de areia ( −28 + 35# Tyler )
com água (20ºC, velocidade superficial 2,04 cm/s),

Distância ao
distribuidor (cm) 10 20 30 40 50 60 70

Porosidade média
na seção (raios gama) 0,821 0,820 0,824 0,822 0,826 0,821 0,830.

Agrupando os parâmetros estruturais ε, k e c nos termos φ1 e φ2 , resulta da equação (12):

µ F φ1 ( ε )U + ρ F φ 2 ( ε )U 2 = (1 − ε )(ρ S − ρ F ) g . (16)

onde φ1 ( ε ) = ε / k , φ 2 ( ε ) = cε 2 / k .

101
 A fluidização conduzida com líquido de viscosidade elevada (velocidade relativa U
reduzida) leva a φ1 ; conhecido este termo, a fluidização com água leva φ2. As funções φ1 e φ
2 para sistemas específicos constituidos por cilindros equiláteros, esferas de vidro e areia
estão representadas na figura (2).

Cabe ainda ressaltar que a fluidodinâmica em meios expandidos fornece valiosas

informações sobre a dependência na porosidade do parâmetro β da equação de Kozeny-
Carman apresentada no capítulo anterior,

( D p φ) 2 ε 3
k= (17)
36β(1 − ε ) 2

( D p φ) 2 ε 2
β= φ1 ( ε ) , (18)
36(1 − ε ) 2

onde D p é o diâmetro médio de Sauter relacionado à distribuição granulométrica das

partículas e φ é a esfericidade destas partículas.

Os dados experimentais reunidos por Restini (1977), Santana e Massarani (1974) e

Massarani e d'Ávila (1974) mostram, para partículas com diferentes formas, que o valor de β
mantem-se constante para porosidades de até 75% e cresce de forma acentuada à medida que
a porosidade aumenta.

102
 5
10
φ1(ε) φ2(ε)
-2
(cm ) -1
4 (cm )
6-10 2
10

4
3-10 60

30

4
10

3
6-10
10

3
3-10 6

3
3
10

2
6-10
1
0,5 0,6 0,7 0,8 0,9 ε 0,5 0,6 0,7 0,8 0,9 ε
Figura 2 - Funções φ1 e φ2 para cilindros equiláteros com diâmetro 1,14 mm Ο (Restini,
1977), esferas de vidro com diâmetro 0,58 mm • (Santana e Massarani, 1974) e areia com
diâmetro 0,20 mm (Massarani e d`Ávila, 1974).

3. O Elo Entre a Fluidodinâmica de Partículas e a Teoria de Misturas

O resultado expresso pela equação (16) permite estabelecer o elo de ligação entre a
fluidodinâmica em meios porosos e a fluidinâmica de partículas.

O capítulo 1 fornece para a partícula isolada:

103
 1,18
4 (ρS − ρ F )D p g  24 
0,85
 0,85 
cD = =   + K2 (19)
3 ρFv∞ 2  K1 Re ∞  
 

D p v∞ ρ F
Re ∞ =
µ (20)
K1 = 0,843 log10 ( φ / 0,065), K2 = 5,31 − 4,88φ.

Nestas equações, cD e Re ∞ são o coeficiente de arraste e o número de Reynolds da

partícula com diâmetro volumétrico D p , v∞ a velocidade terminal da partícula isolada
movimentando-se no fluido com densidade ρF e viscosidade µF e φ é a esfericidade da
partícula.

Combinando as equações (16), (19) e (20):

3 µ 2F
µ F φ1 ( ε )U + ρ F φ 2 ( ε )U 2 = (1 − ε ) c Re∞2
3 D
4 ρ F Dp
(21)
2
U U 3
ψ 1 ( ε ) Re∞   + ψ 2 ( ε ) Re∞2   = c D Re∞2
 v∞   v∞  4

sendo

2
ε Dp ε 2 cD p
ψ 1 (ε) = e ψ 2 (ε) = .
1− ε k 1− ε k

O coeficiente de arraste é função do número de Reynolds, equação (19), e portanto,

vem da equação (21) que

U
= f ( Re ∞ , ε) . (22)
v∞

Este resultado, o mesmo apresentado no capítulo 1 como correlação empírica do tipo

Richardson e Zaki, estabelece o elo de ligação entre a fluidodinâmica da partícula e a
fluidodinâmica em meios porosos.

Exemplo

A influência da porosidade no valor do parâmetro β presente na equação de Kozeny -

Carman, equação (17), pode ser estabelecida com o auxílio de dados obtidos na fluidização
homogênea:
2
ε D p 36β(1 − ε )
ψ1 ( ε ) = = .
1− ε k ε 2φ2

104
Seja a situação em que as partículas são esféricas e Re∞< 0,2. Neste caso,

24 18
cD = , ψ1 = .
Re (U / v∞ )

e, em base aos dados reunidos por Concha e Almendra (Tabela 6C do capítulo 1),

U 0,83ε 3,94 , 0,5 < ε ≤ 0,9

=
v∞ 4,8ε − 3,8 , 0,9 < ε < 1.

Resulta para partículas esféricas:

 −1,94
0,60 ε , 0,5 < ε ≤ 0,9
 1− ε
β=
 ε2
 , 0,9 < ε < 1.
 2(1 − ε )(4,8ε − 3,8)

O valor do parâmetro β cresce com o aumento da porosidade a partir de ε = 0,80 ,

ε 0,6 0,7 0,8 0,9 0,95 0,98

β 4,06 4,01 4,64 7,36 11,9 26,6.

4. Transporte Hidráulico e Pneumático de Partículas

O transporte de partículas sólidas por arraste em fluido conduz, de um modo geral, à

formação de um campo de porosidades heterogêneo na seção transversal de escoamento da
mistura sólido-fluido. Em algumas situações, no entanto, dependendo da natureza do
problema em estudo, a formulação para o transporte de partículas pode ser substancialmente
simplificada considerando que a mistura comporta-se como um fluido homogêneo (Santana,
1982; Gidaspow, 1994):

a) Transporte pneumático vertical em fase densa (fluidização incipiente) ou em fase

diluída (porosidade superior a 95%); transporte hidráulico vertical sem restrições;

b)Transporte hidráulico em qualquer configuração no caso em que as partículas são

pequenas, verificando-se o critério empírico de Newitt

1800gDv ∞
Ne = < 1. (23)
3
VM

Na equação (23) D é o diâmetro do tubo, v∞ a velocidade terminal das partículas no fluido de

arraste e VM a velocidade da mistura sólido-fluido,

105
 QS + QF
VM = , (24)
A

onde QS e QF são a vazão volumétrica de fluido e de sólido e A a área da seção transversal

de transporte.

A diferença entre as formulações para os casos (a) e (b) reflete a dificuldade na

medida das propriedades reológicas da suspensão constituida por partículas relativamente
grandes (caso a) e pelo fato de que nesta situação o valor da velocidade relativa fluido-
partícula no transporte pneumático pode ser significativamente maior que zero. Como
conseqüência, o valor da porosidade no transporte depende da fluidodinâmica do sistema
particulado.

Apesar destas considerações, há o concenso bem cristalizado na literatura de que o

projeto e o estabelecimento das condições operacionais das linhas de transporte hidráulico e
pneumático não podem prescindir de estudos conduzidos em unidade piloto bem
instrumentada (Krauss, 1980; Wasp, 1983).

Transporte vertical homogêneo: partículas "grandes"

Os efeitos causados pela aceleração do sistema particulado não são considerados e o

transporte é, por exemplo, vertical ascendente.

Equação do movimento para a mistura homogênea:

p2 ∆p fV M2 ρ M
− = + ρM g (25)
L 2D

L DVM ρ M
f = f (Re M , e / D), Re M = (26)
µM

p1 ρ M = ερ F + (1 − ε )ρ S = (1 − ε )(ρ S − ρ F ) + ρ F (27)
z Equação do movimento para o fluido no sistema particulado,
equação (16), que permite calcular a porosidade no transporte:

S F µ F φ1 (ε) U + ρ F φ 2 (ε) U 2 = (1 − ε)(ρS − ρ F )g

QF QS (28)
U= − .
Aε A(1 − ε)

Nestas equações, ρM e µM são a densidade e a viscosidade da mistura sólido-fluido e f o fator

de atrito na interação fluidodinâmica entre a mistura e a parede do duto onde ocorre o
transporte.

106
 A comparação entre os valores do gradiente de pressão calculados através da equação
(25) e os valores resultantes da experimentação conduzida no transporte pneumático em fase
densa (Santana, 1982), no transporte pneumático em fase diluida (Ferreira et al., 1996) e no
transporte hidráulico (Restini, 1977) parecem indicar que a viscosidade e o fator de atrito da
mistura podem ser expressos pela viscosidade e o fator de atrito do fluido, este último
representado pela equação clássica

1  e  6,81  
0 ,9

= −2 log10 0,27 +   , (29)
f  D  Re M  
 

onde e/D é a rugoridade relativa do duto.

Do ponto de vista da compatibilidade entre as formulações apresentadas neste

capítulo, cabe ressaltar que a combinação das equações do movimento para as fases fluida e
sólida, equações (3) e (4), leva à parcela da interação sólido-fluido no gradiente de pressão
total,

 ∆p 
−  = (1 − ε)(ρS − ρ F )g + ρ F g = ρ M g ,
 L  SF

resultado este incluído na equação (25).

Exemplo

São comparados na tabela os valores experimentais e calculados do gradiente de

pressão no transporte hidráulico vertical de cilindros equiláteros de plástico
( D p = 1, 14 mm , ρ S = 1, 77 g / cm3 ) com água a 30ºC (Restini, 1977). Tubulação de PVC, 1".
Os resultados mostram, neste pequeno conjunto de dados, que a formulação proposta permite
estimar o gradiente de pressão com erro da ordem de 7%.

Dados e resultados experimentais Resultados calculados

(1)  ∆p  (2) (3) (4)

−  ρM g
WF WS  L  exp fV M2 ρ M  ∆p 
ε exp ε cal − 
( g / s) ( g / s)
3
( dyn / cm )  L  cal
2
( dyn / cm ) 2D
2 2
( dyn / cm ) ( dyn / cm )

137 59,4 0,773 1138 0,780 14,7 1150 1160

984 239 0,875 1260 0,876 208 1075 1283
2004 654 0,845 1677 0,843 422 1100 1522
1310 177 0,920 1246 0,928 305 1035 1340
1348 92,0 0,975 1251 0,962 147 1010 1154

107
(1) Medida experimental da porosidade obtida pela técnica de atenuação de raios gama.
(2) Porosidade calculada com o auxílio da equação (16) e dados experimentais de φ1 e φ 2
apresentados na figura (2).
(3) Fator de atrito calculado através da equação (29).
(4) Gradiente de pressão calculado através da equação (25).

Transporte hidráulico homogêneo

Depedendo das condições fluidodinâmicas, o transporte hidráulico de partículas

"pequenas" (Ne<1, equação 23) pode ser formulado do mesmo modo que o escoamento de
fluidos homogêneos com características não-newtonianas (Massarani e Silva Telles, 1992). O
balanço global de energia entre dois pontos da instalação leva a:

∆p
+ g∆z = W − W A (30)
ρM

f ( ∑ L )V M2
WA = (31)
2D
1  e  6,81  
0,9
DVM ρ M
= −2 log10 0,27 +   , ReM =
f  D  ReM   µ ef
 
Q + QS
VM = F , ρ M = ερ F + (1 − ε )ρ S
A

QF
ε= (velocidade relativa entre as fases nula).
QF + QS

A viscosidade efetiva da suspensão, µ ef , depende da natureza da mistura, através da

relação entre taxa de distensão (λ) e tensão cisalhante (S), e das condições de escoamento,

V
λ* = 6,40 M (taxa de distensão característica) (32)
D

µ ef = S ( λ* ) / λ* . (33)

Nas equações (30) e (31), W e W A são, respectivamente, a energia (por unidade de

massa de suspensão) fornecida pela bomba instalada e a energia dissipada por atrito no
escoamento; (ΣL) é o comprimento equivalente total da instalação, incluindo dutos e
acidentes. Cabe ainda ressaltar que a expressão para a taxa de distensão característica,
equação (32), tem natureza empírica e que o procedimento sugerido para o cálculo da
viscosidade efetiva independe do tipo do “modelo reológico” associado à suspensão.

Finalmente, cabe mencionar que a Samarco (1977) opera no Brasil com a mais extensa
instalação conhecida no mundo para o transporte de minério de ferro. O mineroduto liga a

108
Mina de Germano (MG) à Usina de Pelotização de Ponta do Ubu (ES), numa extensão de
396 km, dutos com 50 cm de diâmetro: 12 milhões ton/ano de concentrados finos são
transportados na forma de uma polpa contendo 65% em massa de sólidos. Em relação à
reologia destas suspensões, Coelho et al. (1981) mostraram que a concentração de sólidos e a
taxa de distensão afetam o valor da viscosidade efetiva, como mostra a figura (3). O
problema (12) da coletânea reunida no final do capítulo trata da estimativa do valor da
potência de bombeamento necessária para a operação do mineroduto da Samarco.
8
µef /µF

6
ε= 0,62

ε= 0,73

2
ε= 0,86

0
0 1000 2000 3000 4000 5000
-1
Taxa de distensão (s )
Figura 3 - Viscosidade efetiva de suspensão de minério de ferro em água (Coelho et al., 1981).

Problemas: Fluidodinâmica em Sistemas Particulados

Expandidos
1. Os seguintes dados resultaram da fluidização com ar de dolomita em tubo com 10 cm de
diâmetro:
Vazão de Ar Queda de Pressão Altura do Leito
(L/min) (cm de coluna de água) (cm)
21,3 74,4 57,3
17,2 73,0 55,3
14,3 71,7 54,3
11,4 69,3 51,2
9,54 67,6 50,2
7,20 51,0 50,2
5,30 37,6 50,1
3,20 22,7 50,0

Propriedades físicas das partículas: densidade 2,6 g/cm3, diâmetro médio de peneira 0,18 mm,
esfericidade 0,60. Massa de partículas sólidas, 5380g.

Fluidização com ar a 20ºC e 1 atm. Pede-se:

109
 a) Determinar através dos dados experimentais a porosidade e a velocidade na
fluidização mínima;

b) Verificar o resultado clássico da fluidização;

− ∆p = W / A

onde ∆p é a queda de pressão no leito, W é peso aparente da fase sólida e A a área de seção de
fluidização;

c) Estimar o valor da velocidade na fluidização mínima a partir da equação válida para

o escoamento darcyano de fluido,

2
φ 2 ε 3fm D p (ρ S − ρ F ) g
q fm = ⋅
170(1 − ε fm ) µ

onde φ é a esfericidade das partículas sólidas, ε fm a porosidade na fluidização mínima e D p

diâmetro médio das partículas.

Resposta:

a) Porosidade na fluidização mínima: 0,475. Velocidade na fluidização mínima: 2,02

cm/s.

b) Valor estimado para a queda de pressão na fluidização: 67,2 cm de coluna de água.

Através dos dados experimentais: 68,5 cm de coluna de água.

c) Valor estimado para a velocidade na fluidização mínima: 2,00 cm/s. Através dos
dados experimentais: 2,02 cm/s
2. Sobreiro (“Um Estudo de Fluidização a Altas Pressões”, Tese de M.Sc., COPPE/UFRJ,
1980) estudou experimentalmente a influência da pressão na fluidização de partículas
esféricas de vidro com ar a 20ºC:

Pressão Porosidade na Velocidade na Fluidização

(atm) Fluidização Mínima Mínima (cm/s)
1 0,502 0,147
5 0,491 0,143
10 0,483 0,146
15 0,483 0,147
20 0,480 0,147
25 0,476 0,145
30 0,476 0,145
35 0,472 0,146
Sabendo-se que o diâmetro médio das partículas é 30,4 µm, estimar pela equação apresentada
no problema anterior os valores da velocidade de fluidização mínima e comparar com os
resultados experimentais. A densidade das partículas de vidro é 2,43 g/cm3.

110
Resposta:

Porosidade na Velocidade na Fluidi- Densidade Velocidade na Fluidi-

Pressão Fluidização zação Mínima (exp.) do ar zação Mínima (cm/s)
(atm) Mínima (cm/s) (g/cm3) (estimada)
(exp.)
1 0,502 0,147 1,2×10−3 0,183
5 0,491 0,143 6×10−3 0,167
10 0,483 0,146 12×10−3 0,156
15 0,483 0,147 18×10−3 0,156
20 0,480 0,147 24×10−3 0,152
25 0,476 0,145 30×10−3 0,146
30 0,476 0,145 36×10−3 0,146
35 0,472 0,146 42×10−3 0,141

3. Deseja-se projetar um sistema de fluidização destinado à secagem de produto químico.

Diâmetro do secador: 30 cm.

Carga de sólido: 39 kg.
Propriedades das partículas: diâmetro médio 90 µm, esfericidade 0,8 e densidade 2,1 g/cm3.
Estimativa do valor da porosidade na fluidização mínima: 0,48.

Para uma velocidade superficial de ar 2 vezes maior que a de fluidização mínima, estimar:

a) A altura do distribuidor formado por esferas de aço com diâmetro 200 µm tal que a
queda de pressão através deste seja 10% da queda de pressão no leito fluidizado; porosidade,
0,38.

b) A potência do soprador para o serviço.

As propriedades do ar devem ser calculadas a 150ºC e 1 atm.

Resposta:

Estimativa do valor da velocidade de fluidização mínima, admitindo o escoamento como

sendo darcyano: 0,56 cm/s.
Cálculo da queda de pressão no leito fluidizado: 54,1 cm de carga de água.
Cálculo da altura do distribuidor: 7,1 cm.
Potência do soprador muito baixa, inferior a 0,1 cv.

4. Estimar o valor da porosidade no transporte hidráulico vertical de dolomita, fluxos

concorrentes: a) ascendente, e b) descendente.
Vazão mássica de água e de sólido: 116 ton/h e 110 ton/h.
Características das partículas: diâmetro médio 254 µm, esfericidade 0,7 e densidade 2,7
g/cm3.

111
Diâmetro da tubulação: 5".
Temperatura do fluido: 30ºC.
Resposta:

Cálculo da velocidade terminal da partícula isolada: cD Re∞2 = 504 , Re∞ = 9,67 ,

v∞ = 3,24cm / s .

Porosidade no transporte hidráulico, empregando a expressão de Richardson e Zaki:

- Ascendente,

254 89, 3
− = 3, 46 ε 2,50 → ε = 0, 739 .
ε 1− ε

- Descendente,

89, 3 254
− = 3, 46 ε 2,50 → ε = 0, 741.
1− ε ε

Pode-se notar que a velocidade relativa fluido-partícula é, neste caso, praticamente nula pois

89, 3 254
− = 0 → ε = 0, 740 .
1− ε ε

5. Seja o transporte pneumático vertical ascendente de alumina em tubo liso de 1,27 cm de

diâmetro interno. Calcular o gradiente de pressão no transporte sabendo que a vazão mássica
das fases fluida e sólida é de respectivamente 0,0514 g/s e 8,42 g/s. O transporte ocorre em
fase densa com porosidade da ordem daquela de fluidização mínima, no caso 0,48 (Santana,
1982).

Propriedades das partículas sólidas: densidade 3,97 g/cm3, diâmetro médio 0,20 mm e
esfericidade 0,7.

O gás de arraste tem as propriedades do ar a 20ºC e 3,3 atm.

Resposta:

Gradiente de pressão no transporte pneumático, equação (25): 2,03×103 dyn/cm3.

6. Estimar a faixa de velocidades superficiais do fluido (solução polimérica em xileno) entre

a fluidização mínima e o arraste de partículas.

Propriedades do fluido: densidade 0,87 g/cm3, e relação tensão cisalhante-taxa de distensão

dada por

112
 1000λ
S= dyn/cms2, com λ em s−1,
λ + 442

determinada experimentalmente para λ < 120s −1.

Propriedades das partículas: densidade 2,7 g/cm3, diâmetro volumétrico 0,32 mm e
esfericidade 0,78.
Porosidade na fluidização mínima: 0,47.

Resposta:

Velocidade do fluido na fluidização mínima, escoamento darcyano: 2,05 cm/h, com taxa de
distensão característica 0,491s−1 e viscosidade efetiva 2,26 P.
Velocidade terminal da partícula isolada, regime de Stokes: 148 cm/h, com taxa de distensão
característica 0,75s-1 e viscosidade efetiva 2,26 P.
Conclusão: A velocidade do fluido pode variar entre 2,05 e 148 cm/h.

7. A coluna de resina troca-iônica é lavada por meio de uma corrente ascendente de água que
acarreta uma expansão do leito e o conseqüente arraste das impurezas retidas. Estimar o valor
da velocidade superficial do fluido tal que a altura do leito expandido seja o dobro daquela do
leito fixo. A resina é constituída por partículas esféricas com diâmetro 0,3 mm e densidade
1,12 g/cm3. A porosidade do leito na fluidização mínima é estimada em 44%. A lavagem é
feita a 25ºC.

Resposta:

Cálculo da porosidade do leito na lavagem:

M = (1 − ε fm )ρ s AH fm = (1 − ε )ρ s A( 2 H fm ) → ε = 0, 72 .

Velocidade superficial do líquido na lavagem, usando a Forma Quadrática de Forchheimer:

9,13×10−2 cm/s ( k = 2,52 × 10 −5 cm 2 e c = 0,23).

8. Uma suspensão de minério finamente dividido em água tem o seguinte comportamento

reológico:
Taxa de distensão ( s −1 ) 0 3 10 50 100 200 300 600
2
Tensão cisalhante (g/cms ) 0 1,54 4,90 20,4 33,2 53,9 71,9 116.

Calcular o tempo necessário para carregar com a suspensão um caminhão com 10 m3 de

capacidade. A tubulação é de aço comercial e tem 2" de diâmetro (#40) e comprimento
equivalente total 25 m.
Densidade da suspensão: 1,3 g/cm3.

113
Resposta:

Formulação (equações 30 a 33):

Fator de atrito, f = 2 , 88 × 10−2 ; velocidade da mistura, VM = 429cm / s ; taxa de distensão

caracaterística, λ* = 528s −1 ; S = 1,32λ0,70dyn / cm(λ > 50s −1) ; viscosidade efetiva,
µ ef = 0, 203 P ; vazão de suspensão, Q M = 0, 56m3 / min.

Tempo para carregar o caminhão: 18 min.

9. Seja o transporte hidráulico de dolomita, 65/100 # Tyler, densidade 2,8 g/cm3 e

esfericidade das partículas 0,59. O transporte é feito a 30ºC em tubulação de aço, 4" de
diâmetro: 1500 m na horizontal e 150 m na vertical ascendente. A perda de carga nos
acidentes pode ser estimada em 20% das perdas nos dutos. Vazão mássica de sólido, 8 ton/h.

Calcular:
a) A vazão de água sabendo que a velocidade da mistura deve ser 20% superior àquela de
deposição das partículas;
b) A potência da bomba para o serviço.

Transporte vertical (equações 25 a 29).

Transporte horizontal (Santana, 1982):

Velocidade crítica da mistura, abaixo da qual ocorre o depósito de partículas,

0, 46 0,077
ρ   Dp 
VMC = 6,34 c1V/ 3 gD  S − 1 
 D 
 (34)
 ρF   

Gradiente de pressão,

114
  ∆p   ∆p 
−  − −  2  −3 / 2  D 0, 23  1,38
 L T  L  F  VM p ρ 
= 385  
 D 
  S − 1 (35)
 ∆p   gD     ρ 
cV  −    F
 L F

onde cV é a concentração volumétrica de sólidos.

Resposta:

Velocidade crítica da mistura, V MC = 189cm / s ; velocidade da mistura no transporte,

V M = 227cm / s ; vazão de água, QF = 63, 4m3 / h ; vazão de mistura Q M = 66, 2m3 / h ; carga
da bomba, 300 m de coluna de suspensão com densidade ρ M = 1, 08 g / cm3 ; potência da
bomba (eficiência 0,7), 115 cV.

10. Calcular a vazão de água e a potência de bombeamento requeridas para o transporte

hidráulico de 40 ton/h de areia na instalação abaixo esquematizada. Os dutos são de aço com
diâmetro 5" e o sistema deve operar com uma velocidade de mistura 20% maior que a
velocidade crítica de deposição. A perda de carga nos acidentes pode ser estimada em 25%
daquela proporcionada pelos dutos. Temperatura no bombeamento: 30ºC. Densidade e
esfericidade da areia: 2,6 g/cm3 e 0,78.

Distribuição granulométrica da areia:

# Tyler Fração Retida

−35+48 0,30
−48+65 0,40
−65+100 0,30

Resposta:

A solução deste problema é obtida através da formulação indicada no problema anterior.

Conclusões: vazão de água, 126 m3/h; vazão da mistura areia-água, 151 m3/h; concentração
volumétrica de sólido no transporte, 10,2%; potência da bomba, 50 cv (eficiência 0,6).

115
11. Especificar o diâmetro da tubulação e estimar a potência de bombeamento no transporte
de 80 m3/h de uma suspensão 2% em massa de polpa de papel em água. A instalação deve ter
2,5 km de dutos e a descarga encontra-se a 30 m acima do nível de alimentação.

Propriedades da suspensão: densidade, 1,02 g/cm3; dados reológicos.

Taxa de distensão ( s −1 ) 0 5 10 50 100 300 500

2
Tensão cisalhante (dyn/cm ) 0 1,39 2,65 10,2 16,0 32,6 45,4.

A velocidade de carga nos acidentes pode ser estimada em 15% da perda de carga nos dutos.

Resposta:

Diâmetro da tubulação: 5", o que leva à velocidade de mistura de 175 cm/s.

Cálculo da taxa de distensão característica e da viscosidade efetiva: respectivamente,
87,1 s−1 e 0,168 P (S = 0,8λ0,65 dyn/cm2, 50<λ<500 s−1).

Potência da bomba: carga da bomba, 130 m de suspensão; potência, 60 cv (eficiência 0,65).

12. Estimar a potência de bombeamento no transporte hidráulico de 12 milhões de

toneladas/ano de minério de ferro. A polpa tem 65% de minério, em massa.

Densidade do minério de ferro: 4,8 g/cm3.

Tubulação: 400 km de dutos de aço, 48 cm de diâmetro interno; a descarga está 1000 m
abaixo da alimentação. A instalação funciona 340 dias/ano.
Dados reológicos da suspensão:

Taxa de distensão ( s −1 ) 0 3 7 10 30 50 100 200

Tensão cisalhante (dyn/cm2) 0 0,151 0,348 0,488 1,33 2,11 3,96 7,45

(Dados semelhantes aos do mineroduto da Samarco que opera entre a Mina de Germano e o
Terminal de Ponta de Ubú).

Resposta:

Velocidade da mistura: 169 cm/s.

Cálculo da taxa de distensão característica e da viscosidade efetiva: respectivamente,
22,3 s−1 e 4,54×10−2 P (S = 0,060 λ0,91 dyn/cms2, 5 < λ < 200s −1 ).

Potência da bomba: carga da bomba, 674 m de carga de suspensão com densidade 2,07 g/cm3;
potência, 7600 cv (eficiência 0,75).

116
 Bibliografia

Coelho, G.L.V., Santana, C.C. e Massarani, G., "Reologia de Suspensões de Minério de

Ferro", Anais do IX ENEMP, Salvador, BA, vol.2, 27-38 (1981).

Ferreira, M.C., Silva, E.M.V. e Freire, J.T., "Mean Voidage Measurements and Fluid
Dynamics Analysis in a Pneumatic Bed with a Spouted Bed Type Solid Feeding System",
Brazilian Journal of Chemical Engineering, vol. 13, nº 01, 29-39 (1996).

Gidaspow, D., "Multiphase Flow and Fluidization: Continuum and Kinetic Theory
Descriptions", Academic Press, Boston, 467 p. (1994).

Gubulin, J.C. e Massarani, G., "Emprego da Radiação Gama em Sistemas Particulados",

2º Congresso Brasileiro de Eng. Química, S. Paulo (1977).

Kaus, M.N., "Pneumatic Conveying of Bulk Materials", McGraw-Hill, N.Iorque, 2º edição,

371 p. (1980).

Massarani, G., "Efeito da Aceleração no Escoamento Através de Meios Porosos", Anais do II

ENEMP, Rio Claro, SP, 47-55 (1974).

Massarani, G. e d'Ávila, J.S., "Transporte Vertical de Partículas Sólidas", Anais do II

ENEMP, Rio Claro, S.P., 33-47 (1974).

Massarani, G., "Aspectos da Fluidodinâmica em Meios Porosos", Revista Brasileira de

Engenharia, Número Especial, 96 p. (1989).

Massarani, G. e Silva Telles, A., "Escoamento de Soluções e Suspensões Não Newtonianas

em Dutos", Anais do XX ENEMP, S. Carlos, S.P. 27-33 (1992).

SAMARCO, "O Maior Mineroduto do Mundo", Construção Pesada, ano 7, julho, 19-
22(1977).

Santana, C.C. e Massarani, G., "Força Resistiva no Escoamento em Meios Porosos de Alta
Porosidade", Unimar, vol. 1, 49-58 (1974).

Santana, C.C., "Transporte Hidráulico e Pneumático de Partículas", em "Tópicos Especiais de

Sistemas Particulados" (Freire, J.T. e Jubulin J.C. editores), UFSCar, S.Paulo, 217-278
(1982).
Restini, C.V., "Transporte Vertical de Partículas Sólidas II: Análise Experimental", Tese de
M. Sc., PEQ/COPPE, Programa de Engenharia Química, COPPE/UFRJ, 57 p. (1977).

Wasp, E.J., "Slurry Pipelines", Scientific American, novembro, 48-55 (1983).

117
 Capítulo 5
Escoamento em Meios Porosos Deformáveis

1. Equações da Continuidade e do Movimento

As operações de filtração e espessamento de suspensões levam à formação de tortas e

de sedimentos que se caracterizam por exibirem uma variação de porosidade ao longo de sua
estrutura, causada pela percolação de fluido.

Tal como nos capítulos anteriores, a formulação para a fluidodinâmica em meios

deformáveis pode ser estabelecida a partir das equações da continuidade e do movimento para
cada fase e mais as equações constitutivas do sistema (Silva Telles e Massarani; 1989;
Massarani et al., 1993):

Equações da continuidade

∂
( ερ F ) + div( ερ F v F ) = 0 (1)
∂t

∂
[(1 − ε )ρ F ] + div[(1 − ε )ρ S vS ] = 0 (2)
∂t

Equações do movimento

 ∂v 
ερF  F + (grad vF )vF  = −div p − m + ρF g (3)
 ∂t 

 ∂v 
(1− ε )ρ S  S + (grad vS )vS  = grad TS + m + (1 − ε )(ρ S − ρ F ) g (4)
 ∂t 

Equações constitutivas

µF
m= εU , U = v F − v S (5)
k

TS = T (FS ) . (6)

Nestas equações, ε é a porosidade em um ponto da matriz (fração volumétrica ocupada pelo

fluido), ρ F e ρ S a densidade do fluido e do sólido, vF e vS a velocidade intersticial das fases
fluida e sólida, p a pressão no fluido, m a força resistiva que o fluido exerce sobre a matriz
sólida (por unidade de volume de sistema particulado), g a intensidade do campo exterior, TS
a tensão exercida sobre a fase sólida e FS o gradiente de deformação.

118
 A equação da continuidade para a fase sólida pode ser escrita em termos do gradiente
de deformação

(1 − ε*)ρ*S = (1 − ε )ρ S det(FS ), (7)

onde ε* e ρ*S denotam, respectivamente, a porosidade e a densidade na configuração de

referência. Sendo a densidade do sólido constante,

1− ε *
det FS = . (8)
1− ε

A força resistiva expressa pela equação (5) é válida nas condições em que o
escoamento é lento, situação que prevalece tanto na filtração quanto no espessamento. Nesta
equação, µF é a viscosidade do fluido, k(ε) a permeabilidade na matriz porosa e U a
velocidade relativa fluido-partícula.

A forma constitutiva expressa pela equação (6) considera que a matriz porosa
comporta-se como um material elástico não-linear. Nos casos comuns de filtração,
espessamento ou adensamento pode-se considerar que o meio poroso se submeta a uma
deformação plana segundo um dada direção, por exemplo ao longo do eixo-z,

x = X

y = Y (9)
z = Z + γ ( Z , t ),


sendo x, y e z a posição de uma partícula que na configuração de referência ocupava a posição

X, Y e Z. Resulta das equações (8) e (9), admitindo que o meio seja inicialmente homogêneo,

 
1 0 0 

FS ( z , t ) = 0 1 0  . (10)
 1− ε * 
0 0 
 1 − ε( z, t ) 

Portanto, os componentes de Ts na deformação plana são:

Txx = Tyy ≠ Tzz

 1 − ε *
Tzz = f   (11)
 1− ε 
Tij = 0 , i ≠ j .

A função f da equação (11) define a pressão sobre os sólidos pS

f = − pS ( ε ) (12)

119
e, como conseqüência, a equação do movimento para a fase sólida pode ser escrita da seguinte
maneira quando os efeitos de aceleração não são considerados e a deformação do meio é
plana:

d µ
0=− pS ( ε ) + F εU z + (1 − ε )(ρ S − ρ F ) g z . (13)
dz k

Os estudos relacionados às aplicações clássicas envolvendo a filtração com formação

de torta, a expressão mecânica e a sedimentação contínua defrontam-se com dois grandes
desafios: de um lado, a dificuldade experimental no levantamento do perfil de porosidades e
da dependência entre a porosidade, pressão nos sólidos e permeabilidade, de outro lado, a
solução de um requintado sistema de equações diferenciais não lineares ao qual estão
associados condições de salto e contornos móveis.

Do ponto de vista tecnológico, o projeto do equipamento de separação industrial não

pode ainda prescindir de ensaios de bancada conduzidos diretamente com a suspensão a ser
tratada, e o “scale-up”, estabelecido através de teorias simplificadas, é basicamente uma
regra-de-três entre capacidade e área de separação.

Exemplo

Damasceno (1992) desenvolveu uma técnica para a determinação das relações entre
porosidade, pressão nos sólidos e permeabilidade nas condições que prevalecem no
espessamento contínuo. Nesta técnica, o sedimento formado no interior de um frasco de
Mariotte se submete a um processo de deformação causado pela percolação lenta do líquido.
Uma vez conhecido o perfil de porosidades (método da atenuação de raios gama) é possível
calcular pS ( ε ) e k ( ε ) a partir da equação (4), em duas etapas,

1ª etapa:

Uz = 0
ε = ε ( z ) (experimental)
água L
pS ( z ) = (ρ S − ρ F ) g ∫ [1 − ε ( z )]dz (14)
z
pS = pS ( ε );
sedimento
2ª etapa:
L
z − εU z = QF / A
placa porosa
ε = ε ( z ) (experimental)
(15)
µ F (QF / A)
k=−
dpS dε
Frasco de ⋅ + (1 − ε )(ρ S − ρ F ) g
x Mariotte dε dz
k = k ( ε ).

120
Nestas equações, QF é a vazão de líquido que percola o sedimento e A é a área da seção
transversal de percolação. Considera-se na equação (14), por falta de melhor informação, que
a pressão nos sólidos na interface líquido-suspensão seja nula, isto é, pS ( L) = 0 .

Estão reunidos na figura (1) alguns resultados obtidos com uma suspensão aquosa de
carbonato de cálcio (diâmetro médio das partículas da ordem de 18µm).

0,25
1-ε

0,20

0,15

0,10

0,05 -4
QF/A = 1,25x10 cm/s

0
0 10 20 30 Z (cm)

1500 10
-8
Ps (Pa)

1-ε )
[( 0,08 8,4
]
ps = 7,0 -1 Pa
k = 1,15x 10
-8 1-ε )
( 0,08 -4,1
cm
2
k(cm2)

1250

1000

-9
750 10

500

250
-10
0 10
0 0,05 0,10 0,15 0,20 0,30 0,1 0,14 0,18 0,22
1-ε 1-ε

Figura 1 - Suspensão aquosa de CaCO3: relações entre porosidade, pressão nos sólidos e
permeabilidade.

121
2. Teoria da Filtração com Formação de Torta

A denominada “teoria científica da filtração” foi desenvolvida nos últimos 40 anos

pelas escolas ilustres de Frank M. Tiller (Universidade de Houston, Estados Unidos) e
Mompei Shirato (Universidade de Nagoya, Japão). Esta teoria leva à “teoria simplificada”,
base para o projeto e análise do desempenho de filtros.

Seja a filtração plana como esquematizada na figura (2)

P(-lm, t) p(0,t) p(l,t)

meio torta
filtrante

suspensão
filtrado
sólido-líquido

lm l(t)
z

Figura 2 - Filtração plana com formação de torta

Considerando que o líquido e as partículas sólidas sejam incompressíveis, a

combinação das equações (1) e (2) da continuidade leva ao resultado

∂
(q F + qS ) = 0
∂z (16)
q F = εv F , q S = (1 − ε )vS

que permite, por integração, relacionar as velocidades das fases líquida e sólida,

q S ( z , t ) = q F (0, t ) − q F ( z , t ) , (17)

sendo qS (0, t ) = 0 na interface torta-meio filtrante.

Sendo a velocidade superficial da fase sólida positiva em qualquer ponto da torta

compressível, a equação (17) informa que a velocidade superficial do fluido é maior junto ao
meio filtrante do que na interface suspensão torta,

q F (0, t ) > qF(l,t),

122
resultado que, apesar de correto, suscitou discussões na literatura quando apresentado pela
primeira vez (Tiller, 1961).

A combinação das equações (3) e (13) do movimento leva à relação entre as pressões
nas fases líquida e sólida. Desconsiderando a aceleração das fases e a influência do campo
exterior,

∂
( p + pS ) = 0
∂z (18)
pS (z, t ) = p(l , t ) − p(z, t )

sendo pS ( l , t ) = 0 na interface torta-suspensão. Portanto, dentro das hipóteses consideradas,

a pressão no sólido na filtração resulta apenas da queda de pressão no líquido.

Equacionamento da filtração plana com formação de torta

• Equação da continuidade para a fase líquida:

∂ε ∂q F
= (19)
∂t ∂z

• Equação do movimento para a fase líquida e equações (5) e (17):

∂p µ F ∂p µ  q q (0, t ) − q F 
0=− + ε(v F − vS ) = − + F ε F − F  (20)
∂z k ∂z k  ε 1− ε 

sendo a viscosidade do fluido conhecida.

• Equação que correlaciona as pressões nas fases:

pS (z, t ) = p(l , t ) − p( x , t ) . (18)

• Equações constitutivas referentes ao sistema particulado:

pS = pS ( ε ) (21)

k = k (ε) . (22)

• Relação entre a espessura da torta e o tempo de filtração:

l
ρS ∫ (1 − ε)dz
C= 0 , (23)
t l
ρF ∫ q F (0, t )dt + ρF ∫ εdz
0 0

onde C é a concentração de sólidos na suspensão a ser filtrada,

123
 ~
massa de sólido na suspensao massa de sólido na torta
C= ~ = .
massa de líquido na suspensao massa de filtrado + massa de líquido na torta

• Equação do movimento do fluido no meio filtrante:

p(0, t ) − p(− l m, t ) µF
= q F (0, t ) (24)
lm km

sendo o valor do comprimento e da permeabilidade do meio filtrante conhecidos.

• Relação vazão-queda de pressão fornecida pelo sistema de bombeamento de suspensão:

q F (l , t ) = f [p( l , t ) − p(−l m, t )] . (25)

• Condições limites:

l (0) = 0
p(l , t ) = f1( t ) ou p(−l m, t ) = f 2 ( t ).

Exemplo

Silva Telles et al. (1973), fazendo uso da mesma formulação aqui apresentada,
estudaram a operação de filtração de uma suspensão aquosa de caulim. A técnica utilizada
nesta simulação é a de similaridade, resultando uma equação diferencial ordinária não linear
para a pressão nos sólidos, resolvida por técnica numérica. Não é considera na análise a
resistência oferecida pelo meio filtrante.

Concentração da suspensão: 333g de sólido/L de água.

Condições operacionais: a filtração é conduzida a 30ºC, sob a queda de pressão

constante com valor 4 atm.

Propriedades da torta (Shirato et al., 1964):

ε = 0,72 , pS ≤ 104 dyn / cm2



ε = 1 − 0,131 p 0,082 , p ≥ 104 dyn / cm2
 S S

k = 4,2 × 10−11 cm2 , pS ≤ 104 dyn / cm2



k = 2,21 × 10−9 p −0,43cm2 , pS ≥ 104 dyn / cm2 .
 S

Os resultados apresentados nas figuras (3) fornecem as distribuições de pressão nos

sólidos, porosidade e permeabilidade ao longo da torta compressível formada na operação de

124
 filtração. As figuras mostram também o resultado óbvio de que as maiores taxas de produção
de filtrado são alcançadas no inicio do processo, quando a torta ainda é delgada. A estratégia
da operação com formação de tortas delgadas é utilizada, por exemplo, na filtração industrial
conduzida no filtro rotativo a vácuo.

1 -11
ps(z,t)/p s(0,t)

torta incompressível 4x10

k(cm )
2
0,8 0,75
-11
3x10
torta incompressível
0,6 0,70
-11
torta incompressível 2x10
0,4 0,65

-11
1x10
0,2 0,60

0,55 0
0 0,2 0,4 0,6 0,8 1 0 0,2 0,4 0,6 0,8 1
z /l
z/l

0,5 0,5

0,4 0,4
Volume de filtrado 3 2
por unidade de área v, cm/cm

Espessura da torta l,cm

0,3 0,3

0,2 0,2

0,1 0,1

0 0
0 2 4 6 8 10
Tempo de filtração t,min

Figura 3 - Filtração de suspensão aquosa de caulim ( ∆p = 4atm) .

125
A teoria simplificada da filtração

O ponto de partida para o estabelecimento da “teoria simplificada da filtração” é

admitir que a velocidade superficial do sólido na torta seja substancialmente menor do que
aquela do líquido, resultando da equação (17),

q F = q F (t ) , (26)

isto é, que a velocidade do líquido é apenas função do tempo de filtração. Esta hipótese,
aceitável quando a torta é moderadamente compressível, permite reduzir a complexidade da
solução do problema.

A equação do movimento para o fluido na torta toma a seguinte forma a partir das
equações (18), (20) e (26):

∂ ∂ µ
p( z , t ) = − pS ( z , t ) = F q F (t ). (27)
∂z ∂z k

Na próxima etapa procura-se correlacionar os resultados da filtração com as condições

operacionais expressas pela queda de pressão no filtro. Seja M a massa de sólido seco que
compõe a torta,

dM = ρ S (1 − ε ) A dz ,

onde A é área da superfície de filtração. Vem da equação (27),

1 µ q µ q
− dpS = ⋅ F F dM = α(ε) F F dM ,
ρS (1 − ε)k A A

sendo α a resistividade local (L/M),

1
α= .
ρ S (1 − ε ) k ( ε )

Integrando:

0 dpS p(l , t ) − p(0, t ) µ Fq F

−∫ = = M, (28)
p (l , t ) − p (0, t ) α <α> A

p(l , t ) − p(0, t )
< α >= , a resistividade média relativa à queda de pressão p(l , t ) − p(0, t ) .
p(l , t ) − p(0, t ) dpS
∫0 α

A torta e o meio filtrante são meios porosos percolados em série pelo fluido. A
expressão para a queda de pressão no filtro pode ser estabelecida combinando as equações
(24) e (28),

126
  M 
p(l , t ) − p(−l m, t ) = ∆p = < α > + R m  µ Fq F , (29)
 A 

R m =l m / k m , a resistência do meio filtrante (1/L).

A velocidade superficial do fluido e a massa de sólido seco na torta estão relacionadas

ao volume de filtrado V, ao tempo t, à área de filtração A e à concentração de sólidos na
suspensão, C:

1 dV
qF = (30)
A dt

~ m
C− . (31)
ρ FV

A equação da filtração resulta da combinação das equações (29) a (31),

dt µ F  < α > Vρ F C 
=  + Rm , < α >= f ( ∆ P ) . (32)
dV A( ∆p)  A 

Na maioria das situações de interesse industrial a filtração é conduzida sob queda de

pressão constante:

t µ F  < α > Vρ F C 
=  + Rm  . (33)
V A( ∆p)  2A 

Exemplo

A equação (33) permite calcular a resistividade média da torta e a resistência do meio

filtrante a partir das medidas de volume de filtrado e tempo de filtração obtidas na operação
sob queda de pressão constante.

A figura (4) refere-se à filtração de uma supensão aquosa de talco conduzida em

unidade de bancada do filtro COPPE (Massarani, 1985). Confirmando um resultado clássico,
a resistividade média da torta aumenta com a queda de pressão e a resistência do meio
filtrante é praticamente constante:

< α >= 2,17 × 1010 ( ∆p)1,05 cm / g , ∆p em atm

Rm = 4,11 × 109 cm−1.

127
t/v (s/l)

11
∆p=5atm
10
∆p=8atm
9
∆p=11atm
8

4
0 5 10 15 V(l) 20
Figura 4 - Desempenho do filtro COPPE: suspensão aquosa de talco, A = 670cm2
(Massarani, 1985).

3. A Sedimentação Contínua

Do ponto de vista da separação sólido-líquido, o projeto do sedimentador contínuo

está basicamente relacionado ao cálculo da área da seção de sedimentação e da altura do
equipamento.

Tal como na filtração, face às dificuldades e incertezas no estabelecimento das

equações constitutivas para as lamas compressíveis, o projeto do sedimentador acaba se
baseando nos ensaios em batelada conduzidos diretamente com a suspensão a ser tratada.
Neste sentido, o procedimento clássico desenvolvido por Kynch leva a resultados satisfatórios
quando a lama é moderadamente compressível (Damasceno, 1992; Damasceno e Massarani,
1993). A conexão entre o ensaio de proveta e o projeto de sedimentadores contínuos é
abordada nos problemas (12) a (15) reunidos no final do capítulo.

A utilização de unidade piloto contínua pode ser de grande valia na determinação

direta da capacidade máxima de sedimentação do sistema em estudo. Avalia-se, neste caso, a
relação entre a vazão de alimentação e a altura da região de compactação que se forma no
fundo do
sedimentador. A figura (5), referente à sedimentação de suspensão aquosa de carbonato de
cálcio, mostra que a capacidade de sedimentação pode ser aumentada quando o ponto de
alimentação é deslocado da supefície do sedimentador (configuração tradicional) para o fundo
do equipamento (França, 1996). O aumento registrado na capacidade de sedimentação é de
0,39 a 0,52 m3/m2h.

128
 Altura da região de compactação (cm)

140 Alimentação no topo

Alimentação a 12 cm do fundo
120

100

80

60

40

20

0
0,15 0,20 0,25 0,30 0,35 0,40 0,45 0,50
3 2
Vazão de alimentação (m /m h)

Figura 5 - Sedimentação contínua de suspensão aquosa de carbonato de cálcio.

Concentração volumétrica de sólido na alimentação e no lodo: 1,3 e 6,7% (França, 1996).

Cabe mencionar que a relação entre a altura da região de compactação e a capacidade

do sedimentador pode ser estabelecida pela combinação das equações da continuidade e do
movimento, equações (1), (2) e (13), uma vez conhecidas as equações constitutivas para a
região mencionada (Tiller e Chen, 1988):

pS [ ε ( 0)] dpS
L=∫ , (34)
0 µq S  1 1 
(ρ S − ρ F )(1 − ε ) g − −
k ( ε ) 1 − ε 1 − ε (0) 

conhecendo-se

ε (0), pS = pS ( ε ) e k = k ( ε ) .

129
 alimentação

extravasante

z
L

lama

Nesta equação, qS é o fluxo volumétrico de sólido. A capacidade máxima do sedimentador

está associada ao comprimento infinito da região de compactação.

Do ponto de vista científico, os problemas relacionados à sedimentação contínua vão

muito além das dificuldades relativas à solução do sistema de equações diferenciais não
lineares estabelecido a partir das equações de conservação e das equações constitutivas para
as condições que prevalecem no sedimentador. A operação de protótipos contínuos parece
indicar, por exemplo, que o estado de equilíbrio atingido no regime permanente, com as
interfaces entre as diferentes regiões estáticas, depende do procedimento adotado na partida
da operação. Em algumas situações, o equilíbrio rompe-se sem causa aparente e a região de
compactação atinge rapidamente o extravasante que antes se apresentava límpido (França,
1996).

130
Problemas: Escoamento em Meios Porosos Deformáveis

1. Deseja-se filtrar uma suspensão aquosa de carbonato de cálcio (23,5 g de sólido/L de

água) com o auxílio da bomba centrífuga Worthington 1¼ FT214, 5 cv. Pede-se: o volume de
filtrado e a espessura da torta em função do tempo de filtração. A área filtrante é 1 m2.

a) Propriedades físicas do sistema: densidade do fluido, 1 g/cm3; densidade do carbonato de

cálcio, 2,7 g/cm3; viscosidade da água, 0,90 cP.
b) Propriedades da torta e do meio filtrante:

Queda de Pressão Porosidade Média Resistividade Média Resistência do

(atm) da Torta da Torta Meio Filtrante
(cm/g) (cm−1)
0,46 0,73 1,11x1010 8,52×108
1,10 0,72 1,47×1010 9,02×108
1,92 0,69 1,61x1010 10,3×108
2,47 0,67 1,74x1010 10,2×108
3,34 0,65 1,85×1010 10,4×108
4,60 0,60 2,23×1010 10,1×108
5,70 0,59 2,36×1010 10,3×108

c) Curva característica da bomba:

Q(m3/h) 2,0 3,6 8,4 11 13 15 20

C(m água) 68 67 66 63 60 55 35.

Resposta:

Curva característica da bomba (Q em m3/h):

C = 68m, Q ≤ 9 m3/h
C = 64,7+1,52Q - 0,149Q2 m de coluna de água Q > 9 m3/h. (35)

Propriedades da torta em função da queda de pressão de filtração:

< α >= 1, 38 × 1010 ( ∆p) 0,29 cm/g (∆p em atm), (36)

< ε >= 0, 744 − 0, 0292 ∆p (∆p em atm). (37)

Resistência do meio filtrante: 109 cm−1.

Formulação:

Resulta da equação (33) da filtração que

131
 A  A∆p 
V=  − µRm  , ∆p = ρ F gC . (38)
µ < α > ρF C  Q 

O tempo de filtração e a espessura da torta correspondentes ao volume V de filtrado são dados

por

V dV
t=∫ , (39)
0 Q

cρFV
l= . (40)
[1− < ε >]ρSA
Q C ∆p <α> <ε> 1/Q V t l
(m3/h) (eq.35) (atm) (eq.36) (eq.37) (h/m3) (eq.38) (eq.39) (eq.40)
(m) (cm/g) (l) (s) (mm)
17,7 44,9 4,49 2,13×1010 0,61 0,0565 0 0 0
15,0 54,0 5,40 2,25×1010 0,59 0,0667 7,8 1,7 0,2
10,0 65,0 6,50 2,37×1010 0,55 0,100 27,8 7,7 0,5
7,0 68,0 6,80 2,41×1010 0,55 0,143 49,6 17,3 1,0
5,0 68,0 6,80 2,41×1010 0,54 0,200 77,4 30,4* 1,5
3,0 68,0 6,80 2,41×1010 0,54 0,333 140 81,1* 2,7
1,0 68,0 6,80 2,41×1010 0,54 1,00 458 788* 8,7
0,3 68,0 6,80 2,41×1010 0,54 3,33 1570 9180* 29,7

A tabela mostra que após apenas 17,3 s, correspondendo ao volume de filtrado de 49,6L, a
filtração prossegue a queda de pressão constante.

Nesta situação, desprezando a resistência oferecida pelo meio filtrante,

t − t0 =
µ < α > ρF C
2
2 A ∆p
(V 2
)
− V02 ,

onde t0 = 17,3s e V0 = 49,6 L . Os valores do tempo de filtração assinalados na tabela foram

calculados através desta última equação.

2. Foram obtidos os seguintes resultados na filtração de uma suspensão aquosa de carbonato

de cálcio (50 g de sólido/L de água) em filtro-prensa piloto operando com apenas um quadro,
dimensões nominais 6×6×1¼ in (valores efetivos 15,1×15×3,2 cm), a 25ºC e queda de pressão
2,72 atm. Determinar a resistividade e porosidade médias da torta, <α> e <ε>, e a resistência
do meio filtrante Rm a partir dos dados reunidos na tabela que se segue. A densidade do
carbonato de cálcio é 2,7 g/cm3 e a relação entre massa de torta e massa de torta seca é 1,60.

132
 Tempo de Filtração Volume de Filtrado
(s) (L)
18,0 0,700
40,7 1,70
108 3,70
160 4,70
321 7,70
467 9,70
550 10,7
638 11,7
833 13,7
943 14,7
1084 15,7
1215 16,7
1425 17,7
1702 18,7
2344 19,7

Resposta:

A representação gráfica de t/V em função do volume de filtrado

120,0
t/V(s/L)

80,0

40,0 2
Área de filtração: 456 cm
∆p= 2,72 atm

0,0
0,0 4,0 8,0 12,0 16,0 20,0
V(L)

mostra que a teoria simplificada da filtração é válida e que:

A resistividade média da torta é 7,34×109 cm/g;

A resistência do meio filtrante é 2,78x109 cm−1;

O volume de filtrado e o tempo de filtração correspondentes ao quadro cheio (final da

reta) são respectivamente 14,5 L e 920s .

133
A porosidade média da torta pode ser calculada a partir da relação entre a massa de torta no
final da filtração e a massa de torta seca: 0,62. A relação entre os volumes de filtrado e de
torta, que depende da natureza do sistema, concentração de sólidos na alimentação e da
pressão de operação, pode ser calculada através do gráfico ou por meio da expressão para a
concentração de sólidos na alimentação:

C=
[1− < ε >]vt ρS =
[1− < ε >]ρS .
Vρ F + < ε > vt ρ F V
ρF + < ε > ρF
vt

Neste problema, V/vt é aproximadamente 20.

3. Foram obtidos os seguintes dados em filtro rotativo de laboratório com 3000 cm2 de
superfície filtrante, operando a uma queda de pressão de 0,73 atm.

Rpm do tambor Vazão de filtrado

(L/min)
1,17×10−2 0,370
5,00×10−2 0,719
12,0×10−2 0,897
36,7×10−2 1,30
57,0×10−2 1,47

Concentração de sólidos na alimentação: 55g de sólido/L de líquido.

Ângulo de imersão do tambor: 80º.
Relação entre massa de torta e massa de torta seca: 1,85.
Propriedades do filtrado: densidade 1 g/cm3 e viscosidade 1,1 cP.
Densidade do sólido: 3,1 g/cm3.

Determinar a resistividade e a porosidade médias da torta e a resistência do meio filtrante.

Resposta:

Resistividade média da torta: 1,97x1010 cm/g.

Resistência do meio filtrante: 1,37x109 cm−1.
Porosidade média da torta: 0,725.

4. Especificar o filtro-prensa com quadros de metal para a filtração de 10 m3/h da suspensão

do problema (2). Condições operacionais: 25ºC e 2,72 atm.

1º caso: a torta não requer lavagem.

2º caso: a lavagem deve ser efetuada com volume 2 vezes maior que o volume de torta.

Considerar nas duas situações que o tempo de desmantelamento, limpeza e montagem do

filtro seja de 20 min.

134
Dados provenientes do Catálogo 59 da T. Shriver & Company (Harrison, N.J., Estados
Unidos):

Dimensões recomendadas para placas e quadros

Área total de filtração Dimensão nominal

(ft2) dos elementos (in)
5−35 12
30−100 18
75−250 24
150−450 30
250−700 36
500−1100 43 ¼
>1000 48 e 56

Área filtrante efetiva por quadro

Dimensão nominal Área filtrante efetiva por quadro (ft2)

dos elementos (in) Metal Madeira
12 1,7 0,9
18 3,9 2,3
24 7,0 4,8
30 10,5 7,3
36 15,6 10,5
43¼ 22,2 15,1
48 28,8 19,7
56 - 28,4

Resposta:

O índice 1 denota a unidade de laboratório e 2 a unidade industrial.

Produção de filtrado na unidade industrial:

(V f )2
P=
(t f )2 + tl + t d . (40)

Tempo de filtração na unidade industrial (espessura do quadro, e):

t f = BV f2 / A2

135
 Vf  2Vf   2Vf  (Vf / A )2 e
=  =  , = 2 (42)
v t  Ae 1  Ae 2 (Vf / A )1 e1

 (V / A )2   e2 
2
(t f )2 =  f 2  (t f )1 =   (t f )1 . (43)
 (Vf / A )12   e1 

Tempo de lavagem da torta na unidade industrial aparelhada com placas de lavagem (placas
com "três botões"):

V 1  dV 
Q = l =  
l tl 4  dt  final da filtração

v 
t = 8β t  (t f )2 , β = V / v t . (44)
l  Vf l l

Volume de filtrado por ciclo completo na unidade industrial:

(Vf )2 = [(t f )2 + tl + t d ]P . (45)

Área de filtração da unidade industrial:

A2 = A1
(V f )2 e1 . (46)
(V f )1 e2
Resulta do problema (2) na filtração a ∆p = 2 , 72 atm e 25ºC: espessura do quadro e1 = 3, 2 cm,
área de filtração A1 = 456 cm2, tempo de filtração, (t f )1 = 920 s; volume de filtrado,
(Vf )1 = 14,5 L; relação entre volume de filtrado e volume de torta, (Vf / v t )1 = 20 .

Produção de filtrado na unidade industrial: P = 9820 L / h .

136
a) A torta não requer lavagem:

e2 (t f ) 2 , eq. (43) (t f ) 2 + td (V f ) 2 , eq. (45) A2 , eq. (46)

(in) (min) (min) (L) (m2)
1 9,79 29,8 4877 19,2
1¼ 15,3 35,3 5777 18,2
1½ 22,0 42,0 6874 18,0
1¾ 30,0 50,0 8183 18,4
2 39,2 59,2 9689 19,0
3 88,1 108 17676 23,2

Solução possível: área filtrante, 18m2; espessura dos quadros, 1½ in; dimensão nominal dos
elementos, 30 in; número de quadros, 19.

b) A torta requer lavagem

e2 ( t f ) 2 , eq. (43) tl , eq. (44) ( t f ) 2 + tl + t d (V f ) 2 , eq. (45) A2 , eq. (46)

(in) (min) (min) (min) (L) ( m2 )
1 9,79 7,83 37,6 6154 24,2
1¼ 15,3 12,2 47,5 7774 24,4
1½ 22,0 17,6 59,6 9755 25,6
1¾ 30,0 24,0 74,0 12111 27,2
2 39,2 31,4 90,6 14828 29,1
3 88,1 70,5 179 29296 38,4

Solução possível: área filtrante, 25,6m2; espessura dos quadros, 1½ in; dimensão nominal
dos elementos, 30 in; número de quadros, 27.

5. Especificar o filtro rotativo a vácuo a partir dos dados obtidos em filtro-folha de

laboratório com suspensão aquosa de carbonato de cálcio, 50 g de sólido/L de suspensão.
Densidade do carbonato de cálcio: 2,7 g/cm3. Queda de pressão no filtro: 600 mm Hg.
Temperatura de operação: 28ºC. Produção de filtrado: 10000 L/h.

Resultados obtidos no filtro-folha operando com a mesma suspensão, nas mesmas condições
operacionais indicadas e área filtrante 133 cm2:

Tempo de filtração para se obter uma torta com 6 mm de espessura (volume de filtrado
950 cm3), 163 s;
Tempo de lavagem da torta (volume de água de lavagem 160 cm3), 130 s;
Tempo de secagem (obtém-se um produto com 81% de sólido em massa), 150 s;
Tempo estimado para a descarga da torta e limpeza do meio filtrante, 10 s.

137
 Dimensões padrões de filtros a vácuo Dorr-Oliver

Diâmetro Área da superfície do filtro, ft2

do tambor, Comprimento, ft
ft 4 6 8 10 12 14 16 18 20 22 24
6 76 113 151 189 226
8 200 250 300 350 400
10 310 372 434 496 558 620
12 456 532 608 684 760 836 912
(Perry, H.R. e Chilton, C.I., "Manual de Engenharia Química", 5ª edição traduzida,
Guanabara Dois, Rio de Janeiro, p. 19-72, 1980)

Resposta:

Sendo o tempo de um ciclo completo 453 s, resulta que a rotação do tambor deve ser
0,132 rpm.

A fração submersa é 163/453 e, portanto, o ângulo de imersão é 130º.

Produção de filtrado por unidade de área filtrante: 570 (L/m2h).

Especificação do filtro considerando um fator de segurança de 20% no cálculo da área:

diâmetro do tambor, 6 ft; comprimento do tambor 12 ft.

6. Deseja-se filtrar uma suspensão altamente viscosa em filtro-folha, a pressão constante. A

adição do solvente MK-X acarreta uma diminuição na viscosidade do filtrado levando ao
valor

µ = 0,63F −4 ,48 P,

onde F é a fração volumétrica do solvente no filtrado. Filtrando em cada ciclo o volume

original de suspensão, determinar o valor de F que conduz ao ciclo mais curto. Desprezar a
resistência do meio filtrante e admitir que o solvente não afete a densidade do filtrado.

Resposta:

<α>ρ
Equação da filtração: t = BµcV 2 , B = .
2 A2 ∆p

Volume de filtrado por ciclo: V = V0 / (1 − F ) , onde V0 é o volume de filtrado isento de

solvente.

Concentração de sólido na suspensão a ser filtrada: C = (1 − F ) M / V0 , onde M é a massa de

sólido seco na torta.
Valor de F que conduz ao ciclo mais curto: 0,82.

138
7. A usina de beneficiamento de caulim Cellini de Mar de Espanha, MG, opera com uma
bateria de filtros-prensa constituída por 210 quadros de 30 in (área filtrante efetiva por quadro
10,5 ft2) e 2 in de espessura. A produção de torta é de 5,8 ton/h, sendo, para cada ciclo, o
tempo de filtração 55 min e o de desmantelamento, limpeza e montagem 20 min. O Sr.
Cellini deseja aumentar a produção em 30% e, ao mesmo tempo, melhorar a qualidade de seu
produto através da lavagem da torta com um volume de líquido de lavagem igual ao de torta.
Sabe-se ainda, através de experiências conduzidas em laboratório, que nas condições
operacionais a relação entre os volumes de filtrado e de torta é 9. Considerar que nas novas
instalações o tempo de desmantelamento, limpeza e montagem é também de 20 min e que a
operação do sistema aumentado se fará à mesma queda de pressão que no sistema inicial.

Resposta:

Usando o mesmo procedimento apresentado no problema (4), resulta que o número de

quadros adicionais é 242.

8. O ferro-velho de Maria da Graça dispõe de um filtro-prensa Shriver de metal, completo:

placas e quadros de 30", 20 quadros de 2" de espessura (área filtrante por quadro, 10,5 ft2).
Determinar a capacidade do filtro operando com uma suspensão aquosa de carbonato de bário
(70 g de sólido/L de água) a 30ºC e com uma queda de pressão de 65 psi. A lavagem da torta,
realizada nas mesmas condições que na filtração, deve empregar um volume de água de
lavagem 1,5 vez o volume de torta. O filtro está aparelhado com placas de “3 botões”. O
tempo de desmantelamento, limpeza e montagem é estimado em 20 min. Testes de
laboratório conduzidos a 30ºC e 65 psi em um único quadro com 1¼" de espessura e área
filtrante de 456 cm2 levaram aos seguintes resultados:

Tempo de filtração e volume de filtrado na situação de quadro cheio, 18 min e 14L;

Relação entre massa de torta e massa de torta seca, 1,5.

A densidade do carbonato de bário é 4,1 g/cm3.

Resposta:

Porosidade média da torta: 0,67.

Relação volume de filtrado-volume de torta: 19,2.
Volume de filtrado produzido em um ciclo completo da unidade industrial: 9,52×103 L.
Tempo de filtração por ciclo: 46,1 min.
Tempo de lavagem da torta: 28,8 min.
Produção de filtrado: 6020 L/h.
Produção de sólido seco: 421 kg/h.

9. Deseja-se filtrar uma suspensão em filtro de tambor rotativo a vácuo de 2 m de diâmetro e

2 m de comprimento. A suspensão apresenta 18% em massa de sólido e as propriedades
físicas do sistema são: densidade do fluido 1 g/cm3, viscosidade do fluido 0,8 cP e densidade
do sólido 2,3 g/cm3.

139
O filtro industrial deve operar com uma queda de pressão de 25 in Hg, com um ângulo de
submersão de 145º. A espessura final da torta deve ser de 5 mm. Pede-se:

a) O número de rotações por minuto do cilindro;

b) A produção de filtrado.

Ensaios de laboratório conduzidos com esta mesma suspensão em filtro-folha de 200 cm2 de
área filtrante, queda de pressão de 25 in Hg, levaram a tortas com 63% de porosidade e à
seguinte relação entre tempo de filtração e volume de filtrado:

Tempo de filtração (min): 1 2

3
Volume de filtrado (cm ): 300 440.

Resposta:

Volume de filtrado por ciclo: 297L.

Tempo de filtração por ciclo: 345s.
Rpm do tambor: 0,174.
Produção de filtrado: 3100 L/h.

10. Calcular a área do filtro gravitacional para o tratamento de uma suspensão aquosa de
rejeito industrial de pequena fábrica da Baixada Fluminense. O filtro tem 1m de altura.

Critério: O tempo de filtração deve conduzir à produção máxima,

Vf dP
P= , = 0.
t f + td dt f

Capacidade: 500 kg de sólido/dia (10 horas de operação).

Propriedades do sistema: densidade do fluido 1 g/cm3, viscosidade do fluido 0,8 cP,
densidade do sólido 1,9 g/cm3.
Concentração da suspensão: 70 g de sólido/L de água.
Propriedades da torta: porosidade 0,6 e permeabilidade 8×10-9 cm2.
Tempo de desmantelamento do filtro, retirada da torta e montagem: 30 min.
A torta é pouco compressível e a resistência do meio filtrante pode ser desprezada em relação
à resistência da torta.

140
Resposta:

O tempo ótimo de filtração é igual ao tempo de desmantelamento limpeza e montagem do

filtro. Portanto, o tempo ótimo do ciclo é 60 min.
Volume de filtrado por ciclo: 714 L.
Área de filtração: 1,14 m2.
Espessura da torta no final do ciclo: 6 cm.

11. No filtro concentrador contínuo e pressurizado Hirondelle a superfície filtrante fica

submersa na suspensão e a torta é dela retirada por meio de raspador rotativo de 2 facas.
Sabendo-se que o filtro opera a pressão constante, estabelecer a relação entre a concentração
da suspensão espessada e: a capacidade do filtro, a pressão de operação, propriedades da
alimentação, características da torta, dimensões do elemento filtrante e velocidade do
raspador.

141
Resposta:

1 µ  < α > Pρ F C 
=  + Rm  ,
2 P A∆p  2 AN 

( )
onde P é a produção de filtrado L3 / θ , N o número de rotações por unidade de tempo e C a
concentração média da suspensão no interior da carcaça do filtro.

12. Calcular o diâmetro e a altura do sedimentador Dorr-Oliver para operar com 30 m3/h de
suspensão aquosa de cal. Dados: concentração de sólido na alimentação 0,08 g/cm3 de
suspensão, concentração de sólido no lodo 0,25 g/cm3 de suspensão, densidade da cal
2,2 g/cm3 e temperatura 25ºC.

Ensaio de aproveta a 25ºC (0,08 g/cm3 de suspensão):

Distância da interface
clarificada ao fundo da 40 32,8 25,5 18,8 14,2 11,2 9,6 6,6 5,2 4,0
proveta, z (cm)

Tempo de
sedimentação, 0 5 10 15 20 25 30 35 40 45
θ (min)

142
Cálculo da capacidade (método de Kynch) e da altura do sedimentador (Damasceno e
Massarani, 1993):

Ensaio de proveta

z0

ponto crítico

z0
ca z0
z(t) cl

t >0 tmin t
tr

Capacidade de projeto

 F z
  = 0 .
 A  proj. tmin

Altura da região de compactação (balanço de massa)

4 Fca t r ρS − ρF
Hc = .
3 AρS ρ − ρF
l

Nomenclatura

A - Área da seção transversal do sedimentador ( L2 ) .

c - Concentração de sólido na suspensão (M sólido/L3 suspensão).
F - Vazão de suspensão na alimentação (L3/θ).
Hc - Altura da região de compactação ( L ) .
t - Tempo de sedimentação (θ) .
tmi - Tempo assinalado na curva de sedimentação ( θ) .
tr - Tempo de resistência da partícula sólida na região de compactação (θ) .

143
 z - Distância entre a interface inferior da região clarificada e a base da proveta ( L ) .
ρ - Densidade da fase ( M / L3 ) .
Índices:

a - Alimentação do sedimentador;
F - Líquido;
l - Lodo que deixa o sedimentador;
S - Sólido.

Resposta:

Diâmetro do sedimentador: 6,1 m.

Altura do sedimentador: 1,3 m, incluindo 0,6 m para a separação do “extravasante” e 0,45 m
levando em conta a inclinação do fundo.

13. A indústria de papel Bananal Paulista estuda a possibilidade da utilização de um

sedimentador Dorr-Oliver com 23 m de diâmetro e 3 m de altura para o tratamento de licor
negro. Calcular a capacidade do sedimentador para as seguintes condições operacionais:
6,7 g/L a concentração de sólidos na alimentação e 19 g/L a concentração de sólidos no lodo.
Densidade do sólido, 2,8 g/cm3. Temperatura: 25ºC.

Ensaio de proveta a 25ºC (6,7 g/L de suspensão):

z(cm) 30 26,5 23,2 16,6 13,5 12,4 11,2 10,4 10,2

θ(min) 0 2,5 5 10 15 20 30 50 70.

Resposta:

A altura do sedimentador não é o fator limitante; a capacidade recomendada é da ordem de

160 m3/h de alimentação.

14. Deseja-se estimar a capacidade de um sedimentador lamelado no tratamento de uma

suspensão floculenta de hidróxido de alumínio: concentração inicial 4,5 g/L e concentração
final 22 g/L. O sedimentador funciona em contacorrente com 30 lamelas ativas, 1,80x2,00 m,
espassamento 6 cm e inclinação de 60º com a horizontal.

Ensaio de proveta (4,5 g/L de suspensão):

z(cm) 35 32,2 27,4 20,6 16,2 11,2 8,5 6,4 5,0 4,2
θ(min) 0 3 7 13 18 25 30 35 40 45.

144
O procedimento para o cálculo da área de sedimentação do sistema lamelado é o mesmo do
sedimentador Dorr-Oliver, devendo-se considerar a soma das áreas projetadas das lamelas
ativas na horizontal.

Resposta:

Área de sedimentação: 54 m2.

Capacidade do sedimentador lamelado: 35 m3/h de alimentação.

Bibliografia

Damasceno, J.J.R., “Uma Contribuição ao Estudo do Espessamento Contínuo”, Tese de

D.Sc., Programa de Eng. Química, COPPE/UFRJ, 176 p. (1992).

Damasceno, J.J.R. e Massarani, G., “Projeto e Análise do Desempenho de Sedimentadores

Contínuos”, Ciência & Engenharia/UFU, ano 2, nº 2, 61-76 (1993).

França, S.C.A., “Operação de Espessadores Não-Convencionais”, Tese de M.Sc., Programa

de Engª Química, COPPE/UFRJ, 53 p. (1996).

Massarani, G., “Filtração”, Revista Brasileira de Tecnologia, número especial, 69 p. (1985).

Massarani, G., Silva Telles, A. e Damasceno, J.J.R., “Evaluation of the Compression-

Permeability Behavior of Sediments Subjected to Small Deformations”, Anais do 6th World
Filtration Congress, Nagoya, Japão, maio, 91-95 (1993).

Shirato, M., Sambuichi, M. e Murase, T., “Hydraulic Pressure Distribution in Filter Cakes
Under Constant Pressure Filtration”, Memoir of the Faculty of Engineering, Nagoya
University, vol. 16, nº 1 e 2, 68-79 (1964).

Silva Telles, A., Cohen, B.M., Massarani, G. e Leite, M.S., “Filtração a Pressão Constante”,
Anais do I ENEMP, Rio de Janeiro, 7 p. (1973).

Silva Telles, A. e Massarani, G., “Compactação de Meios Deformáveis: Uma Descrição

Segundo a Mecânica do Contínuo”, Anais do X ABCM, Rio de Janeiro, dezembro, vol. 2, 13-
16 (1989).

Tiller, F.M., “Theory of Filtration”, IEC, vol. 53, nº 7, 529-537 (1961).

Tiller, F.M. e Chen, W., “Limiting Operating Conditions for Continuous Thickners”, Chem.
Engng. Sci., vol. 43, nº 7, 1695-1704 (1988).

145
Índice Onomástico

Allen, T. 11
Almeida, O.P. 18, 29 Ergun, S. 71
Almendra, E.R. 21, 23 Fernandes, R.C. 62, 65
Ávila, d’, J.S. 102 Ferrandon, J. 65
Aziz, K. 20, 24 Ferreira, M.C. 107
Barnea, E. 21 França, S.C.A. 128, 130
Barrientos, A. 13, 45 Francis, A.W. 18
Bear, J. 65, 75 Freire, J.T. 76, 107
Berker, R. 7, 26, 69 Gauvin, W.H. 28
Bird, R.B. 11, 41, 43 Gidaspow, D. 98, 105
Bloor, M.I.G. 45 Govier, G.W. 20, 24
Bolger, J.C. 32 Grace, J.R. 26
Brauer, H. 38 Green, D.W. 44, 74
Caswell, B. 28 Grodzinski, P. 34
Chen, W. 129 Gubulin, J.C. 101
Christiansen, A. 13 Hackenberg, C.M. 26
Christiansen, E.B. 12 Haider, A. 13
Churchill, S.W. 13 Inghan, D.B. 45
Clark, N.N. 28 Kelly, P.D. 62
Clift, R. 26 Krauss, M.N. 106
Coelho, G.L.V. 109 Lapple, C.E. 14
Coelho, R.M.L. 14 Larerack, S.D. 45
Cohen, B.M. 124 Laruccia, M.B. 29
Concha, F. 13, 21, 45 Leite, M.S. 124
Coulson, J.M. 72, 85 Leith, D. 45
Crochet, M.J. 62 Levenspiel, O. 13
Damasceno, J.J.R. 93, 120, 128, 143 Licht, W. 45
Darcy, H. 64 Lightfoot, E.N. 11, 41, 43
Löffler, F. 45

146
Marchildon, E.K. 28 Silva Telles, A. 28, 62, 64, 65, 88, 108,
Massarani, G. 14, 21, 22, 28, 46, 64, 65, 118, 124
66, 70, 71, 73, 88, 97, 98, Smith, G.F. 65
100, 101, 102, 108, 109, Sobreiro, L.E.L. 110
118, 124, 127, 128, 143 Stewart, W.E. 11, 41, 43
Michaels, A.S. 32 Stokes, C.R. 6
Mizrahi, J. 21 Svarovski, L. 45, 55
Mothes, H. 45 Tiller, F.M. 123, 129
Müller, I. 62 Tobinaga, S. 76, 77
Munroe, B. 18 Truesdell, C. 62
Murase, T. 124 Turton, R. 28
Naghdi, P.M. 62 Wasp, E.J. 106
Norton Chem. Proces. Products 78 Weber, M.E. 26
Pereira, C.M.S. 46 Zaki, W.N. 20, 25
Perry, R.H. 44, 74
Pettyjohn, E.S. 12
Politis, T. 21
Polubarinova-Kochina, P.Ya. 66
Renganathan, K. 28
Restini, C.V. 102, 107
Richardson, J.F. 20, 25, 72, 85
Samarco 108
Sambuichi, M. 124
Santana, C.C. 22, 70, 102, 105, 107, 109
Scheidegger, A.E. 65, 77
Schwarz, W.H. 28
Shepherd, C.B. 14
Shirato, M. 124
Silva, E.M.V. 107

147