Elabore um esquema detalhado de uma unidade experimental projetada para a determinação da velocidade terminal em sistemas, avaliando a influência da concentração das partículas e da presença de contornos rígidos.

Passo 1

O que o enunciado pede: projetar uma unidade experimental detalhada para determinar a velocidade terminal de partículas em suspensão, de modo a avaliar (i) a influência da concentração de partículas e (ii) a influência da presença de contornos rígidos (paredes/confinamento).
Informações e variáveis relevantes que vamos controlar/medir (lista):
  • tipo de partículas: diâmetro \;d\;, densidade \;\rho_p\;, polidispersidade;
  • fluido: densidade \;\rho_f\;, viscosidade dinâmica \;\mu\;, temperatura \;T\;;
  • gravidade \;g\; (constante);
  • coluna experimental: diâmetro interno \;D_c\;, altura útil \;H\;, material transparente;
  • concentração volumétrica de partículas (fração volumétrica) \;\varphi\;;
  • velocidade terminal média medida \;V(\varphi,D_c)\;;
  • outras grandezas auxiliares: número de Reynolds \;Re\;, área projetada da partícula \;A\;, volume da partícula \;V_p\;.
    O objetivo experimental final: obter V como função de \varphi e de parâmetros de confinamento (por exemplo razão \;D_c/d\;), quantificar parâmetros empíricos (ex.: expoente da lei de hindered settling) e caracterizar a correção de parede.

Passo 2

Descrição geral do equipamento (componentes principais e por que cada um é necessário):
  • coluna cilíndrica transparente (material: vidro ou acrílico) com altura útil \;H\gtrsim 1{,}0\\ \mathrm{m}\; para garantir região de regime estabelecido; razão maior que \;10 vezes a altura necessária para atingir velocidade terminal. Ter pelo menos quatro colunas/intercambiáveis com diâmetros internos diferentes para estudar confinamento: escolher \;D_c/d\approx 5,\;10,\;20,\;50\;.
  • suporte rígido e nivelado para evitar vibração.
  • sistema de aquecimento/resfriamento com circulação para controlar \;T\; dentro de \pm 0{,}2^\circ\mathrm{C} (viscosidade é sensível à temperatura).
  • câmera de alta velocidade ou câmera científica (resolução e fps escolhidos conforme V e d); exemplo prático: \;200\\ \mathrm{fps} para partículas que caem a algumas dezenas de \mathrm{mm/s}, ou \;1000\\ \mathrm{fps} para partículas menores/mais rápidas. Lente com campo de visão vertical, marcação de escala (régua) acoplada na frente. Iluminação de fundo (LED) para alto contraste.
  • sistema de injeção homogênea da suspensão na parte superior (seringa larga, placa perfurada) para evitar correntes iniciais. Tampas com telas finas que deixam passar partículas mas evitam arraste de fluxo durante a entrada.
  • agitador magnético/ultrassom para preparar suspensão homogênea antes da injeção.
  • balança analítica e provetas para preparar \varphi com precisão.
  • termômetro e viscosímetro (ou tabela viscosidade versus mistura de glicerol-água) para determinar \mu.
  • software de tracking de partículas (ex.: Trackpy, ImageJ + plugin) para extrair trajetórias z(t) e velocidades instantâneas.
  • opcional: PIV/LDV/ultrassom caso o fluido seja opaco ou para medir campo de velocidade do fluido e verificar fluidização local.

Passo 3

Escolha das partículas e do fluido — critérios práticos e exemplos:
  • partículas: escolher esferas quase monodispersas (ex.: esferas de vidro ou poliestireno). Exemplos práticos de diâmetros: \;d = 0{,}2\\ \mathrm{mm},\;0{,}5\\ \mathrm{mm},\;1{,}0\\ \mathrm{mm}\;. Densidade típica do vidro \;\rho_p\approx 2{,}5\\ \mathrm{g/cm^3}. Evitar partículas porosas ou levíssimo desgaste.
  • fluido: escolher água para baixos \mu ou misturas água–glicerol para aumentar \mu e baixar Re. Medir ou controlar \rho_f e \mu com instrumentos ou tabela.
  • motivos da escolha: variar d e \mu permite explorar regimes de Re diferentes e verificar validade de modelos (Stokes vs regime inertial).
    Como preparar a fração volumétrica \varphi: usar massas medidas \;m_p\; e \;m_f\; e densidades para calcular \varphi por
\varphi = \frac{V_p}{V_p+V_f} = \frac{m_p/\rho_p}{m_p/\rho_p + m_f/\rho_f}.
Exemplo: para preparar \varphi=0{,}05 em volume total \;V_\text{tot}, calcular m_p=\varphi V_\text{tot}\rho_p e m_f=(1-\varphi)V_\text{tot}\rho_f. Misturar e homogeneizar antes de injetar.

Passo 4

Aula rápida sobre velocidade terminal e balanceamento de forças (teoria necessária para interpretar medições):
  • ideia física: na queda de uma partícula existe peso aparente para baixo (gravidade menos empuxo) e arrasto viscoso para cima; em regime permanente essas forças se equilibram.
  • no caso de esferas muito pequenas e fluxo laminar (número de Reynolds \;Re\ll 1\;), o arrasto é dado pela Lei de Stokes. A velocidade terminal isolada de uma esfera em fluido newtoniano é:
V_0 = \frac{2(\rho_p-\rho_f)\,g\,d^2}{9\mu},
onde \;d\; é o diâmetro da esfera, \;\rho_p\; e \;\rho_f\; são densidades da partícula e do fluido, \;g\; é a gravidade e \;\mu\; a viscosidade dinâmica. Expliquei cada termo: \rho_p-\rho_f representa a diferença de densidade (força de flutuabilidade reduzida), d^2 mostra por que partículas maiores caem muito mais rápido.
  • definição do número de Reynolds para partículas:
Re = \frac{\rho_f V d}{\mu}.
Se Re\gtrsim 1 a Lei de Stokes deixa de ser válida; o arrasto tem termo quadrático e usamos o balanço geral:
\frac{1}{2}C_d\rho_f A V_t^2 = (\rho_p-\rho_f)V_p g,
onde A=\pi d^2/4 é a área projetada e V_p=\pi d^3/6 é o volume da esfera. Para obter V_t nesse caso, precisamos de uma expressão para C_d(Re) (tabelas/empíricas). Ex.: para esferas, para Re moderado pode-se usar correlações padrão (Schiller–Naumann, etc.).
  • conclusão prática: calcular V_0 teórico para o caso diluído e usar como referência; medir Re experimentalmente para saber se Stokes é aplicável.

Passo 5

Hindered settling (influência da concentração) — teoria empírica e como extrair do experimento:
  • em suspensão concentrada as partículas interagem hidrodinamicamente e geometricamente: a velocidade média de sedimentação diminui com a concentração. Uma lei empírica amplamente usada é a de Richardson–Zaki:
V(\varphi) = V_0 (1-\varphi)^n,
onde \;n\; é um expoente empírico que depende do regime de Re. Explicação: V_0 é a velocidade terminal na condição diluída (partícula isolada). O fator (1-\varphi)^n representa a redução do transporte de fluido através da rede de partículas.
  • procedimento prático para obter n: medir V(\varphi) para várias \varphi, e ajustar (em escala log) por
\ln\left(\frac{V(\varphi)}{V_0}\right) = n \ln(1-\varphi).
Ajuste linear por mínimos quadrados na forma acima fornece n. Se não puder medir V_0 experimentalmente com precisão, usar V_0 teórica de Stokes para uma estimativa inicial; preferível medir em \varphi\to 0.

Passo 6

Efeito de paredes / confinamento — como avaliar e quantificar:
  • conceito: paredes próximas alteram o padrão de escoamento em torno da partícula e aumentam o arrasto efetivo, reduzindo V. Uma forma prática de quantificar o confinamento é a razão
\Lambda = \frac{D_c}{d}.
  • estratégia experimental: manter \varphi fixo e variar \Lambda (usar colunas intercambiáveis); medir V(\Lambda). Normalizar pela velocidade em coluna muito larga (ou pelo V_0 teórico) e estudar a influência. Exemplo de curva de interesse:
\frac{V(\Lambda)}{V(\Lambda\to\infty)} \quad\text{vs}\quad \frac{d}{D_c} = \Lambda^{-1}.
  • em vez de depender de fórmulas teóricas complexas, recomendo obter um ajuste empírico do tipo
V(\Lambda) = V_{\infty}\left(1 - k_1\frac{d}{D_c} - k_2\left(\frac{d}{D_c}\right)^2\right),
onde os coeficientes k_1,k_2 são determinados pelos dados. Observação didática: para minimizar efeito de parede em medidas de referência, escolha \Lambda\gtrsim 20. Para estudar efeito intencionalmente, use \Lambda pequenos como \;5 e \;10.

Passo 7

Procedimento experimental passo a passo (prático, o que fazer no laboratório):
  1. medir e registrar \rho_p,\rho_f,\mu,T antes de cada série; calibrar a régua na frente da coluna (pixels \rightarrow metros).
  2. preparar suspensão para o valor desejado de \varphi usando a fórmula do Passo 3; homogeneizar com agitador; evitar bolhas.
  3. encher a coluna com fluido até nível marcado; equilibrar temperatura.
  4. injetar delicadamente a suspensão no topo (usar placa perfurada para distribuir) evitando correntes fortes; esperar até que a região de medição esteja livre de turbulência.
  5. posicionar a câmera de modo a capturar uma janela vertical de altura \;L_{view}\; onde se espera regime estabelecido. Ajustar fps para resolver deslocamentos entre frames de alguns pixels.
  6. filmar as partículas enquanto caem pela janela de visualização. Repetir filmagem por vários instantes para obter muitas trajetórias (ideal: \gtrsim 20 trajetórias independentes por condição).
  7. processar imagens: aplicar filtro para realçar partículas, localizar centróides em cada frame, ligar pontos por algoritmo de tracking e extrair trajetórias z(t).
  8. para cada trajetória, identificar região de velocidade constante (estacionária) — por exemplo, realizar ajuste linear por partes em z(t) e escolher a região onde aceleração média é próxima de zero. O coeficiente angular da regressão linear dá V.
  9. repetir para todas as partículas e calcular média e desvio padrão de V em cada condição ( \varphi,\Lambda,d).

Passo 8

Como converter trajetórias em velocidade e avaliação de incerteza (detalhado):
  • conversão pixel \rightarrow metro: calibrar com régua; fator s em \mathrm{m/pixel}. Se z_{pix}(t) é posição em pixels, então z(t)=s\,z_{pix}(t).
  • cálculo de velocidade média na região estabelecida: ajustar reta z(t)=Vt+z_0 por mínimos quadrados; V é estimador. A incerteza u_V da regressão vem do erro padrão do coeficiente angular.
  • estimar erro total incluindo: erro de calibração u_s, erro de localização do centróide u_{loc}, flutuações hidrodinâmicas \sigma_{hyd}. Combine por raiz quadrada somada:
u_\text{total} = \sqrt{u_V^2 + \left(\frac{u_s}{s} V\right)^2 + \sigma_{hyd}^2 }.
  • quando ajustar a lei de Richardson–Zaki, usar pesos segundo incerteza em V e reportar incerteza em n (intervalo de confiança do ajuste).

Passo 9

Plano experimental (matriz de ensaios recomendada — valores práticos):
  • partículas: \;d = 0{,}2,\;0{,}5,\;1{,}0\\ \mathrm{mm}.
  • concentrações: \;\varphi = 0,\;0{,}001,\;0{,}01,\;0{,}05,\;0{,}10,\;0{,}20,\;0{,}35. (inclui diluído e regimes semidensos).
  • confinamento: \;\Lambda = D_c/d = 5,\;10,\;20,\;50.
  • repetições: pelo menos \;N_p=20 partículas medidas por condição e \;N_r=3 preparos independentes para estatística.
  • assim o número total de trajetórias fica grande o suficiente para estimar média e erro. Ajuste dependendo do tempo disponível.

Passo 10

Interpretação dos resultados e como responder às duas questões pedidas:
  • influência da concentração: plotar \;V(\varphi)/V_0 vs \;1-\varphi\; e ajustar Richardson–Zaki para obter n(\,Re\,). Se a lei falhar em altas \varphi, procurar efeitos de rede de partículas e flutuações coletivas; discutir limites do modelo.
  • influência dos contornos rígidos: para cada \varphi compare curvas V obtidas em diferentes \Lambda. Plotar \;V(\Lambda)/V(\Lambda\to\infty) vs \;d/D_c e ajustar uma lei empírica (linear ou polinomial em d/D_c). Reportar o limiar de \Lambda acima do qual o efeito de parede é desprezível (prático: \Lambda\gtrsim 20).
  • verificar regimes de Re: calcular Re médio; separar resultados em regime de Stokes e regimes com efeitos inerciais. Ajustar modelagem de arrasto conforme necessário (usar correlações de C_d(Re) quando Re não for muito pequeno).

Passo 11

Principais fontes de erro e como mitigá-las (detalhado):
  • polidispersidade de partículas → usar partículas monodispersas ou separar por peneiramento; reportar distribuição de d.
  • correntes convectivas e convecção térmica → controlar temperatura, esperar estabilização; bloquear luz solar.
  • acúmulo de partículas na superfície superior durante injeção → usar método de injeção bem distribuído (placa perfurada).
  • efeitos de imagem: foco, profundidade de campo e sobreposição de partículas → usar iluminação de fundo e limitar concentração para rastreamento claro; usar técnicas de correlação se imagem estiver densa.
  • mudanças locais de concentração durante a medida → medir somente no interior do “jato” com homogeneidade garantida ou usar média temporal curta.

Passo 12

Diagnósticos complementares e extensões úteis:
  • medir campo de velocidade do fluido local com PIV para avaliar se há retorno do fluido que altera assentamento; isso ajuda a entender interações hidrodinâmicas a \varphi maiores.
  • medir perfil vertical de concentração \varphi(z,t) para estudar dinâmica de estratificação.
  • usar ultrassom/LDV para amostras opacas.
  • estudar partículas próximas a parede para quantificar perfil de velocidade e comparar com correções teóricas (Faxén/experimentos conhecidos).

Passo 13

Check-list prático de implementação (rápido) e recomendações finais:
  • antes de começar: calibrar câmera e medir \rho_f,\mu,T.
  • comece por casos simples: \varphi\approx 0 e \Lambda\gtrsim 50 para obter V_0 experimental e checar consistência com Stokes.
  • depois varie \varphi mantendo \Lambda grande para extrair n.
  • finalmente faça variações de \Lambda para quantificar efeito de parede.
    Pequena nota de ânimo: experimento tem muitas etapas, mas seguindo a matriz e os cuidados de homogeneização e controle térmico você conseguirá obter curvas claras e ajustar as leis empíricas; faça uma planilha para cada condição e guarde vídeos brutos para reprojetar análises.

Resposta

Projeto experimental: coluna cilíndrica transparente com colunas intercambiáveis ( \;D_c/d = 5,\;10,\;20,\;50\;), câmera de alta velocidade, iluminação de fundo, controle de temperatura e sistema de preparação de suspensão. Medir trajetórias z(t) de partículas esféricas ( \;d=0{,}2,\;0{,}5,\;1{,}0\\ \mathrm{mm}\;) para concentrações \;\varphi = 0,\;0{,}001,\;0{,}01,\;0{,}05,\;0{,}10,\;0{,}20,\;0{,}35\;. Calcular velocidade terminal por ajuste linear em região estabelecida e incertezas por erro padrão da regressão. Analisar: (i) influência de \varphi ajustando Richardson–Zaki \;V(\varphi)=V_0(1-\varphi)^n\; para obter n; (ii) influência de paredes medindo V(\Lambda)/V_{\infty} vs \;d/D_c\; e ajustando correção empírica (por exemplo polinomial em \;d/D_c). Verificar regime de validade com \;Re=\rho_f V d/\mu\; e usar correlações de C_d(Re) quando Re\not\ll 1. Reportar matriz experimental, médias, desvios e ajustes com intervalos de confiança.

Exercícios relacionados