Dissipation and microstructure in sheared active suspensions of squirmers

Mediante simulaciones de dinámica de Stokes activa, este estudio revela que el cizallamiento en suspensiones de squirmers apolares aumenta la disipación energética y reduce la viscosidad relativa mientras induce un orden nemático y correlaciones de pares anisotrópicas, demostrando que la actividad interna, más que la motilidad, es el factor determinante en la respuesta reológica y microestructural de estos sistemas.

Lo esencial

¡Hola! Imagina que tienes un vaso de agua, pero en lugar de solo agua, está lleno de millones de pequeños "robots" microscópicos que pueden moverse por sí mismos. A esto los científicos lo llaman suspensión activa.

Este artículo es como un informe de laboratorio que estudia qué pasa cuando agitamos (cortan) este vaso de robots con una cuchara gigante. Los autores, Zhouyang Ge y Gwynn Elfring, quieren entender dos cosas:

1. ¿Cuánta energía se gasta? (¿Se cansa el sistema?).
2. ¿Qué tan espeso se vuelve el líquido? (¿Se comporta como miel o como agua?).

Aquí te explico los hallazgos principales usando analogías sencillas:

1. Los protagonistas: Los "Sacudidores" (Shakers)

En lugar de estudiar bacterias reales (que son complejas), los autores crearon modelos matemáticos de esferas perfectas llamadas "squirmers".

• Imagina que tienes dos tipos de robots:
  • Los "Empujadores" (Pushers): Como un barco que empuja el agua hacia atrás para avanzar.
  • Los "Tiradores" (Pullers): Como un barco que jala el agua hacia adelante para avanzar.
• El truco: En este estudio, estos robots no pueden nadar solos. Están "paralizados" individualmente. Solo pueden moverse si se empujan entre ellos o si alguien los agita. Son como una multitud de personas en una habitación que no pueden caminar, pero si se empujan unas a otras, la multitud se mueve.

2. La sorpresa: Agitar hace que sea más "fluido" (Shear-thinning)

Lo normal en la vida diaria es que si tienes una mezcla espesa (como pintura o ketchup) y la agitas fuerte, se vuelve más líquida. Esto se llama adelgazamiento por cizalla.

• Lo que esperaban: Pensaban que, como estos robots son "activos" (tienen energía propia), quizás se comportarían como las bacterias reales que a veces hacen que el líquido sea menos viscoso (como un "superfluido" mágico) o más espeso.
• Lo que descubrieron: ¡Agitarlos hace que la mezcla sea menos viscosa (más fácil de mover) para ambos tipos de robots!
  • La analogía: Imagina una habitación llena de gente intentando cruzar. Si todos están quietos, chocan y se atascan (espeso). Si alguien empieza a agitar la multitud (cortar), la gente se organiza mejor y fluye más rápido, aunque sigan chocando.

3. El gasto de energía: Más agitación = Más gasto

Aunque el líquido se vuelve más fácil de mover (menos viscoso), el sistema gasta mucha más energía cuando lo agitas.

• La analogía: Piensa en una fiesta. Si la gente está quieta, se cansa poco. Si pones música y empiezan a bailar (agitar), aunque se muevan con más fluidez por la pista, el consumo de energía (cansancio) de todos los invitados se dispara.
• Los robots "empujadores" gastan más energía que los "tiradores" cuando la mezcla está muy llena, lo cual es lo contrario de lo que ocurre en algunas bacterias reales.

4. El secreto: La "Boda" de los robots (Microestructura)

¿Por qué pasa esto? Los autores miraron cómo se organizan los robots cuando los agitan. Descubrieron dos cosas extrañas:

• Orden en el caos (Orden nemático): Cuando agitan la mezcla a una velocidad media, los robots dejan de mirar en todas direcciones y empiezan a alinearse, como si formaran filas o parejas.
  • La analogía: Imagina un tráfico caótico. De repente, todos los coches deciden ir en la misma dirección y se alinean perfectamente. Esto crea una estructura que, irónicamente, hace que el líquido se comporte de manera extraña.
• Parejas que chocan: Los robots tienden a formar "duetos" o parejas que chocan de frente o de lado.
  • La analogía: Es como si dos personas en una multitud se agarraran de las manos y empezaran a girar juntas. Estas "parejas" crean un flujo de agua que ayuda a empujar a los demás, facilitando el movimiento general (reduciendo la viscosidad), pero al mismo tiempo, el esfuerzo de mantenerse unidos gasta mucha energía.

5. La lección final: La actividad es más importante que el movimiento

Hicieron una comparación con robots que sí pueden nadar solos (como las bacterias reales).

• El hallazgo clave: Los robots que no nadan pero tienen "actividad interna" (se empujan entre sí) mostraron un efecto mucho más fuerte de reducción de viscosidad que los que solo nadan solos.
• La moraleja: No es solo que los robots se muevan; es que su interacción (cómo se empujan y chocan) es lo que realmente cambia las reglas del juego. La "actividad" (la energía interna) es más importante que la "motilidad" (la capacidad de nadar).

En resumen

Este estudio nos dice que si tienes una mezcla de micro-robots que se empujan entre sí y la agitas:

1. Se volverá más líquida y fácil de mover (como el ketchup).
2. Pero gastará mucha más energía en el proceso.
3. Todo esto sucede porque los robots se organizan en patrones extraños (parejas y filas) debido a cómo se empujan, no porque sepan nadar.

Es un recordatorio de que en el mundo microscópico, cómo interactúas con tus vecinos es tan importante como cómo te mueves tú mismo.

Resumen técnico

Aquí presento un resumen técnico detallado del artículo "Dissipation and microstructure in sheared active suspensions of squirmers" (Disipación y microestructura en suspensiones activas de squirmers cortados), basado en el documento proporcionado.

1. Planteamiento del Problema

Las suspensiones activas difieren de los fluidos complejos pasivos porque sus partículas suspendidas pueden gastar y disipar energía internamente, lo que genera motilidad o tensiones activas. Un fenómeno notable en ciertas suspensiones bacterianas (como E. coli) es la reducción de la viscosidad con el aumento de la concentración de partículas, llegando incluso a una transición de "superfluido" a bajas tasas de cizalla. Sin embargo, otros sistemas, como las suspensiones de algas, muestran un aumento de la viscosidad.

Existe una brecha en la comprensión teórica y numérica de estos fenómenos en regímenes semi-diluidos a concentrados (donde las interacciones entre partículas son significativas). La mayoría de las teorías cinéticas asumen regímenes diluidos o no consideran adecuadamente las interacciones hidrodinámicas en 3D. El objetivo de este estudio es elucidar cómo la actividad interna y el flujo externo (cizalla) interactúan para determinar la disipación de energía y la microestructura en suspensiones de partículas activas esferoidales (squirmers), diferenciando entre efectos de motilidad (natación) y efectos de actividad (fuerzas hidrodinámicas).

2. Metodología

Los autores emplean Dinámica de Stokes Rápida Activa (Active Fast Stokesian Dynamics - FSD) para simular suspensiones de squirmers en un flujo de cizalla simple.

• Modelos de Partículas:
  • Shakers (Agitadores): Squirmers apolares que son individualmente inmóviles ($B_1=0$) pero generan tensiones activas ($B_2 \neq 0$). Se estudian dos tipos: Pullers (traccionadores, $B_2 > 0$, análogos a algas) y Pushers (empujadores, $B_2 < 0$, análogos a bacterias). A pesar de no nadar individualmente, en suspensión generan difusión hidrodinámica colectiva.
  • ABPs (Partículas Brownianas Activas): Squirmers neutros con auto-propulsión ($U_s \neq 0$) pero sin interacciones hidrodinámicas activas de primer orden ($B_2 = 0$). Se utilizan para aislar el efecto de la motilidad pura.
• Parámetros de Simulación:
  • Se simulan $N = 1024$ partículas (con pruebas de tamaño hasta 8192) en una caja periódica.
  • Fracciones de volumen ($\phi$): 10%, 20%, 30%, 40%.
  • Números de Péclet ($Pe$): Varían desde regímenes dominados por la actividad ($Pe < 1$) hasta dominados por la cizalla ($Pe \gg 100$).
• Cálculos: Se calculan el tensor de esfuerzos, la disipación de energía (interna y externa), la viscosidad relativa, las diferencias de esfuerzos normales y las funciones de correlación de pares para caracterizar la microestructura.

3. Contribuciones Clave y Resultados

A. Disipación de Energía

• Efecto de la Cizalla: La cizalla aumenta la disipación total de energía en comparación con el estado sin cizalla.
• Regímenes:
  • $Pe < 1$ (Dominio de Actividad): La disipación es casi constante e independiente de la tasa de cizalla. Curiosamente, en este régimen, los pushers disipan más energía que los pullers a fracciones de volumen altas, lo cual es opuesto a la intuición basada en la difusión hidrodinámica (donde los pullers son más difusivos).
  • $Pe \gg 1$ (Dominio de Cizalla): La disipación total escala con $Pe^2$, comportándose como esferas pasivas.
• Interacciones: Las interacciones hidrodinámicas aumentan la disipación, mientras que las repulsiones de corto alcance (fuerzas estéricas) tienden a reducirla ligeramente al mantener las partículas separadas.

B. Reología (Viscosidad y Esfuerzos)

• Comportamiento de Corte (Shear-Thinning): Contrario a la idea de que los pushers espesan el fluido y los pullers lo adelgazan (común en teorías de partículas alargadas), este estudio muestra que tanto pullers como pushers esféricos exhiben adelgazamiento por cizalla (reducción de la viscosidad relativa $\eta_r$ al aumentar $Pe$).
• Efecto de la Concentración: A diferencia de las suspensiones bacterianas que muestran una transición a superfluido (viscosidad que disminuye con $\phi$), las suspensiones de shakers esféricos muestran un aumento de la viscosidad con la fracción de volumen ($\phi$), similar a suspensiones pasivas pero con una viscosidad intrínseca más alta debido a la actividad.
• Motilidad vs. Actividad: Las simulaciones de ABPs (solo motilidad, sin interacciones hidrodinámicas activas) muestran un adelgazamiento por cizalla mucho más débil. Esto demuestra que la actividad hidrodinámica (tensiones dipolares) es el motor principal de los cambios reológicos, no la mera capacidad de nadar.
• Esfuerzos Normales: Se observan diferencias de esfuerzos normales negativas ($N_1, N_2 < 0$) a altos $Pe$, indicativo de una pérdida de isotropía y alineación inducida por colisiones, similar a suspensiones pasivas.

C. Microestructura

El estudio revela firmas microestructurales inusuales que explican la reología:

• Orden Nemático: Existe un orden nemático global (alineación de direcciones) que es máximo en el rango intermedio $1 \lesssim Pe \lesssim 100$. Este orden es no monótono: es bajo a $Pe$ muy bajos (isotrópico) y muy altos (aleatorizado por cizalla).
• Correlación de Pares:
  • En el régimen de transición ($Pe \approx 10$), se observa una fuerte anisotropía en la distribución de pares de partículas.
  • Se forman dúos (doublets) donde las partículas se alinean de manera que sus flujos inducidos por actividad se oponen al flujo de cizalla externo.
  • Para pullers, la alineación tiende a ser en el cuadrante de extensión; para pushers, en el de compresión.
• Mecanismo de Viscosidad: Esta alineación específica y la formación de estructuras de pares que resisten el flujo son las responsables del aumento de la viscosidad a bajas tasas de cizalla y del comportamiento de adelgazamiento por cizalla.

4. Significado e Implicaciones

• Desafío a Teorías Existentes: Los resultados contradicen predicciones de teorías cinéticas que ignoran interacciones hidrodinámicas (como el modelo de Takatori y Brady), las cuales predecían un engrosamiento por cizalla o viscosidades negativas para ciertos parámetros. El estudio demuestra que las interacciones hidrodinámicas son cruciales para la reología correcta.
• Distinción Conceptual: El trabajo establece claramente la diferencia entre motilidad (movimiento auto-propulsado) y actividad (generación de tensiones hidrodinámicas). En sistemas esféricos, la actividad hidrodinámica es el factor dominante en la reología, no la capacidad de natación individual.
• Física de Materia Activa: Proporciona una explicación mecanicista basada en la microestructura (orden nemático y correlaciones de pares) para el comportamiento reológico complejo en suspensiones concentradas, sugiriendo que la transición a "superfluido" observada en bacterias podría depender fuertemente de la forma de las partículas (alargadas) y no solo de su actividad, ya que las esferas activas no muestran dicha transición en este régimen.

En resumen, el artículo demuestra que en suspensiones concentradas de partículas activas esféricas, la interacción entre la actividad interna y el flujo externo genera un orden nemático transitorio y correlaciones de pares anisotrópicas que aumentan la disipación y la viscosidad a bajas tasas de cizalla, resultando en un comportamiento de adelgazamiento por cizalla para ambos tipos de dipolos (traccionadores y empujadores).