Percolation of a cohesive fine particle in a static bed

Resumo Técnico: Percolação de uma Partícula Fina Coesiva em um Leito Estático

Enunciado do Problema
A percolação de partículas finas (finos) através de um leito estático de partículas maiores é um processo fundamental no manuseio industrial (por exemplo, farmacêuticos, processamento químico) e em fluxos geofísicos (por exemplo, transporte de sedimentos, contaminação de águas subterrâneas). Em sistemas não coesivos, se uma partícula fina passa ou fica presa é governado principalmente por restrições geométricas, especificamente a razão entre os diâmetros das partículas grandes e finas ($R = d_l/d_f$). Existe um limiar geométrico de aprisionamento em $R_t \approx 6,464$; quando $R > R_t$, espera-se teoricamente que os finos passem livremente através dos gargalos de poros formados pelas partículas grandes. No entanto, sistemas do mundo real frequentemente envolvem forças coesivas (van der Waals, eletrostáticas, pontes líquidas) que alteram a dinâmica de contato. Embora estudos anteriores tenham abordado a percolação não coesiva e a reologia de fluxos coesivos, há uma lacuna na compreensão de como a coesão influencia a percolação no regime de "peneiramento livre", onde restrições geométricas sozinhas permitiriam a passagem. Este estudo investiga a percolação de partículas finas coesivas através de um leito estático de partículas grandes onde $R=7$ (acima do limiar geométrico), focando no limite de partícula única onde os finos interagem com as partículas do leito, mas não entre si.

Metodologia
Os autores empregam simulações pelo Método dos Elementos Discretos (DEM) utilizando o software de código aberto LAMMPS. As simulações utilizam um esquema Velocity-Verlet aprimorado adequado para grandes razões de tamanho.

• Configuração do Sistema: Um leito estático de partículas grandes monodispersas ($d_l = 4$ mm) é compactado com uma fração $\phi = 0,6$. Partículas finas ($d_f \approx 0,57$ mm, $R=7$) são liberadas do repouso acima do leito e permitem cair sob gravidade.
• Modelos de Contato: As forças de contato normais seguem o modelo elástico de Hertz. A inelasticidade é modelada via um coeficiente de restituição constante ($e_n$). As forças tangenciais utilizam o modelo de Mindlin com atrito de deslizamento ($\mu$) e resistência ao rolamento.
• Coesão: Interações coesivas são modeladas principalmente usando o modelo Derjaguin-Muller-Toporov (DMT), que assume uma força atrativa constante durante o contato e repulsão de Hertz. A força da coesão é quantificada pelo Número de Bond ($Bo = F_c / W$), definido como a razão entre a força de arrancamento e o peso da partícula fina. Um número limitado de simulações usando o modelo Johnson-Kendall-Roberts (JKR) também foi realizado para comparação.
• Parâmetros: O estudo varia sistematicamente o Número de Bond ($Bo$), o coeficiente de atrito de deslizamento ($\mu$), o coeficiente de restituição ($e_n$) e o Módulo de Young ($E$). A dinâmica de $10^4$ partículas finas independentes é rastreada para garantir precisão estatística.

Principais Resultados

1. Aprisionamento Não Geométrico: Mesmo quando $R=7$ (permitindo geometricamente a passagem livre), coesão e atrito suficientemente grandes levam ao aprisionamento de partículas. À medida que $Bo$ aumenta, a distância característica de percolação diminui, significando que as partículas ficam presas em profundidades mais rasas.
2. Distribuição de Aprisionamento: A distribuição vertical de partículas presas segue um decaimento exponencial. A distância característica de aprisionamento ($\lambda$) serve como uma métrica para o comportamento de percolação. A distância inversa de aprisionamento ($d_l/\lambda$) aumenta monotonicamente com $Bo$.
3. Papel do Atrito e da Coesão:
   • Em baixo $Bo$, o atrito de deslizamento desempenha um papel significativo no aprisionamento.
   • Em alto $Bo$, a coesão torna-se o mecanismo de aprisionamento dominante.
   • Um "joelho" é observado na relação entre a distância inversa de aprisionamento e o atrito. Este joelho se desloca com $Bo$ e corresponde a uma transição entre regimes onde o atrito domina e onde a coesão domina.
   • A rigidez da partícula ($E$) influencia o aprisionamento: partículas mais macias (menor $E$) são mais facilmente presas, particularmente em $Bo$ mais altos.
4. Mecanismo Colisional: O processo de aprisionamento é decomposto em duas etapas sequenciais:
   • Condição de Adesão: Após a colisão, uma partícula ou ricocheteia ou adere. A adesão ocorre se a velocidade de impacto estiver abaixo de uma velocidade crítica ($v_{cr}$), que depende de $Bo$, $e_n$ e $E$.
   • Movimento Pós-Contato: Se a partícula aderir, ela pode deslizar e se desprender ou permanecer presa. Uma condição de equilíbrio estático baseada no balanço de forças determina se uma partícula que não ricocheteia permanece estacionária. Esta condição depende de um coeficiente de atrito reescalado $\tilde{\mu} = \mu\sqrt{4Bo^2 - 1}$. Se $\tilde{\mu} > 1$, a partícula provavelmente permanecerá estacionária; se $\tilde{\mu} < 1$, ela provavelmente deslizará e potencialmente se desprenderá.
5. Modelo Preditivo: Os autores desenvolvem um modelo colisional ligando a dinâmica microscópica à distância macroscópica de aprisionamento. A distância inversa de aprisionamento é prevista por $1/\lambda = -k \ln(1 - P_t)$, onde $k$ é a frequência de colisão e $P_t$ é a probabilidade de aprisionamento por colisão (determinada pela razão entre a velocidade de impacto e $v_{cr}$). Este modelo prevê com precisão as distâncias de aprisionamento no regime dominado por colisões (alto $Bo$).
6. Duplo Contato: Em regimes com menor atrito ou coesão, as partículas frequentemente deslizam e formam contatos duplos (tocando duas partículas grandes simultaneamente). Esta configuração pode levar ao aprisionamento mesmo quando condições de contato único sugerissem desprendimento, ou, inversamente, permitir o desprendimento via um ricochete de um segundo contato.

Significância e Alegações
O artigo estabelece que a percolação de partículas finas coesivas não é determinada exclusivamente pela acessibilidade geométrica. Em vez disso, a dinâmica de interação em escala de partícula — especificamente a inter-relação de coesão, restituição, atrito e velocidade de impacto — pode sobrepor expectativas geométricas, levando ao aprisionamento mesmo quando os gargalos de poros são geometricamente largos o suficiente para a passagem.

Os autores afirmam fornecer uma estrutura fundamental para a compreensão desses mecanismos em leitos estáticos, o que é um pré-requisito para abordar o problema mais complexo de segregação granular em misturas em fluxo de partículas coesivas. Ao identificar as condições específicas (velocidade crítica, critérios de equilíbrio estático) que governam a adesão e o movimento subsequente, o estudo oferece uma capacidade preditiva para a distância de aprisionamento baseada em propriedades materiais fundamentais. Os resultados destacam que, em sistemas coesivos, o "destino" de uma partícula fina é governado por uma sequência de interações locais (colisão, adesão, movimento pós-contato) e não apenas pela geometria estática do leito.