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

Resumen Técnico: Distribución y transporte a escala de poros de partículas activas en una red bidimensional

Planteamiento del Problema
Los micronadadores móviles, como las bacterias y los coloides activos sintéticos, habitan frecuentemente entornos porosos complejos caracterizados por obstáculos sólidos y flujos de cizalla heterogéneos. Aunque la acumulación de partículas activas en las fronteras (acumulación superficial) y su respuesta a la cizalla (reotaxia) están bien documentadas en geometrías simples como canales, los mecanismos de transporte en medios porosos estructurados siguen siendo menos comprendidos. Específicamente, existe una brecha en la comprensión de cómo los parámetros gobernantes —actividad de la partícula, porosidad del medio, polidispersidad y fuerza del flujo de fondo— influyen colectivamente en la evolución temporal y las estadísticas de estado estacionario de suspensiones activas en redes ordenadas. Además, se requiere una investigación sistemática sobre el papel de los defectos topológicos en el campo de polarización durante la transición entre los regímenes de natación contra la corriente (upstream) y a favor de la corriente (downstream) en tales entornos de cizalla complejos.

Metodología
Los autores emplean un enfoque de modelado dual para estudiar una suspensión diluida de Partículas de Brownian Activas (ABP) en una red cuadrada bidimensional de pilares sólidos.

1. Modelado Microscópico: Las partículas se modelan como entidades delgadas, de tamaño puntual y no interactuantes que experimentan autopropulsión a una velocidad $v$ y una difusión rotacional $D_R$. Su dinámica está gobernada por ecuaciones de Langevin sobreamortuadas, incorporando ruido traslacional y rotacional.
2. Modelado Continuo: Se deriva una descripción macroscópica utilizando una ecuación de Fokker-Planck para la función de distribución de probabilidad. Esta se cierra mediante una aproximación de dos momentos (negando el orden nemático, $Q=0$) para resolver la evolución de los campos de densidad numérica ($n$) y polarización ($m$).
3. Simulaciones Numéricas:
   • Pilar Único: Las ecuaciones cinéticas continuas se resuelven numéricamente (Método de Líneas de la Forma) para un pilar aislado en un dominio no acotado para establecer comportamientos base de acumulación y polarización.
   • Red Periódica: Se realizan simulaciones de Dinámica de Brownian (BD) para un gran número de partículas en una red cuadrada doblemente periódica. Este enfoque permite el análisis de geometrías complejas y múltiples pilares donde los cierres continuos pueden fallar.
   • Condiciones de Flujo: Se impone un flujo de fondo impulsado por presión, generando un campo de flujo de Stokes bidimensional computado mediante el Método de Soluciones Fundamentales. La dinámica de las partículas se modifica para incluir la advección del flujo y la ecuación de Jeffery para la orientación en cizalla.

Los parámetros adimensionales clave incluyen el número de Péclet de nado ($Pe_s = v/aD_R$), el número de Péclet de flujo ($Pe_f = u_\infty/aD_R$) y la porosidad del medio ($\epsilon$).

Contribuciones Clave y Resultados

• Acumulación Superficial en Medios Estáticos: En ausencia de flujo, la autopropulsión impulsa a las partículas a acumularse en las superficies de los pilares, creando una capa límite de densidad numérica aumentada y polarización radial (partículas alineándose hacia la superficie). El sistema evoluciona desde un estado homogéico inicial hacia un estado estacionario mediante crecimiento exponencial.
• Oscilaciones Transitorias: En redes confinadas (porosidad $\epsilon < 1$), el acercamiento al estado estacionario se caracteriza por oscilaciones transitorias amortiguadas en la densidad numérica superficial. Estas se atribuyen al movimiento balístico entre pilares, lo que crea ondas de agotamiento que se propagan e interactúan. La frecuencia de estas oscilaciones escala inversamente con la longitud del dominio.
• Efecto de la Porosidad y la Actividad:
  • Una menor porosidad (mayor confinamiento) conduce a tiempos de equilibración más rápidos debido a distancias de viaje más cortas, pero reduce la relación de acumulación superficial en estado estacionario ($\eta = \bar{n}_A/\bar{n}_B$) para una actividad fuerte.
  • Una mayor actividad ($Pe_s$) aumenta la acumulación en estado estacionario y agudiza los gradientes.
  • En redes polidispersas (pilares de dos tamaños), los pilares más pequeños exhiben tasas de acumulación inicial más rápidas pero densidades de estado estacionario más bajas en comparación con los pilares más grandes debido a los efectos de agotamiento de los pilares mayores. Se encuentra que los medios polidispersos son más efectivos para filtrar nadadores del fluido de fondo en comparación con los medios monodispersos de la misma porosidad.
• Dinámica Inducida por el Flujo y Defectos Topológicos:
  • Reotaxia: En presencia de flujo de fondo, las partículas exhiben reotaxia positiva (natación contra la corriente) cerca de las superficies de los pilares, lo que conduce a una acumulación preferencial en la estela del pilar.
  • Transición de Natación: El sistema experimenta una transición de una polarización global contra la corriente a bajos niveles de fuerza de flujo ($Pe_f$) hacia una coexistencia de regiones de natación contra la corriente y a favor de la corriente a altos niveles de $Pe_f$.
  • Defectos Topológicos: Esta transición es mediada por la nucleación de defectos topológicos puramente cinemáticos en el campo de polarización. A un $Pe_f$ moderado, emergen pares de defectos espirales (carga +1) y defectos hiperbólicos (carga -1). Estos defectos son estacionarios en el marco de la red y surgen de la naturaleza bidimensional del flujo de cizalla, distinguiéndose de los defectos en fluidos polares activos impulsados por inestabilidades hidrodinámicas.
  • Fuerza de Nado: La distribución no uniforme de partículas genera una fuerza de nado neta sobre los pilares. Esta fuerza alcanza su máximo en fuerzas de flujo intermedias antes de disminuir conforme la advección agota la superficie.

Significancia
El artículo proporciona un marco controlado para aislar los mecanismos físicos que gobiernan el transporte activo en geometrías complejas. Demuestra que el transporte de coloides activos en medios porosos es fundamentalmente distinto del transporte pasivo, impulsado por la interacción entre la autopropulsión y perfiles de cizalla complejos. Un hallazgo principal es la identificación de los defectos topológicos cinemáticos como el mecanismo que facilita la transición entre los regímenes de natación contra la corriente y a favor de la corriente, un fenómeno específico del entorno de cizalla bidimensional de la red. La red estructurada sirve como un sistema modelo para comprender el transporte en entornos microfluídicos estructurados, ofreciendo información sobre cómo el confinamiento, la actividad y la fuerza del flujo dictan la distribución espacial, el comportamiento colectivo y la distribución de la materia activa. El estudio destaca que, incluso en suspensiones diluidas y no interactuantes, los flujos de cizalla complejos pueden inducir comportamientos dinámicos ricos, incluyendo la nucleación de defectos y oscilaciones transitorias, que son críticos para comprender la colonización biológica y el transporte sintético en entornos porosos.