🔬 physics

A Scalable OpenLB LAMMPS Framework for Fully Resolved Simulations of Hindered Settling of Arbitrary Non-Spherical Particles

Este artículo presenta un marco de trabajo escalable OpenLB-LAMMPS que acopla los métodos de Lattice Boltzmann y de Elementos Discretos para simular la sedimentación impedida totalmente resuelta de hasta 100.000 partículas no esféricas arbitrarias, revelando que las partículas cúbicas exhiben dinámicas de capa de coordinación y correlaciones de velocidad dependientes del tamaño del sistema distintas en comparación con las esferas.

Autores originales: Varghese Babu, Adrian Kummerländer, Mathias J. Krause, Santosh Ansumali, Dilip Gersappe

Publicado 2026-08-27

📖 4 min de lectura☕ Lectura para el café

CC BY 4.0

Autores originales: Varghese Babu, Adrian Kummerländer, Mathias J. Krause, Santosh Ansumali, Dilip Gersappe

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 cubo de arena vertiéndose en agua. Los granos no caen a la misma velocidad que si estuvieran solos; se frenan unos a otros, colisionando y empujando contra el fluido en una danza compleja y concurrida. Este fenómeno, conocido como sedimentación dificultada, es un proceso fundamental en la naturaleza y la industria, que rige desde cómo el sedimento construye los deltas de los ríos hasta cómo se mezclan los productos farmacéuticos en una fábrica. Durante décadas, los científicos han comprendido bien este comportamiento cuando las partículas son esferas perfectas, como pequeñas canicas. Sin embargo, el mundo real rara vez está hecho de esferas perfectas. La mayoría de las partículas son irregulares, dentadas o bloqueadas, y sus formas extrañas crean maneras únicas de tocarse y deslizarse unas sobre otras. Predecir cómo se asientan estas formas extrañas ha sido un desafío computacional obstinado porque las matemáticas de sus colisiones y la forma en que arrastran el fluido a su alrededor son increíblemente difíciles de resolver simultáneamente.

Para abordar esto, un equipo de investigadores ha construido un nuevo marco digital que acopla dos potentes herramientas de software de código abierto para simular estas interacciones complejas. Combinaron un resolvedor de fluidos llamado OpenLB, que modela cómo se mueven los líquidos, con un simulador de partículas llamado LAMMPS, que rastrea cómo los objetos sólidos chocan y rozan entre sí. La clave de la innovación reside en cómo representan los objetos sólidos. En lugar de intentar modelar una forma compleja como un cubo con una única definición matemática rígida, descomponen el objeto en un grupo de esferas superpuestas, de forma muy parecida a un muñeco de nieve hecho de tres bolas. Este "grupo" permite que la computadora calcule fácilmente cómo el objeto toca a sus vecinos. Al mismo tiempo, el resolvedor de fluidos no ve el objeto como una colección de bolas, sino como un bloque sólido con una forma específica, utilizando una técnica llamada voxelización para mapear la superficie del objeto en la rejilla del fluido. Esta separación permite al equipo simular la mecánica de contacto de las partículas y el flujo del fluido con alta precisión, incluso cuando se trata de cantidades masivas de objetos.

Los investigadores probaron su sistema simulando el asentamiento de cubos en un fluido, comparando los resultados con comportamientos conocidos de las esferas. Realizaron simulaciones con hasta 103.823 cubos, una escala que lleva al límite la capacidad de computación actual. Sus hallazgos revelaron una diferencia distintiva entre cómo se comportan los cubos y las esferas mientras se asientan. Mientras que las esferas tienden a agruparse mediante el agrupamiento dominado por el contacto, los cubos formaron capas de coordinación pronunciadas sin contacto paralelo a las caras. Esto sugiere que los bordes afilados y las esquinas de los cubos crean un tipo específico de espaciamiento que evita que se empaqueten tan estrechamente como las esferas. Además, el equipo descubrió que las fluctuaciones de velocidad de las partículas no eran aleatorias ni locales; en cambio, el movimiento de una partícula influía en el movimiento de otras en todo el sistema. Incluso en las simulaciones más grandes, la correlación entre las velocidades de las partículas se extendía por todo el contenedor, indicando que el movimiento colectivo de la multitud es un fenómeno unificado en lugar de una serie de eventos aislados.

Al validar su modelo mediante experimentos con esferas y luego aplicarlo a cubos, los autores demostraron que la forma de las partículas juega un papel crítico en la dinámica de suspensión. Las simulaciones mostraron que la velocidad de asentamiento de los cubos difiere de la de las esferas, y que la estructura interna de la nube de asentamiento es fundamentalmente distinta. El trabajo confirma que ignorar la forma de las partículas puede conducir a predicciones inexactas en campos que van desde la ingeniería geotécnica hasta la ciencia ambiental. El marco del equipo proporciona una herramienta robusta para estudiar estos flujos a gran escala, ofreciendo una forma de ver el orden oculto dentro del caos del asentamiento de las partículas. Aunque las simulaciones actuales no incluyeron la fina película de fluido que lubrica el espacio entre las partículas cuando están extremadamente cerca, el método capturó con éxito la dinámica general del sistema. Los resultados sugieren que, para las partículas no esféricas, la forma en que se asientan está gobernada por una compleja interacción de forma, contacto e interacciones de fluido de largo alcance que los modelos simples de esferas no pueden capturar.

Resumen Técnico: Un marco escalable de OpenLB–LAMMPS para simulaciones totalmente resueltas de sedimentación dificultada de partículas no esféricas arbitrarias

Planteamiento del problema
La sedimentación dificultada, el fenómeno en el que las interacciones partícula-partícula y partícula-fluido retrasan la velocidad descendente de suspensiones de alta concentración, es significativamente menos comprendida para partículas no esféricas que para esferas. Mientras que la relación de Richardson-Zaki describe eficazmente las suspensiones esféricas, los modelos numéricos para formas no esféricas enfrentan desafíos computacionales sustanciales. Estas dificultades surgen principalmente de la necesidad de resolver simultáneamente la mecánica de contacto compleja (trabazón, dinámica dependiente de la orientación) y las interacciones hidrodinámicas. Los marcos existentes de Dinámica de Fluidos Computacional–Método de Elementos Discretos (CFD-DEM) suelen tener dificultades para simular sistemas a gran escala ( $\mathcal{O}(10^5)$ partículas) con formas de cuerpos rígidos arbitrarios debido a las limitaciones en la detección de contactos y la resolución de la interacción fluido-estructura (FSI).

Metodología
Los autores presentan un marco de FSI escalable y totalmente resuelto que acopla dos paquetes de software de código abierto: OpenLB (un solver del Método de Lattice Boltzmann) y LAMMPS (un solver del Método de Elementos Discretos).

Solver de fluido (OpenLB): La fase fluida se modela mediante el Método de Lattice Boltzmann (LBM) con un enfoque homogeneizado basado en las ecuaciones de Brinkman-Navier-Stokes (BNSE). Los cuerpos sólidos arbitrarios se representan en una red Euleriana fija mediante un campo de porosidad dependiente del tiempo, $d(\mathbf{x}, t)$ . La condición de no deslizamiento se impone a través de un término de penalización de volumen (arrastre de Brinkman). Las fuerzas y torques hidrodinámicos se calculan mediante el intercambio de momento a través de las interfaces fluido-sólido.
Solver de partículas (LAMMPS): Los cuerpos rígidos se modelan utilizando el Método de Elementos Discretos (DEM). Para manejar formas arbitrarias, los autores emplean una representación de "clump" (aglomerado), donde geometrías complejas (por ejemplo, cubos) se descomponen en esferas superpuestas. Esta representación se utiliza estrictamente para la detección de contactos y el cálculo de fuerzas dentro de LAMMPS.
Estrategia de acoplamiento: El marco opera bajo un esquema de paso de tiempo secuencial. El solver de fluido calcula las fuerzas hidrodinámicas basadas en las posiciones y orientaciones actuales de las partículas, las cuales se aplican luego como fuerzas externas en la simulación DEM. Inversamente, los estados actualizados de las partículas (posiciones, velocidades, velocidades angulares) se mapean de vuelta a la rejilla de fluido como un campo de porosidad.
- Crucialmente, el solver de fluido resuelve la geometría exacta de la superficie de las partículas mediante voxelización, independientemente de la aproximación de "clump" de múltiples esferas utilizada para los contactos en LAMMPS. Esto distingue este trabajo de métodos previos donde el fluido también "ve" la aproximación de clump, evitando así errores geométricos adicionales en la interacción hidrodinámica.
Configuración de la simulación: El estudio se centra en la sedimentación dificultada en un dominio triplemente periódico. Para evitar la aceleración infinita debida a la gravedad, se aplica una fuerza compensatoria a los nodos del fluido para mantener una velocidad media de fluido de cero.

Contribuciones clave

Marco escalable: Los autores demuestran la capacidad de simular sistemas que contienen hasta $\mathcal{O}(10^5)$ cuerpos rígidos arbitrarios (específicamente cubos), un orden de magnitud mayor que las demostraciones previas de código abierto (por ejemplo, openHFDIB-DEM).
Representación híbrida: La integración de la voxelización de superficie exacta para la resolución del fluido con los clumps de múltiples esferas para la mecánica de contacto permite la simulación eficiente de geometrías complejas, como cubos, sin sacrificar la resolución hidrodinámica.
Validación: El marco se valida contra referencias experimentales y analíticas, incluyendo la sedimentación de una sola esfera (Ten Cate et al.), colisiones esfera-pared, arrastre en arreglos periódicos de esferas y la velocidad de sedimentación de partículas cuboidales (Kwon et al.).

Resultados

Dinámica de sedimentación: Las simulaciones de sedimentación dificultada para cubos y esferas equivalentes en volumen revelan diferencias distintivas. Mientras que las esferas tienden a formar cúmulos dominados por el contacto, los cubos forman capas de coordinación pronunciadas sin contacto paralelo a las caras, lo que indica que la morfología de la partícula altera significativamente la microestructura de la suspensión.
Funciones de Distribución Radial (RDF): La RDF para los cubos muestra capas primera y segunda bien definidas, mientras que las esferas muestran una tendencia a tocarse directamente, y el número de contactos aumenta con la fracción de empaquetamiento.
Correlaciones de velocidad: El estudio investiga la correlación espacial de las fluctuaciones de velocidad. Un hallazgo clave es que la longitud de correlación de las fluctuaciones de velocidad escala con el tamaño del sistema, incluso para los sistemas más grandes estudiados ( $\sim 10^5$ partículas). A diferencia de algunas observaciones experimentales donde las fluctuaciones saturan en una longitud de escala específica, las simulaciones muestran correlaciones que se extienden a través de todo el dominio.
Dependencia del tamaño del sistema: Aunque la velocidad terminal promedio y la disposición local de las partículas permanecen estables a través de diferentes tamaños de sistema (desde $\sim 3,000$ hasta $\sim 100,000$ cubos), la naturaleza de largo alcance de las correlaciones de velocidad persiste, lo que sugiere que las condiciones de contorno periódicas o el mecanismo específico de fuerza de equilibrio pueden influir en las correlaciones observadas.

Significancia y Reivindicaciones
El artículo afirma proporcionar un marco robusto y de código abierto para simular flujos de partículas y geotécnicos a gran escala que involucren formas no esféricas arbitrarias. Al acoplar exitosamente OpenLB y LAMMPS, los autores cierran la brecha entre la mecánica de contacto precisa y la hidrodinámica totalmente resuelta para partículas no esféricas. Sus hallazgos resaltan el papel crítico de la morfología de la partícula en la dinámica de la suspensión, mostrando específicamente que los cubos exhiben comportamientos de agrupamiento diferentes en comparación con las esferas, así como correlaciones de fluctuación de velocidad distintas. Los autores señalan modestamente que, si bien el marco es escalable, se requiere trabajo futuro para abordar las correlaciones de velocidad de largo alcance observadas (potencialmente artefactos de las condiciones de contorno) e incorporar fuerzas de lubricación para interacciones de partículas más cercanas. El trabajo establece una base para el estudio de sistemas de partículas complejos a escalas previamente inaccesibles para los métodos CFD-DEM totalmente resueltos.

Más como este