🔬 condensed matter

Granular thermostat implementation within the soft-sphere Discrete Element Method (DEM) framework, considerations and limitations

Este artículo aborda las limitaciones de los termostatos convencionales en sistemas granulares disipativos mediante la introducción de dos novedosas formulaciones híbridas por pares que controlan eficazmente la temperatura granular preservando al mismo tiempo las correlaciones dinámicas esenciales, permitiendo así el estudio de las propiedades reológicas dependientes de la temperatura en simulaciones del Método de Elementos Discretos.

Autores originales: Marco Previtali, Herbert Eric Huppert, Sergio Andres Galindo Torres

Publicado 2026-09-21

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Marco Previtali, Herbert Eric Huppert, Sergio Andres Galindo Torres

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

La arena, la nieve e incluso los posos de café comparten una extraña dualidad. Bajo una presión lenta y constante, se comportan como un sólido, entrelazándose en una red rígida de fricción. Pero si se agitan, se vierten o se agitan con la suficiente rapidez, de repente fluyen como un líquido. Los científicos que estudian estos materiales, conocidos como materia granular, dependen de simulaciones por computadora para comprender exactamente cómo ocurre este cambio. Rastrean millones de partículas diminutas e individuales mientras colisionan, rebotan y chocan entre sí. Un desafío clave en estas simulaciones es controlar la "temperatura" del sistema. En este contexto, la temperatura no significa calor en el sentido tradicional del fuego o el agua hirviendo; es una medida de qué tan salvajemente vibran y fluctúan las partículas en su movimiento. Si las partículas dejan de vibrar, la simulación se congela en un estado sólido. Si vibran demasiado, el material se dispersa. Para mantener una simulación en un estado específico y estable, los investigadores utilizan una herramienta digital llamada termostato. Al igual que un termostato doméstico enciende o apaga un calentador para mantener una habitación a una temperatura confortable, un termostato digital añade o elimina energía de las partículas para mantener su vibración en un nivel preciso.

Durante décadas, los científicos han utilizado dos tipos principales de termostatos digitales para gestionar estas simulaciones. Un tipo actúa como una resistencia, suave y constante, que frena a las partículas rápidas y acelera a las lentas para mantener un nivel de energía promedio. El otro tipo actúa como un bucle de retroalimentación, midiendo constantemente la energía del sistema y aplicando una fuerza correctiva para mantenerlo en el objetivo. Sin embargo, cuando los investigadores aplicaron estas herramientas estándar a pilas densas y agitadas de material granular, descubrieron que los métodos eran defectuosos. El enfoque basado en la resistencia era demasiado agresivo; suavizaba tanto el movimiento natural y caótico de las partículas que mataba la propia dinámica que los científicos querían estudiar. El enfoque basado en la retroalimentación, aunque era mejor para mantener la energía promedio estable, no lograba romper los patrones antinaturales que se forman cuando las partículas colisionan. En una pila de arena real, las colisiones son desordenadas y aleatorias, pero en estas simulaciones, las partículas comenzaban a agruparse y a moverse en ondas sincronizadas y antinaturales, creando un estado que era matemáticamente inestable y físicamente poco realista.

Para resolver esto, un equipo de investigadores de la Universidad de Westlake y la Universidad de Cambridge desarrolló un nuevo enfoque híbrido. Se dieron cuenta de que el problema era intentar usar una sola herramienta para realizar dos trabajos diferentes: controlar el nivel de energía general y romper los patrones antinaturales. Su solución fue combinar las fortalezas de los dos métodos antiguos eliminando sus debilidades. Crearon dos nuevos algoritmos que trabajan juntos. La primera parte del sistema actúa como una mano firme, asegurando que la cantidad total de energía de vibración se mantenga exactamente donde los investigadores desean. La segunda parte actúa como un suave empujón aleatorio, diseñado específicamente para interrumpir los patrones sincronizados y mantener las partículas moviéndose de una manera caótica y natural. Al separar estas dos tareas, los investigadores pudieron controlar la intensidad de la vibración sin congelar accidentalmente el sistema o forzar a las partículas a formar líneas ordenadas y falsas.

El equipo probó estos nuevos métodos ejecutando miles de simulaciones con diferentes configuraciones. Descubrieron que sus herramientas híbridas podían mantener una temperatura perfectamente estable incluso cuando las partículas colisionaban y perdían energía, una hazaña que los métodos antiguos tenían dificultades para lograr. Más importante aún, descubrieron que la forma en que aplicaban el empujón aleatorio cambiaba la naturaleza de la simulación de maneras sutiles pero significativas. Una versión de su nueva herramienta producía un sistema donde las partículas se movían con una aleatoriedad muy específica y predecible, mientras que la otra versión creaba un patrón estadístico ligeramente diferente. Esto demostró que la "temperatura" de un sistema granular no es solo un número único. Es un estado complejo que depende tanto de cuánta energía tienen las partículas como de cómo se distribuye esa energía entre ellas. Dos sistemas podrían tener exactamente la misma temperatura promedio pero comportarse de manera muy diferente dependiendo de cómo se inyecte la energía.

Para demostrar que su método funcionaba en un escenario realista, los investigadores lo aplicaron a un problema clásico: una caja de arena siendo cizallada, o deslizada, entre dos paredes. En esta configuración, podían mantener la presión constante mientras cambiaban la temperatura de las partículas que vibraban. Observaron que, a medida que aumentaban la temperatura, la fricción entre las partículas disminuía, haciendo que el material fluyera más fácilmente. Esto coincidía con lo que los científicos habían visto en experimentos previos, confirmando que su herramienta digital capturaba con precisión la física del material. Sin embargo, los resultados también mostraron que la relación entre la temperatura y la fricción no era tan simple como una sola fórmula pudiera describir; el modo específico en que se controlaba la temperatura importaba, lo que sugiere que para predecir completamente cómo se comportarán los materiales granulares, los científicos deben tener en cuenta algo más que la energía promedio; también deben comprender las estadísticas detalladas de cómo se mueven e interactúan las partículas.

Este trabajo proporciona una forma nueva y fiable de estudiar la física de la arena, la nieve y otros materiales granulares. Al proporcionar a los investigadores una herramienta que puede mantener la temperatura estable sin distorsionar el movimiento natural de las partículas, abre la puerta a modelos más precisos de deslizamientos de tierra, mezcla industrial y flujos geológicos. El estudio no pretende haber resuelto todos los misterios de la materia granular, pero ha eliminado un obstáculo importante. Demuestra que, al separar cuidadosamente el control de la energía del control del movimiento, los científicos pueden crear un entorno digital estable y realista para probar sus teorías. Esto les permite plantear preguntas más agudas sobre cómo fluyen y se atascan estos materiales, acercándose a una comprensión completa del complejo mundo de la física granular.

Resumen Técnico: Implementación de un Termostato Granular dentro del Método de Elementos Discretos (DEM) de Esferas Blandas

Planteamiento del Problema
Los flujos granulares densos se modelan frecuentemente utilizando el método de elementos discretos (DEM) de esferas blandas, aunque las aplicaciones a gran escala suelen depender de modelos de continuo que requieren relaciones constitutivas que vinculen el esfuerzo y la deformación con las variables de estado. La temperatura granular ( $T$ ), definida como la varianza de las fluctuaciones de velocidad no afines de las partículas, es una variable crítica para mejorar el escalamiento reológico inercial (por ejemplo, $\mu(I, \Theta)$ ). Sin embargo, las relaciones existentes basadas en la temperatura carecen de poder predictivo porque no describen la producción, el transporte y la disipación de $T$ . Además, los métodos de forzamiento experimental (por ejemplo, vibración) alteran la temperatura junto con las estadísticas de velocidad, las correlaciones y la anisotropía estructural, lo que dificulta aislar los efectos específicos de la temperatura.

Aunque los algoritmos de termostato son estándar en la Dinámica Molecular (MD) para controlar la energía cinética, su aplicación a sistemas DEM de esferas blandas disipativos está poco explorada. Los termostatos convencionales enfrentan limitaciones complementarias en los sistemas granulares:

Termostatos de Langevin: Para compensar la disipación colisional continua, requieren un acoplamiento fuerte, lo que sobreamortigua la dinámica de las partículas y suprime la difusión.
Termostatos de Nosé-Hoover: Aunque preservan mejor la dinámica, no logran interrumpir las correlaciones de velocidad y la segregación generadas por colisiones inelásticas repetidas. Esto conduce a estados agrupados (clustering) no ergódicos y a la inestabilidad numérica en el DEM de esferas blandas.

Metodología
Los autores investigan los algoritmos de termostato convencionales y proponen dos nuevas formulaciones híbridas por pares diseñadas para controlar la temperatura granular permitiendo la manipulación independiente de las estadísticas de velocidad. El estudio utiliza la biblioteca de código abierto Mechsys con aceleración por GPU.

Configuración del Sistema: Las simulaciones emplean una distribución de tamaño de partícula cuasi-monodispersa en un dominio cúbico periódico con una fracción de sólido $\phi \approx 0.27$ . Las interacciones siguen una ley de contacto elástica lineal (Hookeana) con un amortiguador viscoso y fricción de Coulomb.
Formulaciones Híbridas:
1. Adición Directa (NA): Combina un bucle de retroalimentación determinista de Nosé-Hoover (que controla la temperatura media) con un término de Langevin estándar (forzamiento estocástico y amortiguamiento).
2. Langevin Escalado (SL): Conserva la retroalimentación de Nosé-Hoover pero elimina el término de amortiguamiento de Langevin explícito. En su lugar, la amplitud del forzamiento estocástico se escala basándose en el estado cinético de par instantáneo ( $T_p$ ). El forzamiento desaparece si la temperatura de par excede el objetivo, imitando el comportamiento de transferencia de energía de los límites mecánicamente vibrantes.
Implementación: Ambas formulaciones utilizan un enfoque de pares (inspirado en la Dinámica de Partículas Disipativas) para asegurar la invarianza de Galileo y conservar el momento del centro de masa de los pares. Un rango de interacción ( $r_t$ ) define los pares del termostato, y una función de suavizado minimiza el ruido numérico en el límite del rango.
Estudios Paramétricos: Los autores realizaron tres estudios paramétricos variando:
1. El rango de interacción ( $R_t$ ) y las fuerzas de acoplamiento ( $\xi_{La}, \xi_{NH}$ ).
2. El tamaño del paso de tiempo ( $\Delta t$ ), la temperatura objetivo ( $\hat{T}$ ) y la fuerza de forzamiento estocástico adimensional ( $\pi_2 = \sqrt{\hat{T}}/(d\xi_{La})$ ).
3. El coeficiente de restitución ( $e$ ) para probar los efectos de la inelasticidad.
Aplicación: El marco se aplicó a una prueba de cizalladura simple controlada por presión para explorar el espacio de estados $(I, T)$ .

Resultos Clave

Limitaciones de los Termostatos Convencionales: Los termostatos de Langevin puros no lograron imponer la temperatura objetivo en sistemas disipativos sin un acoplamiento fuerte, lo que suprimió la difusión y llevó a la cristalización. Los termostatos de Nosé-Hoover puros mantuvieron la temperatura objetivo pero resultaron en estados altamente no ergódicos con grandes picos de temperatura, segregación de velocidad e inestabilidad numérica.
Desempeño de las Formulaciones Híbridas: Tanto NA como SL lograron imponer la temperatura escalar prescrita evitando las inestabilidades numéricas de Nosé-Hoover puro y el sobreamortiguamiento de Langevin puro.
Firmas Estadísticas Distintas: A pesar de controlar la misma temperatura escalar, las dos formulaciones híbridas produjeron estados cinéticos mensurablemente diferentes:
- NA produjo distribuciones de velocidad mesocúrticas (tipo Gaussiana) y mayor difusividad.
- SL produjo distribuciones platicúrticas (colas más planas) y menor difusividad.
  Esto demuestra que la magnitud de la fluctuación, la decorrelación y las estadísticas de forzamiento son aspectos separados del estado granular.
Rol del Forzamiento Estocástico: El parámetro adimensional $\pi_2$ gobierna la fuerza relativa de la decorrelación estocástica. Variar $\pi_2$ a una temperatura fija alteró las estadísticas de velocidad, la difusividad y el pico de la función de distribución radial $g(d)$ .
Aplicación de Cizalladura: En la cizalladura simple controlada por presión, el aumento de la temperatura granular a un número inercial ( $I$ ) fijo redujo la fricción aparente, lo cual es consistente con tendencias previas. Sin embargo, la relación no siguió estrictamente leyes de escalamiento simples, lo que sugiere que se requieren variables de estado adicionales.

Significancia y Reivindicaciones
El artículo afirma proporcionar un marco para generar estados de referencia controlados en simulaciones granulares. Al desacoplar la magnitud de las fluctuaciones de velocidad (temperatura) del carácter estadístico del forzamiento (mediante $\pi_2$ y la elección entre NA y SL), los investigadores pueden:

Identificar qué propiedades de transporte y reológicas están gobernadas principalmente por la temperatura granular escalar frente a aquellas que dependen de las estadísticas de velocidad y las correlaciones.
Comparar sistemas forzados de manera diferente (experimentales o naturales) sobre una base de estado (coincidiendo en $I$ y $T$ ), en lugar de una base de parámetros.
Determinar las variables específicas requeridas para cerrar los modelos constitutivos dependientes de la temperatura.

Los autores enfatizan que, si bien los termostatos híbridos ofrecen un método para aislar variables, las relaciones empíricas derivadas (por ejemplo, la dependencia de $\pi_2$ y $e$ ) son específicas al rango de parámetros investigado y requieren una validación adicional con variaciones en las fracciones de sólido y los coeficientes de fricción antes de servir como cierres constitutivos generales. El trabajo no pretende reproducir un aparato experimental específico, sino construir estados idealizados para el análisis reológico comparativo.

Más como este