Pore-scale distribution and transport of active particles in a two-dimensional lattice

Resumo Técnico: Distribuição e transporte de partículas ativas em escala de poros em uma rede bidimensional

Definição do Problema
Micro-nadadores móveis, como bactérias e coloides ativos sintéticos, habitam frequentemente ambientes porosos complexos caracterizados por obstáculos sólidos e fluxos de cisalhamento heterogêneos. Embora o acúmulo de partículas ativas em fronteiras (acúmulo de superfície) e sua resposta ao cisalhamento (reotaxia) sejam bem documentados em geometrias simples como canais, os mecanismos de transporte em meios porosos estruturados permanecem menos compreendidos. Especificamente, há uma lacuna na compreensão de como os parâmetros governantes — atividade da partícula, porosidade do meio, polidispersidade e força do fluxo de fundo — influenciam coletivamente a evolução temporal e as estatísticas de estado estacionário de suspensões ativas em redes ordenadas. Além disso, o papel dos defeitos topológicos no campo de polarização durante a transição entre os regimes de nado a montante (upstream) e a jusante (downstream) em tais ambientes de cisalhamento complexos requer investigação sistemática.

Metodologia
Os autores empregam uma abordagem de modelagem dupla para estudar uma suspensão diluída de Partículas Brownianas Ativas (ABPs) em uma rede quadrada bidimensional de pilares sólidos.

1. Modelagem Microscópica: As partículas são modeladas como entidades pontuais, delgadas e não interagentes que sofrem autopropulsão a uma velocidade $v$ e difusão rotacional $D_R$. Sua dinâmica é governada por equações de Langevin superamortecidas, incorporando ruído translacional e rotacional.
2. Modelagem Contínua: Uma descrição macroscópica é derivada usando uma equação de Fokker-Planck para a função de distribuição de probabilidade. Esta é fechada usando uma aproximação de dois momentos (negligenciando a ordem nemática, $Q=0$) para resolver a evolução dos campos de densidade numérica ($n$) e polarização ($m$).
3. Simulações Numéricas:
   • Pilar Único: As equações cinéticas contínuas são resolvidas numericamente (Método das Linhas) para um pilar isolado em um domínio ilimitado para estabelecer comportamentos de base de acúmulo e polarização.
   • Rede Periódica: Simulações de Dinâmica Browniana (BD) são realizadas para um grande número de partículas em uma rede quadrada duplamente periódica. Esta abordagem permite a análise de geometrias complexas e múltiplos pilares onde os fechamentos contínuos podem falhar.
   • Condições de Fluxo: Um fluxo de fundo impulsionado por pressão é imposto, gerando um campo de fluxo de Stokes bidimensional computado via Método das Soluções Fundamentais. A dinâmica das partículas é modificada para incluir advecção de fluxo e a equação de Jeffery para orientação em cisalhamento.

Os principais parâmetros adimensionais incluem o número de Péclet de nado ($Pe_s = v/aD_R$), o número de Péclet de fluxo ($Pe_f = u_\infty/aD_R$) e a porosidade do meio ($\epsilon$).

Principais Contribuições e Resultados

• Acúmulo de Superfície em Meios Estáticos: Na ausência de fluxo, a autopropulsão leva as partículas a se acumularem nas superfícies dos pilares, criando uma camada limite de aumento de densidade numérica e polarização radial (partículas alinhando-se em direção à superfície). O sistema evolui de um estado inicial homogêneo para um estado estacionário via crescimento exponencial.
• Oscilações Transientes: Em redes confinadas (porosidade $\epsilon < 1$), a aproximação ao estado estacionário é caracterizada por oscilações transientes e amortecidas na densidade numérica de superfície. Estas são atribuídas ao movimento balístico entre os pilares, criando ondas de depleção que se propagam e interagem. A frequência dessas oscilações escala inversamente com o comprimento do domínio.
• Efeito da Porosidade e Atividade:
  • Menor porosidade (maior confinamento) leva a tempos de equilíbrio mais rápidos devido a distâncias de viagem mais curtas, mas reduz a razão de acúmulo de superfície em estado estacionário ($\eta = \bar{n}_A/\bar{n}_B$) para atividade forte.
  • Maior atividade ($Pe_s$) aumenta o acúmulo em estado estacionário e intensifica os gradientes.
  • Em redes polidispersas (pilares de dois tamanhos), pilares menores exibem taxas iniciais de acúmulo mais rápidas, mas densidades de estado estacionário menores comparadas a pilares maiores devido aos efeitos de depleção dos pilares maiores. Constatou-se que meios polidispersos são mais eficazes em filtrar nadadores do fluido de fundo em comparação com meios monodispersos de mesma porosidade.
• Dinâmica Induzida por Fluxo e Defeitos Topológicos:
  • Reotaxia: Na presença de fluxo de fundo, as partículas exibem reotaxia positiva (nado a montante) próximo às superfícies dos pilares, levando ao acúmulo preferencial na esteira do pilar.
  • Transição de Nado: O sistema passa por uma transição de uma polarização global a montante em baixas forças de fluxo ($Pe_f$) para uma coexistência de regiões de nado a montante e a jusante em altas forças de fluxo ($Pe_f$).
  • Defeitos Topológicos: Esta transição é mediada pela nucleação de defeitos topológicos puramente cinemáticos no campo de polarização. Em $Pe_f$ moderado, pares de defeitos espirais (carga +1) e hiperbólicos (carga -1) emergem. Estes defeitos são estacionários no referencial da rede e surgem da natureza bidimensional do fluxo de cisalhamento, distinguindo-se de defeitos em fluidos polares ativos impulsionados por instabilidades hidrodinâmicas.
  • Força de Nado: A distribuição não uniforme de partículas gera uma força de nado líquida sobre os pilares. Esta força atinge o pico em forças de fluxo intermediárias antes de diminuir conforme a advecção depletar a superfície.

Significância
O artigo fornece um quadro controlado para isolar os mecanismos físicos que regem o transporte ativo em geometrias complexas. Demonstra que o transporte de coloides ativos em meios porosos é fundamentalmente distinto do transporte passivo, impulsionado pela interação entre autopropulsão e perfis de cisalhamento complexos. Uma descoberta primária é a identificação de defeitos topológicos cinemáticos como o mecanismo que facilita a transição entre os regimes de nado a montante e a jusante, um fenômeno específico do ambiente de cisalhamento bidimensional da rede. A rede estruturada serve como um sistema modelo para compreender o transporte em configurações microfluídicas estruturadas, oferecendo insights sobre como o confinamento, a atividade e a força do fluxo ditam a distribuição espacial e o comportamento coletivo da matéria ativa. O estudo destaca que mesmo em suspensões diluídas e não interagentes, fluxos de cisalhamento complexos podem induzir comportamentos dinâmicos ricos, incluindo nucleação de defeitos e oscilações transientes, que são críticos para entender a colonização biológica e o transporte sintético em ambientes porosos.