Percolation of a cohesive fine particle in a static bed

Resumen Técnico: Percolación de una Partícula Fina Cohesiva en un Lecho Estático

Enunciado del Problema
La percolación de partículas finas a través de un lecho estático de partículas más grandes es un proceso fundamental en el manejo industrial (por ejemplo, en farmacéutica y procesamiento químico) y en flujos geofísicos (por ejemplo, transporte de sedimentos y contaminación de aguas subterráneas). En sistemas no cohesivos, si una partícula fina pasa o queda atrapada está gobernado principalmente por restricciones geométricas, específicamente la relación entre los diámetros de las partículas grandes y finas ($R = d_l/d_f$). Existe un umbral de atrapamiento geométrico en $R_t \approx 6.464$; cuando $R > R_t$, se espera teóricamente que las finas pasen libremente a través de los cuellos de poro formados por las partículas grandes. Sin embargo, los sistemas del mundo real a menudo involucran fuerzas cohesivas (fuerzas de van der Waals, electrostáticas, puentes líquidos) que alteran la dinámica de contacto. Aunque estudios anteriores han abordado la percolación no cohesiva y la reología de flujos cohesivos, existe una brecha en la comprensión de cómo la cohesión influye en la percolación en el régimen de "tamizado libre", donde las restricciones geométricas por sí solas permitirían el paso. Este estudio investiga la percolación de partículas finas cohesivas a través de un lecho estático de partículas grandes donde $R=7$ (por encima del umbral geométrico), centrándose en el límite de partícula única donde las finas interactúan con las partículas del lecho pero no entre sí.

Metodología
Los autores emplean simulaciones del Método de Elementos Discretos (DEM) utilizando el software de código abierto LAMMPS. Las simulaciones utilizan un esquema mejorado de Velocity-Verlet adecuado para grandes relaciones de tamaño.

• Configuración del Sistema: Un lecho estático de partículas grandes monodispersas ($d_l = 4$ mm) se empaqueta con una fracción $\phi = 0.6$. Las partículas finas ($d_f \approx 0.57$ mm, $R=7$) se liberan desde el reposo sobre el lecho y se dejan caer bajo la gravedad.
• Modelos de Contacto: Las fuerzas de contacto normales siguen el modelo elástico de Hertz. La inelasticidad se modela mediante un coeficiente de restitución constante ($e_n$). Las fuerzas tangenciales utilizan el modelo de Mindlin con fricción por deslizamiento ($\mu$) y resistencia a la rodadura.
• Cohesión: Las interacciones cohesivas se modelan principalmente utilizando el modelo Derjaguin-Muller-Toporov (DMT), que asume una fuerza atractiva constante durante el contacto y repulsión de Hertz. La fuerza de la cohesión se cuantifica mediante el número de Bond ($Bo = F_c / W$), definido como la relación entre la fuerza de despegue y el peso de la partícula fina. Se realizó también un número limitado de simulaciones utilizando el modelo Johnson-Kendall-Roberts (JKR) para comparación.
• Parámetros: El estudio varía sistemáticamente el número de Bond ($Bo$), el coeficiente de fricción por deslizamiento ($\mu$), el coeficiente de restitución ($e_n$) y el módulo de Young ($E$). Se rastrean las dinámicas de $10^4$ partículas finas independientes para garantizar la precisión estadística.

Resultados Clave

1. Atrapamiento No Geométrico: Incluso cuando $R=7$ (permitiendo geométricamente el paso libre), una cohesión y fricción suficientemente grandes conducen al atrapamiento de partículas. A medida que aumenta $Bo$, la distancia característica de percolación disminuye, lo que significa que las partículas quedan atrapadas a profundidades más someras.
2. Distribución de Atrapamiento: La distribución vertical de las partículas atrapadas sigue un decaimiento exponencial. La distancia característica de atrapamiento ($\lambda$) sirve como métrica para el comportamiento de percolación. La distancia inversa de atrapamiento ($d_l/\lambda$) aumenta monótonamente con $Bo$.
3. Rol de la Fricción y la Cohesión:
   • A bajo $Bo$, la fricción por deslizamiento juega un papel significativo en el atrapamiento.
   • A alto $Bo$, la cohesión se convierte en el mecanismo de atrapamiento dominante.
   • Se observa un "codo" en la relación entre la distancia inversa de atrapamiento y la fricción. Este codo se desplaza con $Bo$ y corresponde a una transición entre regímenes donde domina la fricción y donde domina la cohesión.
   • La rigidez de la partícula ($E$) influye en el atrapamiento: las partículas más blandas (menor $E$) se atrapan más fácilmente, particularmente a valores más altos de $Bo$.
4. Mecanismo de Colisión: El proceso de atrapamiento se descompone en dos pasos secuenciales:
   • Condición de Adhesión: Tras la colisión, una partícula rebota o se adhiere. La adhesión ocurre si la velocidad de impacto está por debajo de una velocidad crítica ($v_{cr}$), que depende de $Bo$, $e_n$ y $E$.
   • Movimiento Post-Contacto: Si la partícula se adhiere, puede deslizarse y desprenderse o permanecer atrapada. Una condición de equilibrio estático basada en el balance de fuerzas determina si una partícula que no rebota permanece estática. Esta condición depende de un coeficiente de fricción reescalado $\tilde{\mu} = \mu\sqrt{4Bo^2 - 1}$. Si $\tilde{\mu} > 1$, es probable que la partícula permanezca estática; si $\tilde{\mu} < 1$, es probable que se deslice y potencialmente se desprenda.
5. Modelo Predictivo: Los autores desarrollan un modelo de colisión que vincula la dinámica microscópica con la distancia macroscópica de atrapamiento. La distancia inversa de atrapamiento se predice mediante $1/\lambda = -k \ln(1 - P_t)$, donde $k$ es la frecuencia de colisiones y $P_t$ es la probabilidad de atrapamiento por colisión (determinada por la relación entre la velocidad de impacto y $v_{cr}$). Este modelo predice con precisión las distancias de atrapamiento en el régimen dominado por colisiones (alto $Bo$).
6. Doble Contacto: En regímenes con menor fricción o cohesión, las partículas a menudo se deslizan y forman dobles contactos (tocando dos partículas grandes simultáneamente). Esta configuración puede conducir al atrapamiento incluso cuando las condiciones de contacto único sugerirían el desprendimiento, o inversamente, permitir el desprendimiento mediante un rebote desde un segundo contacto.

Significado y Afirmaciones
El artículo establece que la percolación de partículas finas cohesivas no está determinada únicamente por la accesibilidad geométrica. Por el contrario, las dinámicas de interacción a escala de partícula —específicamente la interacción entre cohesión, restitución, fricción y velocidad de impacto— pueden anular las expectativas geométricas, conduciendo al atrapamiento incluso cuando los cuellos de poro son geométricamente lo suficientemente anchos para el paso.

Los autores afirman proporcionar un marco fundamental para comprender estos mecanismos en lechos estáticos, lo cual es un prerrequisito para abordar el problema más complejo de la segregación granular en mezclas en flujo de partículas cohesivas. Al identificar las condiciones específicas (velocidad crítica, criterios de equilibrio estático) que gobiernan la adhesión y el movimiento posterior, el estudio ofrece una capacidad predictiva para la distancia de atrapamiento basada en propiedades materiales fundamentales. Los resultados destacan que en sistemas cohesivos, el "destino" de una partícula fina está gobernado por una secuencia de interacciones locales (colisión, adhesión, movimiento post-contacto) y no solo por la geometría estática del lecho.