🔬 condensed matter

A Smoluchowski equation for a sheared suspension of frictionally interacting rods

Este artículo desarrolla una ecuación de Smoluchowski y un tensor de esfuerzos para suspensiones densas de varillas sometidas a cizallamiento mediante la aplicación del método variacional de Onsager de Doi para incorporar tanto la fricción sólida como la lubricada, al tiempo que propone un modelo generalizado para el número promedio de contactos entre varillas.

Autores originales: Genevieve Quiñones, Peter D. Olmsted

Publicado 2026-07-20

📖 8 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Genevieve Quiñones, Peter D. Olmsted

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 bebemos, pintamos o imprimimos no es solo un líquido simple, sino una danza caótica de diminutos palos rígidos flotando en un fluido. Este es el reino de las suspensiones, una rama de la física que estudia cómo se mueven las partículas cuando se mezclan en un líquido. Durante más de un siglo, los científicos han comprendido cómo se comportan estos palos cuando están alejados, derivando perezosamente como hojas en una corriente suave. También saben cómo actúan cuando están amontonados pero aún se deslizan unos sobre otros sobre una fina capa de fluido resbaladizo, como patinadores sobre hielo en un estanque congelado. Pero hay un tercer escenario, más desordenado: ¿qué sucede cuando la multitud se vuelve tan densa que el "hielo" se rompe? Cuando los palos chocan entre sí, pierden su recubrimiento resbaladizo y comienzan a triturarse unos con otros con fricción? Esta es la pregunta que impulsa la reología moderna (el estudio del flujo). Comprender esto es crucial porque explica por qué algunos materiales se convierten repentinamente de un líquido fluido en un bloque sólido cuando se los agita demasiado rápido, un fenómeno que puede obstruir tuberías industriales o atascar impresoras 3D.

En este artículo, Genevieve Quiñones y Peter D. Olmsted abordan la matemática desordenada de estos choques friccionales. Construyen un nuevo conjunto de reglas, llamado ecuación de Smoluchowski, para predecir cómo se comporta una multitud densa de varas friccionales cuando son sometidas a cizallamiento (agitación). Piensa en ello como una actualización de una simulación de tráfico, pasando de un modelo donde los coches solo conducen por carreteras suaves y sin fricción a uno donde realmente pueden derrapar, triturarse y bloquearse. Los autores sugieren que, cuando estas varas se amontonan, no solo se deslizan; a veces se pegan y se raspan, creando un "torque friccional" que cambia la forma en que rotan. Su trabajo sugiere que esta fricción actúa como un freno, suprimiendo los movimientos erráticos de vueltas que las varas suelen realizar y obligándolas a alinearse más ordenadamente con el flujo. También proponen una nueva forma de contar cuántas veces una vara choca con sus vecinas, encontrando que cuanto más ordenada es la multitud, menos choques ocurren. Aunque aún no han construido una máquina física para probar esto, sus simulaciones por computadora muestran que añadir fricción cambia los pasos de baile de las varas, evitando que giren fuera de control y explicando potencialmente por qué algunos fluidos espesos se vuelven aún más espesos cuando se les presiona más fuerte.

La historia de los palos que trituran

Imagina que estás en una pista de baile abarrotada. Si la sala es espaciosa, todos pueden girar y dar vueltas libremente. Esto es como una suspensión diluida de varas, donde las partículas están alejadas y solo chocan entre sí ocasionalmente. Los físicos han sabido durante mucho tiempo cómo predecir este baile; utilizan ecuaciones (como la ecuación de Jeffery) que describen cómo una sola vara rota en un líquido en movimiento.

Ahora, imagina que la pista de baile se llena. Los bailarines ya no solo están girando; están chocando entre sí. En una multitud semi-diluida, es posible que todavía se deslicen unos sobre otros debido a una fina capa de fluido (como sudor o aceite) que actúa como lubricante. Pero en una multitud concentrada, esa capa de lubricante puede romperse. De repente, los bailarines no se están deslizando; están triturando sus zapatos contra los de los demás. Esto es fricción sólida.

Durante mucho tiempo, los científicos tuvieron dificultades para escribir la matemática de este triturado. Los intentos anteriores asumían que, incluso cuando las varas estaban triturándose, seguían los mismos "pasos de baile" (ecuación de Jeffery) como si estuvieran en una multitud dispersa, simplemente con un "empujón" adicional añadido. Pero cuando los autores de este artículo intentaron calcular el efecto promedio de estos empujones, encontraron que la matemática se cancelaba a cero. Era como si la fricción fuera invisible. Esto no tenía sentido para ellos, porque sabemos por la vida real que la fricción ralentiza las cosas y cambia cómo se mueven.

El nuevo enfoque: La "factura de energía"

Quiñones y Olmsted decidieron probar una estrategia diferente. En lugar de adivinar cómo se mueven las varas y luego añadir la fricción, comenzaron con la energía. Utilizaron un método llamado Rayleighiano, que es esencialmente una forma de calcular la "factura de energía" de todo el sistema.

Piensa en el Rayleighiano como un calculador de costos totales. Suma dos cosas:

La Energía Libre: Cuánto quieren las varas organizarse de cierta manera (como imanes alineándose).
La Disipación: Cuánta energía se pierde en forma de calor debido a la fricción y la resistencia.

Los autores se dieron cuenta de que, para encontrar el movimiento real de las varas, hay que minimizar esta "factura de energía" total. Cuando incluyeron la energía perdida por la fricción sólida (el triturado) y la fricción lubricada por frontera (el deslizamiento con una fina película) en esta factura, algo mágico sucedió. La matemática produjo naturalmente un nuevo "torque" (una fuerza de torsión) que no era cero.

Descubrieron que la fricción actúa como un freno rotacional. Cuando las varas se trituran entre sí, se vuelve más difícil para ellas girar. Esto cambia la ecuación de Smoluchowski, que es la ecuación maestra que predice cómo cambia la orientación de todas las varas a lo largo del tiempo.

Lo que encontraron: El baile cambia

Los autores realizaron simulaciones por computadora para ver qué sucede cuando se agitan estas varas friccionales. Esto es lo que descubrieron:

El "kayaking" se detiene: En una multitud sin fricción, las varas a menudo realizan un movimiento llamado "kayaking", donde se inclinan y giran en una órbita compleja fuera del plano. Los autores descubrieron que, cuando se activa la fricción, este giro salvaje se suprime. Las varas dejan de hacer kayaking y, en cambio, tienden a alinearse con el flujo, como un banco de peces nadando en línea recta.
La conexión con el "atascamiento": Sugieren que este frenado friccional es un ingrediente clave en el Espesamiento por Cizallamiento Discontinuo (DST). Este es el extraño fenómeno por el cual un líquido (como el almidón de maíz y el agua) de repente se convierte en un sólido cuando se le golpea con fuerza. El modelo de los autores sugiere que, a medida que se agita más rápido, la fricción entre las varas aumenta, los "frenos" se bloquean y el material se atasca.
Contando los choques: También idearon una nueva fórmula para estimar cuántas veces una vara toca a sus vecinas. Descubrieron que el número de contactos depende de qué tan ordenada sea la multitud. Si las varas están perfectamente alineadas (alto orden), tocan a menos vecinos. Si son desordenadas y caóticas, tocan a más. Esto es una generalización de una idea previa que solo funcionaba para multitudes desordenadas.

Los límites del modelo

Es importante señalar lo que este artículo no afirma. Los autores tienen cuidado de decir que su modelo asume que las varas siempre se están deslizando unas sobre otras. Todavía no pueden modelar el momento en que las varas se quedan tan trabadas que forman un grupo gigante e inmóvil (fricción estática). En el mundo real, estos grupos trabados son probablemente la razón principal por la que los materiales se atascan por completo. Los autores admiten que su modelo omite esta parte de "estar trabado" porque la fricción estática no genera calor (disipación), por lo que no aparece en su factura de energía.

Además, sus resultados se basan en simulaciones y derivaciones matemáticas, no en nuevos experimentos físicos. Sugieren que sus hallazgos explican por qué la fricción importa, pero aún no la han medido en un laboratorio. También asumen que la suspensión está perfectamente mezclada y es uniforme, ignorando la posibilidad de que la fricción pueda causar que las varas se agrupen en puntos específicos.

La conclusión

En términos simples, Quiñones y Olmsted han construido una nueva caja de herramientas matemáticas que finalmente tiene en cuenta el "triturado" entre las partículas en un fluido espeso. Demostraron que si se ignora la fricción, se pierden los frenos que detienen el giro salvaje de las varas. Al incluir esta fricción, encontraron una forma autosuficiente de predecir cómo se mueven estas varas amontonadas y cuánta tensión ejercen sobre el recipiente que las contiene. Aunque no han resuelto el misterio de por qué los fluidos se atascan por completo (eso requiere comprender el estado de "trabado"), han proporcionado una pieza crucial del rompecabezas: el torque friccional que ralentiza el baile antes de que ocurra el atasco.

Resumen Técnico: Una ecuación de Smoluchowski para una suspensión sometida a cizalla de varillas con interacción friccional

Planteamiento del Problema
La reología de suspensiones densas de partículas rígidas en forma de varilla bajo cizalla no es capturada completamente por las teorías existentes, particularmente en regímenes que exhiben engrosamiento discontinuo de la cizalla (DST), bloqueo por cizalla (shear jamming) y rotación de la vorticidad (vorticity tilting). Mientras que los modelos clásicos (p. ej., Jeffery, Doi, Advani-Tucker) describen con éxito la dinámica orientacional en regímenes diluidos y semidiluidos mediante interacciones hidrodinámicas y efectos de volumen excluido, fallan al no poder contabilizar los contactos sólidos entre partículas. Evidencia experimental y computacional reciente sugiere que los contactos friccionales entre las varillas son el motor principal de estos comportos reológicos no explicados. Sin embargo, los intentos previos de incorporar la fricción en la dinámica de las varillas (p. ej., Sandstrom & Tucker, Bounoua et al.) arrojaron resultados nulos para el torque friccional promedio, debido principalmente a que asumieron que el movimiento de las varillas seguía la ecuación de Jeffery y simplemente añadieron un término de torque, lo que provocó cancelaciones al promediar.

Metodología
Los autores desarrollan un marco constitutivo para suspensiones densas de varillas con interacción friccional aplicando el método variacional de Onsager tal como fue formulado por Doi. El núcleo de la metodología consiste en construir un Rayleighano ( $R$ ), definido como la suma de la tasa de cambio de la energía libre ( $\dot{A}$ ) y una función de disipación ( $\Phi$ ) dependiente de la velocidad angular ( $\omega$ ) de la varilla.

Funciones de Disipación: Los autores derivan funciones de disipación para dos mecanismos friccionales distintos:
- Fricción de Lubricación de Frontera (BL): Escala linealmente con la velocidad de contacto tangencial relativa ( $v_c$ ), modelada como $F \propto v_c$ .
- Fricción de Amontons-Coulomb Sólida: Escala con la magnitud de la fuerza normal ( $|F_N|$ ) y depende únicamente de la dirección de $v_c$ , modelada como $F \propto |F_N|$ .
Enfoque de Campo Medio: Para manejar la naturaleza de muchos cuerpos de los contactos, los autores emplean una aproximación de campo medio. Derivan una expresión para el número promedio de contactos por varilla, $\langle c \rangle$ , como una función de la fracción de volumen ( $\phi$ ), la relación de aspecto ( $L/D$ ) y el parámetro de orden nemático ( $S$ ). Esto generaliza la ecuación de contactos aleatorios de Philipse para empaquetamientos desordenados para incluir el orden orientacional.
Derivación Variacional: A diferencia de otros enfoques que asumían una ecuación de movimiento predeterminada, los autores minimizan el Rayleighano total con respecto a $\omega$ . Esto produce una ecuación de movimiento autoconsistente donde el torque friccional es naturalmente una función de la velocidad angular misma.
Ecuación de Smoluchowski y Esfuerzo: La minimización del Rayleighano produce una nueva ecuación de Smoluchowski para la función de distribución orientacional $\psi(\hat{u})$ y un tensor de esfuerzo $\sigma$ correspondiente. El modelo incorpora una función de transición $\Theta(p)$ para transicionar entre los regímenes de fricción lubricada y sólida basándose en la presión reducida.

Contribuciones Clave

Torque Friccional No Nulo: El artículo demuestra que, al derivar la dinámica a partir de una función de disipación en lugar de añadir un torque a la ecuación de Jeffery, se obtiene un torque friccional promedio no nulo. Esto resuelve el "resultado nulo" encontrado en la literatura previa (p. ej., Bounoua et al.), el cual surgió al promediar potencias impares de las variables de posición de contacto.
Número de Contacto Generalizado: Los autores proponen una forma para el número promedio de contactos $\langle c \rangle$ que depende del parámetro de orden nemático $S$ , específicamente $\langle c \rangle \approx \phi (L/D) (1 - S^2)$ . Esto vincula la microestructura de los contactos directamente con el orden orientacional macroscópico.
Dinámica Autoconsistente: La ecuación de Smoluchowski resultante es autoconsistente; el coeficiente de arrastre rotacional efectivo depende de la velocidad de contacto, la cual a su vez depende de la velocidad angular. Esta no linealidad es inherente al modelo de fricción de Coulomb.
Formulación del Tensor de Esfuerzo: Se deriva un nuevo tensor de esfuerzo que incluye contribuciones tanto del arrastre hidrodinámico como de la disipación friccional, donde la componente friccional escala con el torque friccional promedio.

Resultados y Estudio Numérico
Los autores realizan un estudio numérico de la ecuación de Smoluchoso derivada (aproximada mediante una expansión de serie truncada y dinámica de Langevin) para un flujo de cizalla simple.

Límite Lubricado: En ausencia de fricción sólida, el modelo se reduce a la teoría de Doi-Edwards con una constante de arrastre rotacional renormalizada, lo que ralentiza efectivamente la rotación de las varillas.
Límite de Fricción Sólida: Cuando la fricción sólida está activa ( $\mu_k = 1$ $μ_{k} = 1$ ), la dinámica cambia significativamente. El modelo predice una supresión del movimiento del director en el plano flujo-gradiente.
- Alineación con el Flujo: A diferencia del caso sin fricción donde el director puede exhibir "kayaking" (balanceo fuera del plano) u órbitas periódicas a números de Peclet intermedios, la inclusión de la fricción estabiliza al director en un estado alineado con el flujo.
- Condiciones Iniciales: El modelo friccional muestra que incluso partiendo de una distribución isotrópica, el sistema tiende a alinearse con el flujo y permanecer allí, con solo excursiones esporádicas hacia órbitas de tipo kayaking.
- Inclinación de la Vorticidad: A números de Peclet bajos, la componente de la vorticidad del director deriva cerca de cero, contrastando con el vagabundeo errático visto en el caso sin fricción.

Significancia
El artículo afirma proporcionar una descripción teórica de las suspensiones densas de varillas que contabiliza explícitamente los contactos friccionales, un factor que ha sido previamente ignorado o tratado incorrectamente en los modelos constitutivos. Al utilizar el principio variacional de Onsager, los autores establecen un marco autoconsistente donde la disipación friccional genera naturalmente los torques necesarios para alterar la dinámica orientacional. El trabajo sugiere que los contactos friccionales son un mecanismo crítico para la supresión de movimientos orientacionales complejos (como el kayaking) y pueden ser responsables del inicio del engrosamiento por cizalla y el bloqueo en suspensiones de varillas densas. Los autores señalan que, aunque su modelo asume homogeneidad espacial y fricción por deslizamiento (excluyendo cúmulos de adherencia estática), ofrece una base rigurosa para comprender las anomalías reológicas observadas en suspensiones de varillas concentradas.

La historia de los palos que trituran

El nuevo enfoque: La "factura de energía"

Lo que encontraron: El baile cambia

Los límites del modelo

La conclusión

Más como este