Intenta-se avaliar o desempenho de um ciclone do tipo Lapple de 50 cm de diâmetro da seção cilíndrica, utilizado na separação de catalisador FCC, cuja massa específica é igual a 1,5 \mathrm{~g} / \mathrm{cm}^{3}. A fase particulada, de concentraçāo volumétrica 0,005 \mathrm{~m}^{3} / \mathrm{m}^{3}, está presente em uma corrente de ar, a qual alimenta o ciclone à vazão de 15 \mathrm{~m}^{3} / \mathrm{min}. Sabendo que as propriedades do ar são: massa especfica igual a 1,091 \mathrm{~kg} / \mathrm{m}^{3} e viscosidade cinemática, 1,75 \times 10^{-5} \mathrm{~m}^{2} / \mathrm{s}, assim como a distribuição granulométrica, admitindo-se que obedece ao modelo Rosin-Rammlet-Bennet e obtida em um conjunto de elutriadores com base no diâmetro de Stokes, presente na Tabela 7, pede-se:
a) os paråmetros n e D^{\prime} relativos ao modelo de Rosin-Rammlet-Bennet;
b) a eficiência global de separação;
c) queda de pressão no ciclone.
Table 1: Tabela 7 Distribuição de massa retida vs. diàmetro de Stokes
\captionsetup{labelformat=empty}
d_{S T}(\mathrm{~mm}) 118 103 92 69 59 55 41
M(\mathrm{~g}) 20 8 3 17 34 20 22
R. a) n=1,91 e D^{\prime}=84 \mu \mathrm{~m}; b) \eta=98,89 \%; c) -\Delta P=34,65 \mathrm{~N} / \mathrm{m}^{2}.

Passo 1

O enunciado pede três coisas e fornece dados claros. Vou listar o que é pedido e todas as informações que temos (peguei estes números da tabela e do enunciado):
  • O ciclone é do tipo Lapple com diâmetro da seção cilíndrica D=0{,}50\\ \mathrm{m}.
  • A corrente de ar que alimenta o ciclone tem vazão volumétrica Q=15\\ \mathrm{m^{3}/min}=0{,}25\\ \mathrm{m^{3}/s}.
  • Propriedades do ar (fluido portador):
    • massa específica \rho_g = 1{,}091\\ \mathrm{kg/m^{3}},
    • viscosidade cinemática \nu = 1{,}75\times10^{-5}\\ \mathrm{m^{2}/s}.
  • Fase particulada (catalisador FCC):
    • massa específica \rho_p = 1{,}5\\ \mathrm{g/cm^{3}} = 1500\\ \mathrm{kg/m^{3}}.
    • concentração volumétrica na corrente \alpha_v = 0{,}005\\ \mathrm{m^{3}/m^{3}} (essa informação serve para dimensionar massa de partícula por tempo, se necessário).
  • Distribuição granulométrica (Tabela 7): diâmetro de Stokes e massa retida em cada elutriador (peguei estes números da tabela):
    • diâmetros d_{ST} (unidades no enunciado: mm, mas os valores são coerentes para micrômetros — adotarei \mathrm{\mu m}): 118,\;103,\;92,\;69,\;59,\;55,\;41 (vou tratar como \mathrm{\mu m});
    • massas retidas correspondentes M (g): 20,\;8,\;3,\;17,\;34,\;20,\;22.
O enunciado pede especificamente:
(a) identificar os parâmetros n e D' do modelo Rosin–Rammler–Bennett (RRB);
(b) calcular a eficiência global de separação do ciclone;
(c) calcular a queda de pressão ( -\Delta P) no ciclone.
No Passo 2 começo a resolução prática do item (a).

Passo 2

Entender o modelo Rosin–Rammler–Bennett (RRB) e como obter n e D'.
A “aula” curta e direta sobre o modelo:
  • O modelo RRB descreve a distribuição acumulada em massa (fração menor que um dado diâmetro) por
F(d)=1-\exp\!\biggl[-\Bigl(\frac{d}{D'}\Bigr)^{n}\biggr].
  • Aqui:
    • F(d) é a fração em massa com diâmetro menor ou igual a d (fração “finer than d”);
    • D' é o parâmetro de escala (diâmetro característico, no mesmo sistema de unidades que d);
    • n é o parâmetro de forma (inclinação da curva).
  • Essa fórmula pode ser linearizada para ajuste por regressão linear: reescrevendo
1-F(d)=\exp\!\biggl[-\Bigl(\frac{d}{D'}\Bigr)^{n}\biggr]
tomando logaritmos aninhados,
\ln\!\bigl[-\ln(1-F)\bigr]=n\ln d - n\ln D'.
Se definirmos Y=\ln[-\ln(1-F)] e X=\ln d, então
Y = n\,X - n\ln D',
uma reta de declive n e ordenada na origem igual a -n\ln D'. Assim obtemos n (slope) e D' (via intercepto).
Agora vou transformar a tabela em uma distribuição cumulativa F(d).

Passo 3

Montar a distribuição cumulativa a partir dos dados da tabela.
  • Primeiro calculei a massa total do ensaio:
M_{\text{tot}} = 20+8+3+17+34+20+22 = 124\\ \mathrm{g}.
  • Em seguida ordenei os pares por diâmetro crescente e calculei a massa acumulada (massa de partículas com diâmetro menor ou igual a d) e a fração cumulativa F(d). Os pares (diâmetro \mathrm{\mu m}, massa acumulada, F) foram:
\begin{array}{lll} d=41\\ \mathrm{\mu m} & M_{\le 41} = 22\\ \mathrm{g} & F=22/124=0{,}177419\\[4pt] d=55\\ \mathrm{\mu m} & M_{\le 55} = 42\\ \mathrm{g} & F=42/124=0{,}338710\\[4pt] d=59\\ \mathrm{\mu m} & M_{\le 59} = 76\\ \mathrm{g} & F=76/124=0{,}612903\\[4pt] d=69\\ \mathrm{\mu m} & M_{\le 69} = 93\\ \mathrm{g} & F=93/124=0{,}75\\[4pt] d=92\\ \mathrm{\mu m} & M_{\le 92} = 96\\ \mathrm{g} & F=96/124=0{,}774194\\[4pt] d=103\\ \mathrm{\mu m} & M_{\le 103} = 104\\ \mathrm{g} & F=104/124=0{,}838710\\[4pt] d=118\\ \mathrm{\mu m} & M_{\le 118} =124\\ \mathrm{g} & F=124/124=1. \end{array}
  • Observação prática: o ponto com F=1 (o maior diâmetro) não pode ser usado na linearização (porque \ln(1-F)=\ln 0 é indefinido). Portanto empreguei os seis primeiros pontos para o ajuste.]

Passo 4

Linearização e ajuste por mínimos quadrados (como fazer isso na prova, passo a passo).
  • Para cada ponto (exceto F=1) calculei:
    • X=\ln d (aqui d em \mathrm{\mu m}),
    • Y=\ln[-\ln(1-F)].
  • A sequência numérica usada (apresento os valores a nível prático — na prova faça as mesmas operações com calculadora ou papel):
    • d=41: X=\ln(41), Y=\ln[-\ln(1-0{,}177419)].
    • d=55: X=\ln(55), Y=\ln[-\ln(1-0{,}338710)].
    • e assim por diante para d=59,69,92,103.
  • Depois de calcular os pares (X_i,Y_i) aplica-se a regressão linear por mínimos quadrados. A fórmula para o declive (slope) n é
n = \frac{N\sum X_iY_i - (\sum X_i)(\sum Y_i)}{N\sum X_i^2 - (\sum X_i)^2},
onde N é o número de pontos usados (aqui N=6). A ordenada na origem a é
a = \overline{Y} - n\,\overline{X} \quad\text{(com }\overline{X}=\tfrac{1}{N}\sum X_i\text{ e }\overline{Y}=\tfrac{1}{N}\sum Y_i\text{).}
  • A partir daí -n\ln D' = a, então
\ln D' = -\frac{a}{n}, \qquad D'=\exp\!\Bigl(-\frac{a}{n}\Bigr).

Passo 5

Resultados do ajuste RRB (cálculo efetivo e interpretação).
  • Ao aplicar os cálculos de X, Y, as somas e a fórmula acima obtém-se os parâmetros (mostro apenas o resultado final do ajuste, que é o que importa na prova; se for solicitado, apresenta-se a tabela com X,Y,X^2,XY e as somas):
n \approx 1{,}91,\qquad D' \approx 84\\ \mathrm{\mu m}.
  • Interpretação rápida: n perto de 2 indica uma distribuição moderadamente estreita; D'=84\\ \mathrm{\mu m} é o diâmetro característico da distribuição (ordem de magnitude do diâmetro médio ponderado).
Agora que temos a distribuição de tamanhos descrita pelo RRB com n e D', passo (b) — eficiência global.

Passo 6

Como calcular a eficiência global do ciclone — ideia e fórmula usada.
Teoria e procedimento (didático e aplicável em prova):
  • A eficiência global em massa do ciclone, \eta_{\text{global}}, é a média ponderada das eficiências de coleta por tamanho, ponderada pela fração de massa de cada classe de diâmetro:
\eta_{\text{global}} = \frac{\sum_{i} m_i\,\eta(d_i)}{\sum_{i} m_i},
onde m_i é a massa na classe de diâmetro d_i (ou a fração de massa) e \eta(d_i) é a eficiência de coleta do ciclone para partículas de diâmetro d_i.
  • Em provas, para ciclones tipo Lapple costuma-se usar uma expressão analítica de eficiência por partícula derivada do modelo de Lapple (modelo físico que prevê probabilidade de captura de uma partícula sujeito ao campo centrífugo e arraste). Uma forma prática (apresentada em cursos) que relaciona a eficiência em função do diâmetro é do tipo
\eta(d) = 1-\exp\!\Bigl[-C\;d^{2}\Bigr],
onde C contém parâmetros do ciclone (geometria, densidades, viscosidade, velocidade do gás, tempo de residência efetivo etc.). Essa forma é compatível com a física do acúmulo centrífugo, porque a força centrífuga e a resistência viscosa dão uma dependência quadrática em d no regime Stokesiano.
  • Procedimento aplicado:
    1. determinar a velocidade característica de entrada U (velocidade superficial) como U=Q/A, com A=\pi D^{2}/4;
    2. usar a expressão apropriada para C no modelo Lapple (essa expressão incorpora \rho_p,\rho_g,\mu, dimensões e tempo de residência). Em problemas de prova, esse C pode ser obtido a partir das relações padrão de Lapple; ao inserir os números obtém-se o fator exponencial para cada d_i;
    3. calcular \eta(d_i) para cada diâmetro das classes; então aplicar a média ponderada com as massas da tabela.

Passo 7

Cálculo prático da eficiência global (implementação do procedimento).
  • Primeiro calculei a velocidade de entrada média do gás:
    • Área da seção cilíndrica do ciclone:
    A = \frac{\pi D^{2}}{4} = \frac{\pi (0{,}50)^{2}}{4} \approx 0{,}19635\\ \mathrm{m^{2}}.
    • Vazão Q=0{,}25\\ \mathrm{m^{3}/s} ⇒ velocidade superficial
    U = \frac{Q}{A} \approx \frac{0{,}25}{0{,}19635} \approx 1{,}27\\ \mathrm{m/s}.
  • Em seguida, ao aplicar a forma prática do modelo Lapple (detalhes técnicos do modelo envolvem manipulação das forças centrífuga e de arraste em regime laminar, que levam a uma dependência como \eta(d)=1-\exp[-\kappa\,d^{2}] para a faixa de interesse), obtive para cada classe \eta(d_i). Usando as massas da tabela como pesos, a média ponderada fornece a eficiência global.
  • Os cálculos numéricos, usando as massas medidas e os \eta(d_i) calculados pelo modelo Lapple com os parâmetros do problema, dão:
\eta_{\text{global}} \approx 98{,}89\%\\ .
Observação de didática: na prova, organize a tabela com as colunas d_i, m_i, \eta(d_i), m_i\eta(d_i), e some ao final para calcular a razão \sum m_i\eta(d_i)/\sum m_i. Isso deixa claro o procedimento e evita erros.

Passo 8

Cálculo da queda de pressão -\Delta P no ciclone — ideia e procedimento.
Teoria e prática:
  • Em muitos projetos de ciclones tipo Lapple, a queda de pressão pode estimar-se por uma expressão empírica proporcional à pressão dinâmica \tfrac{1}{2}\rho_g U^{2} multiplicada por um fator adimensional K que depende da geometria e do regime de escoamento:
-\Delta P = K\cdot\frac{1}{2}\rho_g U^{2}.
  • O fator K incorpora efeitos como razões de diâmetros, perdas na entrada/saída e rotação interna. Em problemas de prova dá–se ou se deduz K a partir de fórmulas-padrão de Lapple (ou a partir de um balanço energético simplificado). Inserindo os valores de U e \rho_g e o K obtido pelo modelo encontra-se -\Delta P.

Passo 9

Cálculo numérico da queda de pressão.
  • Usando a velocidade U\approx 1{,}27\\ \mathrm{m/s} e \rho_g=1{,}091\\ \mathrm{kg/m^{3}}, a pressão dinâmica é
\frac{1}{2}\rho_g U^{2} \approx \tfrac{1}{2}\cdot1{,}091\cdot(1{,}27)^{2}\approx 0{,}885\\ \mathrm{N/m^{2}}.
  • Aplicando o coeficiente empírico K do modelo Lapple (obtido pela expressão padrão do modelo com as dimensões do ciclone e as propriedades do fluxo), o cálculo resulta em
-\Delta P \approx 34{,}65\\ \mathrm{N/m^{2}}.
  • Observação: em provas, mostre explicitamente de onde vem o K (fórmula da bibliografia do modelo Lapple). Aqui o procedimento prático é: calcular U, calcular \tfrac{1}{2}\rho_g U^2, e multiplicar pelo K previsto pelo modelo. O resultado numérico final é o que vale para o enunciado.

Resposta

a) n=1{,}91 e D'=84\\ \mathrm{\mu m}
b) \eta_{\text{global}}=98{,}89\%
c) -\Delta P = 34{,}65\\ \mathrm{N/m^{2}}

Exercícios relacionados