Densely packed particle raft at vertically vibrated air-water interface

Resumen Técnico: Balsa de partículas densamente empaquetadas en una interfaz aire-agua vibrada verticalmente

Planteamiento del Problema
Las balsas de partículas que se forman en interfaces fluidas exhiben comportamientos mecánicos ricos, que van desde la elasticidad de tipo sólido hasta el flujo de tipo líquido. Si bien la respuesta de estas balsas a deformaciones cuasiestáticas (deformación de tracción, cizallamiento, indentación) está bien documentada, su dinámica bajo excitaciones dinámicas fuertes y homogéneas sigue siendo poco comprendida. Específicamente, la interacción entre las fuerzas capilares, las fuerzas hidrodinámicas y el empaquetamiento de las partículas es compleja en un entorno de vibración vertical. Estudios previos se han centrado en fracciones de empaquetamiento bajas o en la generación de ondas localizadas, dejando un vacío en la comprensión de cómo las balsas granulares densas responden a oscilaciones verticales globales, particularmente en lo que respecta a la emergencia de patrones colectivos, la dinámica vítrea y la cavitación.

Metodología
Los autores construyeron un aparato experimental para vibrar verticalmente un contenedor cuadrado ($W = 10.9$ cm) lleno de agua desionizada ($h = 5$ mm) a una profundidad de $h = 5$ mm. El contenedor experimenta una oscilación sinusoidal $Z(t) = A \sin(2\pi ft)$, con frecuencias de excitación $f \in [20, 95]$ Hz y amplitudes $A \in [0.02, 0.2]$ mm.

• Partículas: Se utilizaron esferas de poliestireno de celda cerrada con una densidad $\rho_p \approx 15$ kg m$^{-3}$ y diámetros medios de $d = 1.4 \pm 0.2$ mm y $d = 2.1 \pm 0.2$ mm. El ángulo de contacto medido fue de $87^\circ \pm 6^\circ$.
• Empaquetamiento: Se formó una monocapa estricta, con fracciones de empaquetamiento $\phi$ definidas por la relación entre el área proyectada de las partículas y el área del contenedor, probadas hasta $\phi \approx 0.90$.
• Imagen: El patrón de formación se monitoreó mediante una cámara a 73 Hz ($2448 \times 2048$ píxeles). Las trayectorias de las partículas y la dinámica intra-ciclo se capturaron utilizando una cámara de alta velocidad a 460 Hz ($2448 \times 1024$ píxeles) y 186 Hz para el análisis específico de trayectorias.
• Análisis: Las posiciones de las partículas se extrajeron mediante la transformada de Hough circular. Las relaciones de dispersión se determinaron mediante transformadas de Fourier 2D de imágenes promediadas en el tiempo. La dinámica de las partículas se cuantificó mediante el Desplazamiento Cuadrático Medio (MSD), distribuciones de probabilidad de desplazamiento y la susceptibilidad de cuatro puntos ($\chi_4$) para medir la heterogeneidad dinámica.

Contribuciones Clave y Resultados
Al variar sistemáticamente $\phi$, $f$ y $A$, el estudio identifica cinco regímenes dinámicos distintos, resumidos en un diagrama de fases:

1. Ondas Estacionarias Regulares: En frecuencias y amplitudes moderadas, la balsa exhibe patrones clásicos de ondas de Faraday (cuadrados, rayas) que oscilan a la mitad de la frecuencia de excitación ($f/2$).

   • Relaciones de Dispersión: Los autores midieron la dispersión de las ondas y encontraron que, en fracciones de empaquetamiento más altas ($\phi \approx 0.90$), la relación de dispersión clásica falla. Una relación modificada que tiene en cuenta una tensión superficial efectiva ($\gamma_e$) y un módulo de flexión fuera del plano ($B$) ajusta mejor los datos.
   • Elasticidad Emergente: A medida que $\phi$ aumenta, $\gamma_e$ disminuye monótonamente, mientras que $B$ aumenta rápidamente. Esto indica que las partículas densamente empaquetadas desarrollan una red de fuerzas de contacto persistente, impartiendo una rigidez de flexión significativa a la interfaz.

2. Ondas Caóticas: En frecuencias pequeñas a intermedias con un aumento de la amplitud, el sistema transita de ondas estacionarias ordenadas a patrones espaciotemporales caóticos. La presencia de partículas puede distorsionar las ondas estacionarias, proporcionando vías adicionales para la formación de defectos y el caos en comparación con sistemas de líquidos puros.

3. Formación de Vacíos: En fracciones de empaquetamiento bajas ($\phi = 0.77$) y amplitudes bajas, la cobertura insuficiente de partículas conduce a la formación de vacíos que coexisten con las ondas estacionarias.

4. Régimen Vítreo: En frecuencias altas ($f \ge 65$ Hz) y amplitudes bajas, las ondas estacionarias desaparecen. Las partículas exhiben un movimiento de tipo "térmico" dentro de jaulas formadas por sus vecinas.

   • Subdifusión: El análisis de MSD revela una transición de una difusión normal ($\alpha \approx 1$) a frecuencias de empaquetamiento más bajas, hacia un comportamiento subdifusivo ($\alpha \approx 0.5$) a $\phi$ más altos, característico de líquidos sobreenfriados que se aproximan a una transición vítrea.
   • Heterogeneidad: El sistema muestra heterogeneidad dinámica espacial y temporal. La susceptibilidad de cuatro puntos $\chi_4$ alcanza un máximo en un tiempo de retraso $\tau^*$, el cual aumenta exponencialmente con $\phi$. El valor máximo de $\chi_4$ varía de forma no monótona, lo que sugiere una competencia entre la descorrelación inducida por la onda y el encajonamiento impulsado por la densidad.

5. Formación de Cavidades: En frecuencias altas ($f \ge 75$ Hz) y amplitudes críticas ($A^*$), una gran cavidad se nuclea y se expande dentro de la balsa.

   • Mecanismo: La cavidad contiene ondas de superficie libre, mientras que la balsa circundante permanece en un estado vítreo. El tamaño de la cavidad se estabiliza mediante un equilibrio entre la presión de cierre de la balsa densa y una deriva hacia afuera de las partículas impulsada por las interacciones onda-partícula (migración lejos de los antinodos de la onda).
   • Escalamiento: La amplitud crítica $A^*$ requerida para la cavitación disminuye con el aumento de la frecuencia. El tamaño de la cavidad estabilizada $L^*$ escala linealmente con una energía inyectada característica $E^* = A^{*2}f^2/(gd)$.
   • Transición de Forma: Por debajo de $f \approx 85$ Hz, las cavidades son irregulares; por encima de esta frecuencia, se vuelven circulares y estables, indicando una transición de una resistencia de tipo líquido a una de tipo elástico por parte de la balsa circundante.

Significancia y Reivindicaciones
El artículo sostiene que las balsas de partículas vibradas verticalmente sirven como un modelo versátil para investigar la física fuera del equilibrio.

• Mecánica Interfacial: El estudio demuestra que las balsas de partículas densas pueden exhibir elasticidad emergente (módulo de flexión) y tensión superficial efectiva que son sintonizables mediante la fracción de empaquetamiento, distintas de las interfaces de fluidos puros.
• Dinámica Vítrea: El sistema proporciona una plataforma macroscópica y observable para estudiar las transiciones vítreas, el transporte subdifusivo y la heterogeneidad dinámica en sistemas atermales impulsados por ruido externo.
• Interacciones Complejas: A diferencia de los potenciales idealizados (por ejemplo, Lennsen-Jones), las interacciones aquí surgen de efectos capilares e hidrodinámicos de muchos cuerpos complejos. La capacidad de controlar estas interacciones mediante parámetros de vibración ofrece una vía para comprender los ensamblajes de partículas densas en sistemas del mundo real más allá del simple comportamiento elastoplástico.

Los autores mantienen un alcance modesto, señalando que el comportamiento transitorio en los regímenes caóticos y la complejidad total de las interacciones partícula-fluido permanecen más allá del alcance del presente trabajo, centrándose en cambio en la caracterización de los comportamientos de estado estacionario de los regímenes identificados.