Mezcla en fluidos viscoelásticos utilizando turbulencia elástica

Información Básica

• ID del Artículo: 2409.06391
• Título: Mixing in viscoelastic fluids using elastic turbulence
• Autores: Reinier van Buel, Holger Stark (Technische Universität Berlin)
• Clasificación: physics.flu-dyn
• Fecha de Publicación: Septiembre de 2024
• Enlace del Artículo: https://arxiv.org/abs/2409.06391

Resumen

Este artículo investiga el efecto de la turbulencia elástica en la mezcla del campo de concentración escalar dentro de fluidos viscoelásticos en una geometría bidimensional de Taylor-Couette, utilizando soluciones numéricas del modelo Oldroyd-B. El estado del flujo se determina mediante un parámetro de orden de flujo secundario, que indica la región de turbulencia elástica. Cuando se comienza desde un estado turbulento y se reduce posteriormente el número de Weissenberg, aparece un flujo débilmente caótico por debajo de Wic. La advección en estados de turbulencia elástica y flujo débilmente caótico induce mezcla, ilustrada mediante la evolución temporal de la desviación estándar de la concentración de soluto respecto a una distribución uniforme. En particular, la mezcla es fuerte en el estado de turbulencia elástica, cuantificada mediante la velocidad de mezcla, el tiempo de mezcla y la eficiencia de mezcla. Estas tres cantidades obedecen leyes de escala. Importantemente, el parámetro de orden se correlaciona fuertemente con la velocidad de mezcla, siendo por lo tanto un buen indicador de mezcla dentro del fluido.

Antecedentes de Investigación y Motivación

Contexto del Problema

1. Desafíos de mezcla a bajo número de Reynolds: En fluidos newtonianos, bajo condiciones de bajo número de Reynolds, el flujo es laminar y el proceso de mezcla es lento, controlado principalmente por difusión. Esta es una limitación importante en dispositivos de microfluídica.
2. Descubrimiento de turbulencia elástica: La adición de polímeros de alto peso molecular al fluido lo hace viscoelástico, permitiendo que exhiba fenómenos de turbulencia elástica incluso cuando la inercia es despreciable.
3. Importancia de la mezcla a microescala: En dispositivos de laboratorio en chip (lab-on-a-chip), la mezcla eficiente de fluidos es crucial, particularmente a escala micrométrica.

Motivación de la Investigación

• Explorar el potencial de la turbulencia elástica para mejorar la eficiencia de mezcla
• Comprender la naturaleza física de los mecanismos de mezcla en fluidos viscoelásticos
• Proporcionar orientación teórica para aplicaciones de microfluídica

Contribuciones Principales

1. Determinación de la naturaleza de la transición: Identificar la transición de flujo laminar a turbulencia elástica como una transición subcrítica y descubrir la existencia de un estado de flujo débilmente caótico
2. Establecimiento de métodos de caracterización de mezcla: Introducir tres indicadores cuantitativos —velocidad de mezcla, tiempo de mezcla y eficiencia de mezcla— que obedecen leyes de escala
3. Descubrimiento de fuerte correlación entre parámetro de orden y mezcla: Demostrar que el parámetro de orden de flujo secundario se correlaciona fuertemente con la velocidad de mezcla, sirviendo como un buen indicador del grado de mezcla del fluido
4. Verificación del mecanismo de mezcla de Batchelor: Observar decaimiento exponencial en el estado de turbulencia elástica, correspondiente a la región de mezcla de Batchelor, consistente con resultados experimentales

Detalles Metodológicos

Configuración Geométrica y Ecuaciones de Control

El estudio utiliza una geometría bidimensional de Taylor-Couette que contiene dos cilindros concéntricos:

• Radio del cilindro externo: $r_o = 10\mu m$
• Radio del cilindro interno: $r_i = 2.5\mu m$
• El cilindro externo rota en sentido antihorario con velocidad angular $\Omega = 2\pi s^{-1}$

El movimiento del fluido se describe mediante la ecuación de Navier-Stokes generalizada: $\rho\left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p + \eta_s \nabla^2 \mathbf{u} + \nabla \cdot \boldsymbol{\tau}$

Modelo Constitutivo Oldroyd-B

El tensor de esfuerzo polimérico $\boldsymbol{\tau}$ se modela utilizando el modelo Oldroyd-B: $\boldsymbol{\tau} + \lambda \stackrel{\nabla}{\boldsymbol{\tau}} = \eta_p[\nabla \otimes \mathbf{u} + (\nabla \otimes \mathbf{u})^T]$

donde $\stackrel{\nabla}{\boldsymbol{\tau}}$ es la derivada convectiva superior.

Modelado de Mezcla

La evolución del campo de concentración escalar $c(\mathbf{r}, t)$ se rige por la ecuación de advección-difusión: $\frac{\partial c}{\partial t} + (\mathbf{u} \cdot \nabla)c = D\nabla^2 c$

En este estudio se establece $D = 0$, ignorando la difusión molecular.

Parámetros Adimensionales Clave

• Número de Reynolds: $Re = \rho\Omega r_o^2/\eta_s = 0.628 \times 10^{-3}$
• Número de Weissenberg: $Wi = \lambda\Omega$
• Razón de viscosidad: $\beta = \eta_p/\eta_s = 1.5$

Configuración Experimental

Método Numérico

• Utilizar el solucionador de volumen finito de código abierto OpenFOAM®
• Emplear el solucionador especializado RheoTool para procesar flujos viscoelásticos
• Malla: malla de radios con 100 elementos en dirección radial y 120 elementos en dirección angular
• Paso de tiempo: $\delta t = 10^{-5}s$

Definición del Parámetro de Orden

Se define la intensidad de flujo secundario como parámetro de orden: $\Phi = \sigma = \sqrt{\langle[\mathbf{u}(\mathbf{r},t) - \mathbf{u}_0(\mathbf{r})]^2\rangle}/u_0^{\max}$

Indicadores de Caracterización de Mezcla

Desviación estándar normalizada: $\Delta c = \sqrt{\frac{\langle[c(\mathbf{r},t) - c_\infty]^2\rangle}{\langle[c(\mathbf{r},0) - c_\infty]^2\rangle}}$

donde $c_\infty$ es la concentración en estado completamente mezclado.

Resultados Experimentales

Transición de Estados de Flujo

1. Transición subcrítica: Se descubre que la transición de flujo laminar a turbulencia elástica presenta naturaleza subcrítica, con existencia de histéresis
2. Número de Weissenberg crítico:
   • Rama ascendente: $Wi_c = 9.95$
   • Rama descendente: $Wi_c^I = 9.26$
3. Estado débilmente caótico: Se observa flujo débilmente caótico por debajo de $Wi_c^I$

Análisis de Características de Mezcla

Estado de Flujo Débilmente Caótico ($5.03 < Wi \leq 8.17$)

• La desviación estándar decae según ley de potencia con exponente cercano a $-1/6$
• El espectro de potencia muestra escalado de ley de potencia con exponente de $-5$ o mayor

Estado de Turbulencia Elástica ($Wi \geq 9.42$)

• Decaimiento inicial según ley de potencia (exponente aproximadamente $1/5$), seguido de decaimiento exponencial
• Forma de decaimiento exponencial: $\Delta c \sim \exp(-\alpha_{mix}t)$
• El espectro de potencia presenta características del espectro de Batchelor (escalado $k^{-1}$)

Relaciones de Ley de Escala

1. Tiempo de mezcla: $t_{mix} \sim (Wi - Wi_c^I)^{-0.52}$
2. Concentración en región II: $c_{II}/c_\infty \sim (Wi - Wi_c^I)^{0.28}$
3. Relación de escala: $c_{II} \sim t_{mix}^{-1/2}$

Correlación entre Parámetro de Orden y Velocidad de Mezcla

Los resultados experimentales muestran una fuerte correlación entre el parámetro de orden $\Phi_d$ y la velocidad de mezcla $\alpha_{mix}$, indicando que el aumento de fluctuaciones de velocidad mejora directamente la eficiencia de mezcla dentro del fluido.

Trabajos Relacionados

Investigación en Turbulencia Elástica

• Primer descubrimiento de turbulencia elástica en soluciones poliméricas por Groisman & Steinberg
• Análisis de estabilidad lineal y no lineal en geometría de Taylor-Couette
• Observación de patrones de flujo helicoidal y en bandas en simulaciones numéricas tridimensionales

Investigación en Mejora de Mezcla

• Mezcla por flujo extensional en molinos de cuatro rodillos
• Inestabilidad elástica en geometría de ranura cruzada
• Mejora de transferencia de calor por convección en canales sinuosos

Conclusiones y Discusión

Conclusiones Principales

1. Naturaleza de la transición: Confirmar la naturaleza subcrítica de la transición de turbulencia elástica en flujo de Taylor-Couette
2. Mecanismo de mezcla: Verificar que la turbulencia elástica puede lograr mezcla eficiente a números de Reynolds extremadamente bajos
3. Región de Batchelor: Observar características de mezcla de Batchelor en estado de turbulencia elástica
4. Indicador predictivo: El parámetro de orden puede servir como indicador confiable de eficiencia de mezcla

Limitaciones

1. Restricción bidimensional: La investigación se limita a geometría bidimensional; los efectos tridimensionales pueden producir resultados diferentes
2. Dependencia de malla: El comportamiento de histéresis desaparece en mallas más finas, indicando el impacto de la precisión numérica
3. Sensibilidad a condiciones iniciales: Las características de mezcla dependen de la elección de la distribución de concentración inicial

Direcciones Futuras

1. Efecto de ondas elásticas: Investigar el impacto de ondas elásticas en eficiencia de mezcla, velocidad de mezcla y tiempo de mezcla
2. Extensión tridimensional: Verificar las leyes de escala encontradas en geometría tridimensional
3. Aplicaciones prácticas: Explorar el potencial de aplicación práctica en dispositivos de microfluídica

Evaluación Profunda

Fortalezas

1. Investigación sistemática: Análisis integral de características de mezcla durante la transición de flujo laminar a turbulencia elástica
2. Caracterización cuantitativa: Introducción de múltiples indicadores cuantitativos para caracterizar grado y eficiencia de mezcla
3. Verificación teórica: Los resultados numéricos son consistentes con la teoría de mezcla de Batchelor y observaciones experimentales
4. Valor práctico: Proporciona orientación teórica importante para aplicaciones de microfluídica

Deficiencias

1. Problemas de precisión numérica: La dependencia de malla del comportamiento de histéresis sugiere limitaciones potenciales en el método numérico
2. Comprensión de mecanismos: La explicación del mecanismo físico del decaimiento según ley de potencia en estado débilmente caótico es insuficiente
3. Rango de parámetros: El rango de parámetros estudiado es relativamente limitado, faltando exploración más amplia del espacio de parámetros

Impacto

1. Contribución teórica: Proporciona perspectivas importantes para comprender mecanismos de mezcla en fluidos viscoelásticos
2. Perspectivas de aplicación: Tiene importancia significativa para orientación en diseño de dispositivos de microfluídica
3. Valor metodológico: Los métodos de caracterización establecidos pueden generalizarse a otras geometrías de flujo

Escenarios Aplicables

• Optimización de mezcla en dispositivos de microfluídica
• Control de mezcla en procesos de procesamiento de polímeros
• Mezcla a microescala en aplicaciones biomédicas

Referencias

El artículo cita 57 referencias relacionadas, abarcando múltiples aspectos incluyendo turbulencia elástica, estabilidad de flujos viscoelásticos y teoría de mezcla, proporcionando una base teórica sólida para la investigación.