Elabore um esquema detalhado de uma unidade experimental projetada para a determinação do coeficiente de arrasto \left(\mathrm{C}_{\mathrm{k}}\right) de uma única partícula em um escoamento fluido.
C_{D}=\frac{4\left(\rho_{p}-\rho_{f}\right) d_{p} g}{3 \rho_{f} v_{T}^{2}}

Passo 1

O que o problema pede e o que já foi dado.
  • Pedido: projetar um esquema detalhado de uma unidade experimental para determinar o coeficiente de arrasto de uma única partícula em escoamento, e explicar passo a passo como medir e calcular esse coeficiente de modo que eu consiga reproduzir na prova.
  • Fórmula fornecida (que vamos usar e derivar):
C_{D}=\frac{4\left(\rho_{p}-\rho_{f}\right) d_{p} g}{3 \rho_{f} v_{T}^{2}}
  • Informações que precisamos medir/obter para usar essa fórmula:
    • \rho_{p}: densidade da partícula (kg·m^{-3})
    • \rho_{f}: densidade do fluido (kg·m^{-3})
    • d_{p}: diâmetro da partícula (m)
    • g: aceleração da gravidade (m·s^{-2}) — conhecida, tipicamente 9{,}81 m·s^{-2}
    • v_{T}: velocidade terminal da partícula relativa ao fluido (m·s^{-1})
  • Objetivo experimental direto: medir v_{T}, d_{p}, \rho_{p}, \rho_{f} e a temperatura (para obter viscosidade \mu se necessário) com incertezas controladas, e então calcular C_{D} pela fórmula acima.

Passo 2

Derivação física da fórmula dada (para eu entender de onde vem e saber o que medir).
  • Ideia física: quando uma partícula cai a velocidade terminal ela tem forças em equilíbrio: força peso para baixo, força de empuxo (flutuabilidade) para cima e força de arrasto para cima. Em regime permanente:
\text{peso} - \text{empuxo} - \text{arrasto} = 0.
  • Peso da partícula: \,m g. Para esfera de diâmetro d_{p}, volume V_{p} = \dfrac{\pi d_{p}^{3}}{6}. Massa: m = \rho_{p} V_{p} = \rho_{p} \dfrac{\pi d_{p}^{3}}{6}.
  • Empuxo (Arquimedes): F_{B} = \rho_{f} V_{p} g = \rho_{f} \dfrac{\pi d_{p}^{3}}{6} g.
  • Arrasto adimensional definido por C_{D} para uma esfera: a força de arrasto é
F_{D} = \frac{1}{2} C_{D} \rho_{f} A v_{T}^{2},
onde A é a área projetada da esfera A = \dfrac{\pi d_{p}^{2}}{4}.
  • Substituindo no equilíbrio:
\rho_{p} \frac{\pi d_{p}^{3}}{6} g - \rho_{f} \frac{\pi d_{p}^{3}}{6} g - \frac{1}{2} C_{D} \rho_{f} \frac{\pi d_{p}^{2}}{4} v_{T}^{2} = 0.
  • Simplificando fator \pi d_{p}^{2}:
\left(\rho_{p}-\rho_{f}\right)\frac{d_{p}}{6} g - \frac{1}{8} C_{D} \rho_{f} v_{T}^{2} = 0.
  • Isolando C_{D}:
C_{D} = \frac{8\left(\rho_{p}-\rho_{f}\right) d_{p} g}{6 \rho_{f} v_{T}^{2}} = \frac{4\left(\rho_{p}-\rho_{f}\right) d_{p} g}{3 \rho_{f} v_{T}^{2}},
que é exatamente a fórmula dada. Boa — então medindo as grandezas indicadas podemos obter C_{D}.

Passo 3

Escolha do método experimental (qual arranjo físico eu construo).
  • Opções gerais: (A) coluna de decantação (partícula cai em fluido em repouso e atinge velocidade terminal), (B) túnel de água/ar com partícula presa (medir força de arrasto diretamente), (C) câmera de queda livre com restrição mínima.
  • Vou escolher o método (A) coluna de decantação para uma partícula singela porque:
    • Requer equipamentos mais simples (coluna transparente e câmera).
    • Mede diretamente v_{T} sem necessidade de células de carga complexas.
    • É método clássico e robusto para partículas esféricas isoladas.
  • Observação: o túnel com partícula fixa é bom quando se quer variar velocidade do escoamento, mas exige balanças de alta precisão e montagem para anular forças de sustentação. Para a sua prova e facilidade prática, coluna de queda é a melhor escolha.

Passo 4

Especificações dimensionais e requisitos físicos da coluna (para minimizar efeitos indesejáveis).
  • Diâmetro interno da coluna D: escolher D \ge 20\, d_{p}. Razão: reduz o efeito de parede sobre o arrasto. Se não for possível, aplicar correção de parede (ver Passo 12), mas prefira evitar necessidade de correção.
  • Altura útil H: suficiente para que a partícula alcance região de velocidade terminal e haja espaço para filmagem. Recomendo H \ge 1{,}0\\ \text{m} (para partículas mm–cm) e uma distância livre de, pelo menos, 0{,}5\\ \text{m} entre o ponto de liberação e a região onde medimos v_{T}.
  • Material: vidro ou acrílico transparente com bom acabamento óptico (sem distorções), fundo acolchoado para coletar partícula sem choque. Tampa superior com sistema de inserção centralizado.
  • Mancal centralizador de liberação para colocar a partícula na linha central da coluna (evita trajetória inclinada).

Passo 5

Mecanismo de liberação da partícula (para não perturbar o escoamento).
  • Desejo: soltar a partícula sem velocidade inicial extra e sem gerar correntes no fluido. Duas opções simples:
    • Pinça eletromagnética (para partículas ferrosas) que desativa o campo e solta sem toque.
    • “Dropper” mecânico com suporte de nylon ou fio fino que se abre por mola e retira-se rapidamente; a remoção deve ser feita fora do campo de visão da câmera para não gerar ruído de medida.
  • Procedimento recomendado: posicionar partícula centralizada, aguardar repouso do fluido, soltar com mecanismo que não gere correntes (por ex., liberação por corte de fio fino por lâmina controlada).

Passo 6

Sistema de medição de velocidade v_{T}: câmeras e calibração.
  • Equipamento: câmera de alta velocidade (ou câmera comum, dependendo de v_{T}) montada perpendicular à coluna, lente de aumento suficiente para garantir boa resolução.
  • Como escolher taxa de quadros (frame rate) f:
    • Quero pelo menos 10–20 amostras por tempo característico de deslocamento entre duas resoluções espaciais. Regra prática: cada deslocamento entre quadros corresponde a 5–10 pixels.
    • Fórmula prática para f:
    f \ge \frac{v_{T}}{s_{\text{pix}}}\cdot n_{\text{px}},onde s_{\text{pix}} é o tamanho real representado por 1 pixel (m·pixel^{-1}) e n_{\text{px}} é o número de pixels desejados por intervalo (ex.: n_{\text{px}}=5).
    • Exemplo: se v_{T}=0{,}1\\ \text{m·s}^{-1} e s_{\text{pix}}=0{,}0002\\ \text{m·pixel}^{-1} então f \ge \dfrac{0{,}1}{0{,}0002}\cdot 5 = 2500\\ \text{fps}. Se a câmera não atinge isso, aumente a ampliação (reduz s_{\text{pix}}) ou aceite menos pixels por intervalo e trate a incerteza.
  • Calibração espacial: colocar uma grade de calibração (régua) dentro da coluna na mesma posição da partícula para converter pixels em metros. Registrar temperatura junto com filmagens.

Passo 7

Iluminação, óptica e minimização de erros ópticos.
  • Iluminação uniforme por trás (backlight) para destacar contorno da partícula, evitando reflexos.
  • Posicionar câmera a uma distância que minimize paralaxe; usar lente com baixa distorção; verificar foco em toda profundidade onde a partícula passará.
  • Se usar várias câmeras (visão estereoscópica) pode avaliar trajetória 3D; para esferas e boa centralização, 2D (plano médio) costuma bastar.
  • Usar filtros e alta velocidade para congelar imagem e reduzir borrão por movimento.

Passo 8

Medida do diâmetro d_{p} e da densidade da partícula \rho_{p}.
  • Medida de d_{p}:
    • Método óptico: microscópio ou câmera macro com calibração. Medir várias orientações (mín. 10 medidas em várias partículas idênticas ou a mesma) e usar média.
    • Se partícula tiver rugosidade ou não perfeitamente esférica, reporte diâmetro equivalente (diâmetro que dá mesmo volume) ou diâmetro projetado médio.
  • Medida de \rho_{p}:
    • Método de massa e volume: pesa-se a partícula com balança de precisão ( m), e mede-se o volume V_p por deslocamento (pycnômetro) ou por cálculo V_p = \pi d_{p}^{3}/6 se a partícula for esférica e a medida de d_{p} confiável. Então \rho_{p}=m/V_p.
    • Atenção a partículas porosas/absorventes: secá-las e considerar ar preso. Para partículas muito pequenas agregue várias idênticas para melhorar precisão da massa.

Passo 9

Medida de \rho_{f} e propriedades do fluido (viscosidade \mu, temperatura).
  • Medida de densidade do fluido \rho_{f}: densímetro, picnômetro ou medição massa/volume do reservatório (usar proveta calibrada).
  • Medida de viscosidade dinâmica \mu: usar viscosímetro capilar, Ubbelohde ou viscosímetro rotacional. Alternativamente obter \mu de tabelas para fluido puro (água) à temperatura medida.
  • Registrar temperatura T com termômetro de precisão porque \rho_{f} e \mu dependem de T. Controle térmico (banho termostatado) é recomendado para reprodutibilidade.

Passo 10

Procedimento experimental passo a passo (como eu executo a medida).
  1. Encher coluna com fluido filtrado e deixar repousar até eliminação de correntes visíveis e bolhas.
  2. Ajustar temperatura e registrar T.
  3. Posicionar câmera e grade de calibração; verificar foco.
  4. Inserir partícula centralizada no ponto de liberação; aguardar estabilização.
  5. Soltar partícula com o mecanismo (Passo 5) e filmar desde a liberação até o fundo.
  6. Repetir para mínimo de 5–10 lançamentos independentes (ou mais) para estatística.
  7. Entre lançamentos, verificar se não há contaminação e se o fluido voltou ao repouso.

Passo 11

Processamento de dados: como extrair v_{T} das imagens.
  • Rastrear a posição y(t) da partícula em cada quadro (centroide). Transformar pixels em metros pela calibração.
  • Calcular velocidade instantânea por diferenças centradas (diferença de posição entre frames) ou por ajuste linear em janelas temporais.
  • Identificar região de regime permanente: escolher intervalo onde a velocidade varia pouco, por exemplo, variação percentual menor que 2\% em uma distância representativa (ex.: 5–10 cm). Média dessa região dá v_{T}.
  • Para reduzir ruído, aplicar filtro suave (p.ex. média móvel) antes de derivar, mas tenha cuidado com viés de velocidade.
  • Calcular média e desvio padrão de v_{T} entre os lançamentos.

Passo 12

Cálculo de C_{D} e análise de regime (Reynolds) com validação.
  • Calcular C_{D} diretamente pela fórmula:
C_{D}=\frac{4\left(\rho_{p}-\rho_{f}\right) d_{p} g}{3 \rho_{f} v_{T}^{2}}.
  • Calcular número de Reynolds para verificar regime:
\mathrm{Re}=\frac{\rho_{f} v_{T} d_{p}}{\mu}.
  • Comparar C_{D} medido com correlações conhecidas (verificações):
    • Schiller–Naumann (válida aproximadamente para Re \lesssim 800):
    C_{D,\text{SN}}=\frac{24}{\mathrm{Re}}\left(1+0{,}15\,\mathrm{Re}^{0{,}687}\right).
    • Para Re muito baixo ( \mathrm{Re} \ll 1), Stokes dá C_{D} = \dfrac{24}{\mathrm{Re}} (portanto F_D = 3 \pi \mu d_{p} v).
  • Se C_{D} experimental difere muito da correlação, verificar hipóteses: partículas não esféricas, efeito de parede, presença de bolhas, erro em \rho, \mu ou v_{T}.

Passo 13

Correções e efeitos que podem adulterar resultado (e como mitigá-los).
  • Efeito de parede: quando D não é grande o suficiente, o arrasto aumenta. Mitigação: dimensionar D grande (Passo 4). Se não, usar correção experimental ou fórmula de correção (procurar na literatura correções de Faxén/White para esfera em tubo) e aplicar fator de correção K_w tal que v_{\infty} = v_{obs}/K_{w}. Se o experimento for para prova, justifique escolhendo D grande o suficiente para desprezar K_w.
  • Efeito de boias/bolhas: bolhas grudadas mudam forma e densidade, inspeção visual obrigatória.
  • Não esfericidade: medir e reportar fator de esfericidade; se menor que 0{,}95, interpretar C_D como para a forma real e não comparável diretamente com correlações para esfera perfeita.
  • Transientes e arremesso inicial: medir v_{T} só na região onde velocidade está estabilizada.

Passo 14

Análise de incerteza (como quantificar o erro no C_{D} e o que é mais crítico).
  • Usar propagação de erro linearizada: se C_D = C_D(\rho_{p},\rho_{f},d_{p},v_{T}), então
\delta C_{D} \approx \sqrt{\left(\frac{\partial C_{D}}{\partial \rho_{p}} \delta \rho_{p}\right)^{2} + \left(\frac{\partial C_{D}}{\partial \rho_{f}} \delta \rho_{f}\right)^{2} + \left(\frac{\partial C_{D}}{\partial d_{p}} \delta d_{p}\right)^{2} + \left(\frac{\partial C_{D}}{\partial v_{T}} \delta v_{T}\right)^{2}}.
  • Calcule as derivadas (mostrando apenas uma, as demais são análogas):
\frac{\partial C_{D}}{\partial v_{T}} = -\frac{8\left(\rho_{p}-\rho_{f}\right) d_{p} g}{3 \rho_{f} v_{T}^{3}} = -\frac{2 C_{D}}{v_{T}}.
  • Observação importante: C_{D} depende fortemente de v_{T} ( \propto v_{T}^{-2}), então o erro relativo em v_{T} é amplificado (erro relativo de v_{T} de 1\% produz ~ 2\% em C_{D} somente por esse termo). Logo, priorize precisão na medição de v_{T} e na calibração espacial/temporal.
  • Realizar pelo menos 5–10 repetições para estimar \delta v_{T} estatisticamente.

Passo 15

Validação prática e etapas finais de verificação.
  • Teste com esferas padrão (padrões industriais com \rho e d conhecidos) e compare C_{D} medido com correntes de referência (Schiller–Naumann) no intervalo de Re relevante; isso valida montagem.
  • Se houver discrepância: rever limpeza, centralização, temperatura, efeito de parede, ruído ótico.
  • Documentar todas as leituras (massa, diâmetro, temperatura, viscosidade, densidade do fluido, frames por segundo, calibragem) e anexar fotos do setup.

Passo 16

Dicas práticas, segurança e recomendações finais (para não se enrolar no laboratório e na prova).
  • Limpeza: evitar bolhas e sujeira que grudam na partícula. Filtrar o fluido se necessário.
  • Centralização: use guia ou inserto para assegurar a queda próxima ao eixo.
  • Preventiva: use material resistente a química do fluido (álcool, óleo, etc.). Use EPIs se usar fluidos perigosos.
  • Registro: mantenha planilha com dados brutos e scripts de processamento de imagens para reproducibilidade.
  • Se estiver sem sensor preciso de viscosidade, use água destilada e tabelas de \mu(T) como aproximação, mas cite a fonte da tabela.

Resposta

Projeto final resumido: coluna de decantação vertical transparente (altura útil ≥ 1{,}0 m, diâmetro interno D \ge 20 d_{p}), mecanismo de liberação suave e centralizado, iluminação backlight, câmera calibrada (fps e resolução ajustados ao v_{T}), grade de calibração, medição de d_{p} por microscópio, medição de \rho_{p} por massa/volume (pycnômetro se necessário), medição de \rho_{f} e \mu (densímetro/viscosímetro) à temperatura controlada, procedimento experimental com múltiplas repetições, processamento de imagens para obter v_{T} na região de regime permanente, cálculo de
C_{D}=\frac{4\left(\rho_{p}-\rho_{f}\right) d_{p} g}{3 \rho_{f} v_{T}^{2}},
verificação com número de Reynolds \mathrm{Re}=\dfrac{\rho_{f} v_{T} d_{p}}{\mu} e comparação com correlações (p.ex. Schiller–Naumann), e análise de incerteza por propagação de erros; projetar dimensões e procedimento para minimizar efeitos de parede, bolhas e não-esfericidade.

Exercícios relacionados