Resumen Técnico: Elasto-hidrodinámica de las Interacciones Gotícula-Pozo

Planteamiento del Problema

Si bien la hidrodinámica de una gota impactando un pozo de líquido está bien caracterizada para fluidos newtonianos —resultando en la formación de cavidades, expansión de la corona, ondas capilares y chorros de Worthington—, el comportamiento correspondiente para fluidos elásticos (o de Boger) permanece en gran medida inexplorado. La literatura existente se centra fuertemente en sistemas newtonianos o en el adelgazamiento de filamentos viscoelásticos en geometrías controladas (por ejemplo, reómetros de ruptura capilar). Existe una brecha significativa en la comprensión de la dinámica inestable y espacialmente evolutiva de la formación, colapso y posterior eyección de chorros cuando tanto la gota que impacta como el pozo objetivo son fluidos poliméricos elásticos. Específicamente, la interacción entre las fuerzas inerciales, capilares y elásticas que rigen la evolución espacio-temporal de estos eventos, así como los distintos regímenes morfológicos que surgen, requiere una investigación sistemática.

Metodología

El estudio emplea un tríptico de enfoques experimentales, teóricos y computacionales:

1. Configuración Experimental:

   • Fluidos: Se utilizaron soluciones acuosas de óxido de polietileno de alto peso molecular (PEO, ~5–6 millones de Da) a concentraciones de 0, 250, 500, 1000 y 2000 ppm para crear fluidos elásticos (de Boger). Los fluidos de control incluyeron agua pura y mezclas de agua-glicerol (con viscosidad igualada a la de los fluidos elásticos).
   • Procedimiento: Las gotas se generaron de forma cuasi-estacionaria y se dejaron caer libremente sobre un pozo profundo. Las velocidades de impacto oscilaron entre 1.5 y 4.3 m/s.
   • Diagnóstico: Imágenes de alta velocidad (4000 fps) capturaron la evolución temporal de la expansión de la cavidad, el colapso, la formación de la corona y la eyección del chorro.
   • Parámetros Adimensionales: La dinámica se analizó a través de un rango de números de Weber ($We$) y de Deborah ($De$), con el número de Bond ($Bo$) fijado en 0.27 para asegurar gotas pre-impacto esféricas. El número de Ohnesorge ($Oh$) se mantuvo bajo (<0.0046), haciendo que los efectos viscosos fueran despreciables.

2. Marco Teórico:

   • Evolución de la Cavidad: Se derivó un modelo de conservación de energía dentro de un marco de flujo potencial. La energía cinética inicial de la gota se equilibró contra la energía potencial gravitatoria, la energía superficial, la energía cinética del pozo y un término para la energía elástica almacenada.
   • Adelgazamiento del Chorro: Se desarrolló un modelo teórico para la evolución del radio del chorro de Worthington utilizando el modelo constitutivo FENE-P (Finitely Extensible Nonlinear Elastic - Peterlin). Este modelo tiene en cuenta la extensibilidad finita de las cadenas poliméricas y el equilibrio entre los esfuerzos elásticos y capilares en el régimen elasto-capilar.

3. Simulaciones Numéricas:

   • Solucionador: Las simulaciones se realizaron utilizando COMSOL Multiphysics.
   • Ecuaciones Gobernantes: Se acopló el modelo de campo de fase Cahn-Hilliard con la ecuación constitutiva FENE-P para rastrear la interfaz y modelar los esfuerzos elásticos.
   • Validación: Las simulaciones se validaron con datos experimentales tanto para casos newtonianos como elásticos, centrándose en la morfología de la cavidad, la altura del chorro y los campos de velocidad.

Principales Contribuciones y Resultados

1. Dominancia de la Elasticidad sobre la Viscosidad

Los experimentos comparativos entre fluidos elásticos (2000 ppm de PEO) y fluidos newtonianos con viscosidad igualada (agua-glicerol) demostraron que los cambios morfológicos observados son impulsados puramente por la elasticidad del fluido, no por la viscosidad.

• Morfología de la Cavidad: Mientras que los fluidos newtonianos producen cavidades hemisféricas, los fluidos elásticos exhiben una forma de cavidad trapezoidal distintiva durante la fase de retracción. Esto se atribuye al estiramiento no uniforme de los polímeros y a la distribución de esfuerzos elásticos resultante, lo que retrasa la propagación de las ondas capilares y aplana el fondo de la cavidad.
• Independencia de la Viscosidad: El estudio confirma que, para el rango probado, la viscosidad juega un papel insignificante en el proceso de ruptura, mientras que la distribución del esfuerzo elástico dicta la forma de la cavidad y la dinámica del chorro.

2. Conversión y Almacenamiento de Energía

El modelo de conservación de energía cuantificó la conversión de la energía cinética en energía elástica.

• Almacenamiento Elástico: El estudio indica que aproximadamente el 35% de la energía de impacto cinética inicial de la gota puede almacenarse como energía elástica mediante el estiramiento de las cadenas poliméricas durante la expansión de la cavidad en casos específicos de alta concentración. Específicamente, para el caso 500D-500P, la reducción del factor de partición de energía ($\beta$) de 0.80 (Newtoniano) a 0.45 (elástico) corresponde a que aproximadamente el 35% de la energía cinética inicial se almacena como energía elástica.
• Impacto en la Dinámica: Este almacenamiento de energía reduce la energía disponible para la expansión de la cavidad, lo que conduce a tamaños máximos de cavidad más pequeños y tasas de crecimiento más lentas en fluidos elásticos en comparación con sus contrapartes newtonianas. El modelo predijo con éxito la evolución temporal del radio medio de la cavidad ajustando el factor de partición de energía ($\beta$).

3. Mapas de Regímenes y Transiciones Morfológicas

El estudio estableció cuatro regímenes distintos en el espacio de parámetros $[De, We]$:

1. Pinch-off (Ruptura): Desprendimiento completo de gotas secundarias (observado principalmente en sistemas newtonianos).
2. Gotícula Satélite + Filamento: Formación de una gotícula satélite conectada por un filamento. Esto incluye:
   • Cuentas en un Cordel (Beads-on-a-String, BOAS): Observado cuando las tasas de crecimiento de la inestabilidad de Rayleigh son comparables o más lentas que el inverso del tiempo de relajación, permitiendo que se desarrolle la asimetría antes de que el estiramiento polimérico domine.
   • Filamentos Delgados/Gruesos: Observado cuando los esfuerzos elásticos equilibran la presión capilar, evitando la ruptura completa.
3. Sin Chorro: Supresión de la formación del chorro a bajos $We$ o altos $De$ debido al retraso en la convergencia de las ondas capilares.
4. Sin Pinch-off: Formación de un chorro de Worthington estable que no se rompe, que ocurre a altos $De$, donde los esfuerzos elásticos contrarrestan continuamente el estrechamiento capilar.

4. Adelgazamiento Elasto-Capilar del Chorro

El análisis teórico y experimental del adelgazamiento del chorro de Worthington reveló:

• Decaimiento Exponencial: En el régimen elasto-capilar, el radio del chorro decae exponencialmente con el tiempo, siguiendo la ley de escala $r(t) \sim \exp(-t/2\lambda)$.
• Diferencia de Escala: Esto difiere de la ley de decaimiento clásica $r(t) \sim \exp(-t/3\lambda)$ observada en el adelgazamiento de filamentos viscoelásticos estándar, atribuido a las condiciones de contorno específicas y a las fuerzas de tensión axial en el chorro de Worthington.
• Estabilidad: Una mayor elasticidad suprime la inestabilidad de Rayleigh-Plateau, lo que conduce a filamentos más gruesos y estables, y retrasa o evita el pinch-off.

5. Perspectivas Numéricas sobre los Campos de Esfuerzo y Velocidad

Las simulaciones proporcionaron la visualización de los campos de flujo internos no accesibles experimentalmente:

• Distribución de Esfuerzos: El esfuerzo elástico ($\tau_{zz}$) es máximo en el eje de simetría del chorro y disminuye hacia la interfaz, lo que corresponde a la mayor velocidad axial y al mayor estiramiento polimérico en el núcleo.
• Supresión de la Velocidad: En fluidos elásticos, la velocidad de expansión de la cavidad es aproximadamente la mitad que la del agua pura debido a los esfuerzos elásticos opuestos. En consecuencia, la entalpía de aire se reduce significamente y no se observa un vórtice dentro de la cavidad durante la retracción (a diferencia de los casos newtonianos).

Significado y Reivindicaciones

Los autores afirman que este trabajo avanza significativamente en la comprensión del paradigma inexplorado de la interacción inercia-capilaridad-elasticidad en las interacciones gota-pozo.

• Física Fundamental: El estudio aísla los efectos elásticos de los efectos viscosos, demostrando que la elasticidad del fluido por sí sola puede alterar fundamentalmente la morfología de la cavidad (trapezoidal vs. hemisférica) y suprimir la ruptura del chorro.
• Perspectiva Cuantitativa: Al derivar un modelo teórico para la evolución del radio de la cavidad y el adelgazamiento del chorro, el artículo proporciona un método para estimar el tiempo de relajación extensional de los fluidos elásticos directamente a partir de estudios de interacción gota-pozo.
• Definición de Regímenes: La identificación de cuatro regímenes morfológicos distintos ofrece un marco predictivo para los resultados del impacto de gotas en fluidos elásticos, lo cual es crucial para aplicaciones que involucran fluidos no newtonianos como pinturas, tintas y fluidos biológicos.

El artículo concluye que la conversión de la energía cinética en energía elástica es el mecanismo principal que gobierna la supresión de la ruptura de la corona, la reducción del tamaño de la cavidad y la estabilización de los chorros de Worthington en fluidos elásticos.