🔬 condensed matter

Rheology of dense suspensions of granular spherocylinders by particle-based simulation

Este artículo presenta un modelo de simulación basado en partículas que predice con éxito la viscosidad y la microestructura de suspensiones densas de esferocilindros granulares bajo flujo de cizalla, revelando cómo estas propiedades evolucionan con la fracción de volumen y la relación de aspecto al tiempo que corrobora datos experimentales limitados.

Autores originales: Alex Dixon, John Hone, Gavin Melaugh, Christopher Ness

Publicado 2026-07-31

📖 4 min de lectura☕ Lectura para el café

CC BY 4.0

Autores originales: Alex Dixon, John Hone, Gavin Melaugh, Christopher Ness

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Imagina un mundo donde lo que viertes, bombeas o mezclas no es solo agua o arena, sino una sopa caótica de diminutos palos rígidos flotando en un líquido. Este es el reino de la reología, la ciencia de cómo fluyen las cosas. Probablemente hayas visto esto en acción: piensa en un río obstruido por árboles caídos que forman un enorme "atasco de troncos", o en la densa y pegajosa pasta utilizada para fabricar papel. En estas situaciones, las partículas están tan amontonadas que chocan constantemente entre sí, y su forma importa tanto como cuántas hay. Si intentas agitar un cubo de canicas, estas ruedan unas junto a otras con facilidad. Pero si intentas agitar un cubo de varillas largas y delgadas, se enredan, se traban entre sí y, de repente, toda la mezcla actúa como un sólido. Los científicos han comprendido durante mucho tiempo cómo se comportan las bolas redondas en estas multitudes, pero las reglas para las partículas largas con forma de varilla han sido un poco misteriosas, especialmente cuando están empaquetadas estrechamente. Comprender esto es crucial porque estas mezclas de "palos pegajosos" aparecen en todas partes, desde la fabricación industrial hasta el flujo de roca fundida dentro de los volcanes.

En este artículo, un equipo de investigadores construyó una simulación computacional superinteligente para actuar como un laboratorio virtual para estas suspensiones de varillas amontonadas. En lugar de mezclar productos químicos reales en un laboratorio, crearon un mundo digital lleno de miles de "esferocilindros": partículas que parecen cápsulas o píldoras (un cilindro con una semiesfera en cada extremo). Los situaron flotando en un fluido virtual y luego comenzaron a "cizallarlos", que es una forma elegante de decir que arrastraron la capa superior del fluido lateralmente mientras mantenían la inferior quieta, tal como se unta la mantequilla en una tostada. El objetivo era ver cómo cambiaba la viscosidad de la mezcla a medida que empaquetaban más varillas en el espacio y a medida que hacían las varillas más largas.

La simulación reveló comportamientos fascinantes. Primero, cuando comenzaron el cizallamiento, la mezcla se volvió increíblemente espesa por un breve instante —un "pico de viscosidad"— antes de estabilizarse. Esto sucedió porque las varillas estaban inicialmente desordenadas en direcciones aleatorias, chocando entre sí como en un mosh pit. Una vez que comenzó el flujo, las varillas empezaron a alinearse, como soldados marchando en formación, y la mezcla se volvió más fácil de agitar. Los investigadores descubrieron que cuanto más varillas empaquetaban (la fracción de volumen, $\phi$ ), más espesa se volvía la mezcla, llegando eventualmente a un punto en el que se trababa y dejaba de fluir por completo. Curiosamente, cuanto más largas eran las varillas (la relación de aspecto, $A$ , que va desde 1 para una esfera hasta 20 para palos muy largos), más pronto ocurría este atasco. Por ejemplo, con varillas de longitud 5, el punto de atasco estaba alrededor de una fracción de volumen de 0.57, pero para varillas de longitud 20, la mezcla se trababa a una densidad mucho menor de aproximadamente 0.41.

El equipo también descubrió que las varillas no solo se alineaban; se alineaban específicamente con el flujo. Midieron un "parámetro de orden" ( $S$ ) para ver qué tan bien se alineaban. A densidades bajas, las varillas se tambaleaban y rotaban en bucles (conocidos como órbitas de Jeffery), pero a medida que la multitud se hacía más densa, se bloqueaban en una línea recta a lo largo de la dirección del flujo. Sin embargo, si la multitud se volvía demasiado densa, las varillas empezaban a frustrarse y ya no podían mantenerse perfectamente alineadas, lo que hacía que el orden disminuyera ligeramente. Los investigadores también desglosaron las fuerzas en juego, mostrando que cuando la mezcla era fluida, la propia resistencia del líquido era lo más importante, pero a medida que se amontonaba, los golpes físicos y la fricción entre las varillas se convertían en la fuerza dominante, asumiendo la tarea de hacer que la mezcla fuera espesa.

Al crear un modelo que equilibra la física del contacto, la fricción y la resistencia del fluido, los autores proporcionaron una nueva herramienta para predecir cómo se comportan estas mezclas complejas. No se limitaron a suponer; realizaron miles de simulaciones con diferentes números de partículas, diferentes longitudes de varilla y diferentes densidades de empaquetamiento. Sus resultados sugieren que, si bien la idea básica de que "más materia equivale a un fluido más espeso" se mantiene, la forma de las partículas cambia drásticamente el punto en el que la mezcla se bloquea. Este trabajo no pretende resolver todos los misterios de la mecánica de fluidos, pero ofrece una base simulada sólida para comprender cómo gestionar estos complicados fluidos llenos de varillas en el mundo real, desde la prevención de atascos de troncos en los ríos hasta el diseño de mejores lodos industriales.

Resumen Técnico: Reología de Suspensiones Densas de Esferocilindros Granulares mediante Simulación Basada en Partículas

Planteamiento del Problema
Las suspensiones densas de partículas granulares con forma de varilla son prevalentes en contextos naturales e industriales, que van desde atascos de troncos en ríos hasta suspensiones cristalinas en el magma y redes de fibras en la fabricación de papel. Aunque la reología de las suspensiones esféricas está bien caracterizada, particularmente a través de la unificación de suspensiones no Brownianas y flujos granulares secos por Boyer et al., el comportamiento de las partículas alargadas permanece incompletamente comprendido. El grado de libertad adicional asociado con la orientación de las partículas introduce fenómenos microestructurales complejos, tales como las órbitas de Jeffery y la alineación inducida por cizalladura, que alteran significativamente la viscosidad efectiva y la respuesta de esfuerzo. Los datos experimentales existentes para los regímenes densos (altas fracciones de volumen sólido, $\phi$ ) son limitados debido a dificultades de preparación y medición. Además, los enfoques de simulación previos enfrentan compromisos: la Dinámica de Stokes (SD) es computacionalmente costosa para sistemas densos debido a las inversiones de matrices requeridas para la hidrodinámica de muchos cuerpos, mientras que los modelos de Método de Elementos Discretos (DEM) a menudo simplifican demasiado la geometría de la partícula (por ejemplo, cadenas de esferas) o descuidan las fuerzas hidrodinámicas necesarias. Existe la necesidad de un modelo basado en partículas que sea computacionalmente tratable y que capture fielmente la física microscópica dominante de las suspensiones de esferocilindros no Brownianas densas bajo cizalladura.

Metodología
Los autores desarrollaron un modelo de simulación basado en partículas para suspensiones densas de esferocilindros no flexibles y no Brownianos (cilindros con extremos hemisféricos) bajo flujo de cizalladura simple. El modelo opera en un régimen donde la inercia de la partícula está presente pero es despreciable.

Ecuaciones de Movimiento: Las trayectorias de las partículas se actualizan utilizando el algoritmo de Verlet-Velocidad, integrando las ecuaciones de movimiento traslacional y rotacional. Estas ecuaciones consideran fuerzas de contacto de pares, fuerzas de lubricación de pares y fuerzas hidrodinámicas de cuerpo único (arrastre, sustentación y torque).
Hidrodinámica: Las fuerzas de cuerpo único se modelan utilizando soluciones analíticas para esferoides prolatos como una aproximación para los esferocilindros. Esto incluye arrastre dependiente de la orientación y sustentación inducida por cizalladura. El torque hidrodinámico reproduce las órbitas de Jeffery.
Contacto y Lubricación:
- Contacto: Modelado como Hookeano con fricción de Coulomb, incluyendo componentes de resorte-amortiguador normales y tangenciales.
- Lubricación: Se desarrolló un esquema novedoso para calcular las fuerzas de lubricación entre esferocilindros. Para abordar la discontinuidad en la curvatura en la unión extremo-eje, el modelo calcula tres componentes de interacción (eje-eje, extremo-eje y extremo-extremo) simultáneamente. Estos se ponderan basándose en la posición relativa de las partículas para asegurar transiciones suaves entre los tipos de interacción, resolviendo las discontinuidades presentes en modelos anteriores (por ejemplo, Yamane et al., Butler y Shaqfeh).
Paso de Tiempo Dinámico: Una innovación metodológica clave es el cálculo de un paso de tiempo dinámico. El paso de tiempo se determina por el tiempo característico más pequeño entre las interacciones de contacto, las fuerzas de lubricación y los torques de lubricación. Esto asegura la estabilidad numérica mientras evita pasos de tiempo innecesariamente pequeños, mejorando la eficiencia computacional.
Configuración de la Simulación: El sistema utiliza fronteras periódicas en las direcciones de flujo y gradiente, y fronteras de Lees-Edwards para la cizalladura. Se utiliza una mezcla bidispersa de radios (1:1.4) para prevenir la cristalización. Las simulaciones cubren relaciones de aspecto ( $A$ ) de 1 a 20 y fracciones de volumen ( $\phi$ ) desde diluido hasta cerca del empaquetamiento máximo (jamming).

Resultados Clave

Comportamiento de la Viscosidad: El modelo predice un aumento sistemático de la viscosidad en estado estacionario ( $\eta_s$ $η_{s}$ ) con la fracción de volumen ( $\phi$ $ϕ$ ) y la relación de aspecto ( $A$ $A$ ).
- Respuesta Transitoria: Al iniciar la cizalladura, se observa un pico de viscosidad, atribuido a la distribución inicial aleatoria de las orientaciones de las partículas. Este pico disminuye a medida que las partículas se alinean, llevando al sistema a un estado estacionario. La magnitud del pico relativo al estado estacionario aumenta con $\phi$ .
- Estado Estacionario: La viscosidad en estado estacionario sigue una divergencia de ley de potencia $\eta_s \sim (1 - \phi/\phi_m)^{-\varsigma}$ . La fracción de volumen de empaquetamiento $\phi_m$ disminuye conforme $A$ aumenta (de $\approx 0.57$ para $A=5$ a valores más bajos para $A$ mayores). El modelo confirma la independencia de la tasa de la viscosidad en estado estacionario, consistente con el marco de trabajo para suspensiones granulares.
- Contribuciones de Fuerza: A $\phi$ bajo, la viscosidad está dominada por el esfuerzo del fluido y las fuerzas de lubricación. A medida que $\phi$ aumenta, las fuerzas de contacto (tanto normales como tangenciales) se vuelven dominantes, excediendo eventualmente las contribuciones de lubricación alrededor de $\phi \approx 0.5$ .
Microestructura y Alineación:
- Parámetro de Orden: El parámetro de orden orientacional ( $S$ ) aumenta con $\phi$ hasta un punto crítico, tras el cual disminuye ligeramente debido a la frustración geométrica que impide una alineación perfecta.
- Dirección de Alineación: Las partículas se alinean preferencialmente a lo largo de la dirección del flujo ( $x$ -axis), con un pequeño ángulo de inclinación ( $\phi_v$ ) respecto al plano de flujo-gradiente. Este ángulo de inclinación aumenta con $\phi$ hasta un umbral y luego disminuye.
- Efecto de la Relación de Aspecto: Relaciones de aspecto más altas conducen a una alineación más fuerte a fracciones de volumen más bajas, pero también resultan en un punto de empaquetamiento más bajo. Para $A$ muy alto, el ángulo de inclinación se vuelve insignificante.
- Transitorio vs. Estado Estacionario: El pico inicial de viscosidad se correlaciona con el tiempo requerido para que el sistema pase de un estado aleatorio a uno de cizalladura alineada.

Significancia y Reivindicaciones
El artículo afirma proporcionar una herramienta fundamental para explorar la mecánica de fluidos, la microestructura y la reología de las suspensiones granulares densas de varillas. Al modelar con éxito la transición de regímenes diluidos a densos y capturar la interacción entre fuerzas hidrodinámicas, contactos friccionales y la geometría de las partículas, el modelo corrobora los limitados datos experimentales de reología disponibles para varillas granulares.

Los autores sostienen que su trabajo extiende la comprensión de la reología de suspensiones granulares (previamente establecida para esferas) a la amplia clase de partículas con forma de varilla. El modelo ofrece una plataforma para estudios futuros para mapear la dependencia del punto de empaquetamiento sobre los detalles microestructurales, investigar esfuerzos normales y explorar el comportamiento bajo condiciones inhomogéneas. El artículo enmarca modestamente su contribución como un paso hacia una descripción completa de la mecánica de fluidos de estos materiales, dejando para el trabajo futuro la cuantificación exhaustiva de la relación entre el punto de empaquetamiento, la alineación y la fricción. No pretende resolver problemas específicos de procesamiento industrial directamente, sino que proporciona la infraestructura de simulación para investigar la física subyacente de tales sistemas.

Más como este