Beyond binary scission: a generalized three-species cascade breakage model for wormlike micellar solutions

Resumen Técnico: Un modelo de ruptura en cascada de tres especies generalizado para soluciones de micelas filamentosas

Planteamiento del problema
Las soluciones de micelas filamentosas (WLM, por sus siglas en inglés) exhiben comportamientos reológicos complejos, incluyendo viscoelasticidad, adelgazamiento por cizalladura y bandeo de cizalladura, impulsados por la ruptura y recombinación continua de redes micelares autoensambladas. Los marcos teóricos existentes, particularmente los modelos de dos especies ampliamente utilizados (por ejemplo, el modelo de Vasquez-Cook-McKinley o VCM), dependen de una representación binaria gruesa de la población micelar (que típicamente distingue únicamente entre cadenas largas entrelazadas y fragmentos cortos). Aunque son exitosos al capturar ciertos rasgos macroscópicos, estos modelos poseen una limitación fundamental: no logran resolver explícitamente la ruptura en cascada multiescala y los regímenes de relajación continua inherentes a las redes micelares reales. Esta aproximación binaria filtra los estados estructurales intermedios, lo que conduce a una incapacidad para describir espectros de relajación amplios, particularmente el aumento de alta frecuencia en los módulos dinámicos observado en experimentos, y limita la interpretación física de las respuestas transitorias y los mecanismos de bandeo de cizalladura.

Metodología
Para abordar estas limitaciones, los autores desarrollan un modelo de ruptura en cascada de tres especies que introduce un estado micelar intermedio en el marco cinético. El modelo asume que la red micelar consiste en tres estados estructurales discretos:

1. Especie A: Una red de gel (primaria, altamente entrelazada).
2. Especie B: Cadenas largas (estado intermedio).
3. Especie C: Cadenas cortas (fragmentos desentrelazados).

La formulación matemática se basa en un enfoque de teoría cinética utilizando distribuciones de espacio de fase continuas. Los autores derivan ecuaciones de Smoluchowski acopladas para las funciones de densidad de probabilidad de las tres especies, incorporando términos de estiramiento convectivo, difusión y reacción. Las características metodológicas clave incluyen:

• Cinética Jerárquica: El modelo define una vía de cascada donde la red de gel (A) se rompe en cadenas largas (B), que posteriormente se rompen en cadenas cortas (C). También se definen vías de recombinación para asegurar la conservación de la masa.
• Escisión Dependiente del Flujo: Las tasas de ruptura para la red de gel y las cadenas largas están acopladas al flujo macroscópico local mediante un mecanismo de convección de hebras con extensión parcial (PEC). Esto crea un mecanismo de ruptura dependiente del esfuerzo donde la deformación acelera la escisión.
• Ecuaciones Constitutivas: El modelo deriva ecuaciones de evolución adimensionales para las densidades numéricas y los tensores de conformación, acopladas a la ecuación de cantidad de movimiento de Cauchy macroscópica. El esfuerzo total es la superposición de los esfuerzos elásticos de las tres especies.
• Implementación Numérica: Los autores analizan el modelo tanto bajo condiciones de flujo homogéneo (para estudiar la viscoelasticidad lineal y la respuesta de cizalladura estacionaria) como en flujo de Couette cilíndrico inhomogéneo (para investigar el bandeo de cizalladura, el arranque transitorio y la relajación).

Contribuciones Clave y Resultados

1. Respuesta Viscoelástica de Tres Modos:
   Bajo flujo homogéneo, el modelo produce analíticamente una estructura de Maxwell de tres modos discretos. A diferencia de los modelos de dos especies, este marco captura con éxito el aumento de alta frecuencia en los módulos de almacenamiento ($G'$) y pérdida ($G''$) dinámicos. Los autores demuestran que el modo de relajación adicional proporcionado por la especie intermedia (B) y las cadenas cortas (C) mejora significativamente el ajuste a las curvas maestras experimentales de soluciones de CPyCl/NaSal, particularmente en el régimen de alta frecuencia.

2. Curvas de Flujo No Monotónicas y Bandeo de Cizalladura:
   El modelo predice una curva constitutiva no monotónica con una región de pendiente estrictamente negativa, un sello distintivo de la inestabilidad constitutiva. En las simulaciones de flujo de Couette inhomogéneo, esta inestabilidad conduce a la formación espontánea de bandas de cizalladura coexistentes (regiones de alta y baja tasa de cizalladura). El modelo de tres especies revela que la especie intermedia (B) forma una capa de transición distinta en la interfaz de la banda, vinculando la banda de baja cizalladura dominada por la red de gel con la banda de alta cizalladura dominada por las cadenas cortas. Esto proporciona una descripción microestructural de la interfaz de la banda más resuelta de lo que permiten los modelos binarios.

3. Dinámica Transitoria e Histéresis:
   El modelo captura comportamientos transitorios complejos:

   • Sobreimpulso de Esfuerzo (Stress Overshoot): Durante el arranque, la competencia entre la carga elástica y la ruptura en cascada acelerada conduce a pronunciados sobreimpulsos de esfuerzo.
   • Histéresis de Flujo: El modelo exhibe histéresis, donde el estado estacionario final depende de la historia del flujo (arranque desde el reposo frente a reducción de paso desde alta cizalladura).
   • Relajación Multietapa: Tras el cese de la cizalladura, el modelo predice un proceso de relajación de múltiples etapas. En lugar de un decaimiento exponencial único, el esfuerzo se relaja a través de etapas secuenciales: liberación rápida de las cadenas cortas (C), ajuste intermedio de las cadenas largas (B) y recuperación lenta de la red de gel (A). Esto resuelve la "relajación prolongada" observada frecuentemente en experimentos pero difícil de capturar con modelos reducidos.

4. Consistencia Mecánica de Continuo:
   El modelo satisface la relación de Lodge-Meissner ($N_1/\sigma_{r\theta} = \gamma_0$) durante deformaciones escalonadas de gran amplitud, confirmando su consistencia tensorial y validez física bajo deformación finita.

Significancia y Reivindicaciones
El artículo afirma que el marco de tres especies propuesto cierra la brecha entre los modelos de red fenomenológicos discretos y las teorías cinéticas continuas (como la teoría original de reptación-reacción de Cates). Al resolver explícitamente la cinética de cascada que gobierna los estados estructurales intermedios, el modelo proporciona un marco físicamente interpretable para el bandeo de cizalladura de micelas filamentosas.

Los autores sostienen que este enfoque captura la conexión intrínseca entre:

• La evolución microestructural en cascada (A $\to$ B $\to$ C).
• La dinámica de relajación amplia (tres escalas de tiempo distintas).
• La localización del flujo macroscópico (bandeo de cizalladura).

La significancia radica en la capacidad del modelo para mantener la tractabilidad computacional macroscópica mientras recupera los modos de relajación intermedios faltantes. Esto ofrece una base física más completa para interpretar la rica dinámica transitoria y de estado estacionario de los fluidos micelares complejos, ayudando potencialmente en la interpretación de futuros experimentos microfluídicos y aplicaciones industriales que involucren WLM. El artículo no propone nuevos protocolos experimentales, sino que ofrece una herramienta teórica para interpretar mejor los fenómenos reológicos existentes.