On three-dimensional flows of thermo-viscoelastic fluids of Giesekus type

Aquí se presenta un resumen técnico detallado del artículo "On three-dimensional flows of thermo-viscoelastic fluids of Giesekus type" de Miroslav Bulíček, Tomáš Los y Jakub Woźnicki.

1. Planteamiento del Problema

El artículo aborda el análisis matemático del flujo de un fluido viscoelástico incompresible, conductor del calor, de tipo Giesekus en un dominio acotado tridimensional $\Omega \subset \mathbb{R}^3$ durante un intervalo de tiempo $(0, T)$.

El sistema acopla la conservación de la masa, el momento lineal, la energía interna y la evolución del esfuerzo elástico (representado por el tensor de Cauchy–Green izquierdo $B$ o el gradiente de deformación $F$). Las ecuaciones gobernantes son:

• Incompresibilidad: $\text{div } v = 0$.
• Balance de Momento: $\rho \partial_t v + \rho \text{div}(v \otimes v) - \text{div } T = \rho f$.
• Evolución del Esfuerzo: Una ecuación de evolución para el tensor $B$ (o $F$) que involucra la derivada de Oldroyd convectada superior y un término disipativo específico del modelo Giesekus.
• Balance de Energía Interna: $\rho \partial_t e + \rho \text{div}(ev) + \text{div } q - T : D = 0$.

Características Clave del Modelo:

• Dependencia de la Temperatura: Todos los parámetros materiales (viscosidad $\nu$, tiempo de relajación $\tau$, conductividad térmica $\kappa$ y la función del módulo elástico $g$) dependen de la temperatura $\theta$.
• Disipación Giesekus: El modelo incluye un término disipativo no lineal en la evolución del esfuerzo, lo que lo distingue del modelo Oldroyd-B.
• Condiciones de Contorno: Condiciones de deslizamiento de Navier para la velocidad ($v \cdot n = 0, (Tn)_\tau = -\gamma^* v_\tau$) y flujo de calor nulo ($\nabla \theta \cdot n = 0$).
• Datos Iniciales: Se consideran datos grandes, lo que significa que no se imponen suposiciones de pequeñez sobre la velocidad inicial, la temperatura o el gradiente de deformación.

2. Metodología

Los autores emplean un marco débil-fuerte sofisticado combinado con herramientas avanzadas de compacidad para manejar las fuertes no linealidades y la falta de regularidad en las soluciones.

A. Consistencia Termodinámica y Reformulación

• Entropía y Energía: Los autores derivan el sistema desde un marco termodinámicamente consistente utilizando la energía libre de Helmholtz. Reemplazan la ecuación de energía interna con una desigualdad de entropía y una desigualdad de energía total. Este enfoque es crucial porque la presión $p$ podría no estar bien definida en sentido débil para ciertas condiciones de contorno, pero la desigualdad de energía permanece válida.
• Transformación de Variables: Para manejar la ecuación de evolución de $B$, que contiene términos no bajo derivadas (lo que dificulta la identificación de límites), los autores introducen una variable $F$ tal que $B = FF^T$. Reformulan la ecuación de evolución del esfuerzo en términos de $F$ (una ecuación de tipo transporte-difusión), siguiendo la estrategia de Masmoudi. Esto permite un mejor control de los términos no lineales.

B. Esquema de Aproximación

La prueba de existencia se basa en un procedimiento de aproximación multietapa que involucra:

1. Aproximación de Galerkin: Proyección de la velocidad, el gradiente de deformación y la energía interna sobre subespacios de dimensión finita.
2. Regularización: Introducción de términos de difusión artificial ($\varepsilon_4 \Delta F$, $\varepsilon_7 \Delta e$) y funciones de corte ($\Lambda_{\varepsilon_3}$, $g_{\varepsilon_1}$) para asegurar la existencia de soluciones aproximadas suaves.
3. Límites Iterativos: Los autores pasan al límite en una secuencia específica de parámetros ($m, n, l \to \infty$ para los índices de Galerkin, seguido de $\varepsilon_7, \varepsilon_4, \varepsilon_2, \varepsilon_1, \varepsilon_5, \varepsilon_6, \varepsilon_3 \to 0$).

C. Herramientas de Compacidad y Convergencia

El núcleo del análisis reside en establecer la convergencia fuerte de la temperatura $\theta$ y del gradiente de deformación $F$, lo cual es necesario para pasar al límite en las relaciones constitutivas no lineales.

• Convergencia de la Temperatura: Los autores utilizan el Lema Div-Curl (Murat-Tartar) combinado con la ecuación de entropía. Mediante la construcción de funciones de prueba específicas que involucran la entropía $\eta$ y la temperatura, demuestran que una subsucesión de temperaturas converge casi por todas partes.
• Convergencia del Gradiente de Deformación: Adaptan la técnica de renormalización (DiPerna-Lions) y una desigualción tipo Gronwall específica para la cantidad $|F|^2 - |F|^2$ (donde $F$ es el límite). Esto implica dividir el campo de velocidad en una parte "tipo Stokes" y un resto para manejar los términos convectivos, demostrando finalmente la convergencia fuerte de $F$ en $L^2$.
• Límites de Mordida (Biting Limits): Se utiliza el lema de mordida de Chacon para manejar la convergencia débil en espacios $L^1$ para cantidades como los términos de disipación.

3. Contribuciones Clave

1. Existencia Global en el Tiempo con Datos Grandes: El artículo establece la existencia de una solución débil global para el modelo completo 3D de Giesekus con dependencia de la temperatura, sin requerir:
   • Pequeñez de los datos iniciales.
   • Difusión artificial de esfuerzos (una regularización común en la literatura anterior).
   • Suposiciones antinaturales sobre el crecimiento de los parámetros materiales (por ejemplo, se permite que la conductividad térmica sea acotada, a diferencia de trabajos anteriores que requerían crecimiento polinomial).
2. Solución Débil Termodinámicamente Consistente: La solución se define en un sentido "débil disipativo", satisfaciendo la desigualdad de energía total y la desigualdad de entropía. Esto asegura que la solución respete la segunda ley de la termodinámica.
3. Manejo de la No Linealidad Giesekus: Los autores tratan exitosamente el término de disipación no lineal específico del modelo Giesekus ($\alpha(B-I)$ y $\tau(B^2-B)$) en presencia de acoplamiento térmico, un problema que permaneció abierto para datos grandes en 3D.
4. Parámetros Materiales Generales: El resultado se mantiene para una amplia clase de funciones materiales continuas (viscosidad, tiempo de relajación, etc.) que satisfacen únicamente condiciones naturales de acotamiento y monotonía, haciendo que el modelo sea aplicable a una amplia gama de escenarios físicos.

4. Resultados Principales

Teorema: Bajo la suposición de que el dominio $\Omega$ es Lipschitz y los datos iniciales $(v_0, F_0, \theta_0)$ pertenecen a los espacios funcionales naturales de energía-entropía (específicamente $v_0 \in L^2$, $F_0 \in L^2$, $\theta_0 \in L^1$ con $\det F_0 > 0$), existe una solución débil disipativa global en el tiempo $(v, \theta, F)$ para el sistema.

La solución satisface:

• Estabilidad Débil Secuencial: La solución es estable bajo límites de secuencias aproximantes.
• Desigualdad de Energía: La energía total (cinética + interna) satisface una desigualdad diferencial que cuenta con la disipación en el contorno.
• Desigualdad de Entropía: La producción de entropía es no negativa, satisfaciendo la segunda ley de la termodinámica.
• Regularidad: La solución posee propiedades de regularidad específicas, incluyendo $v \in L^\infty(0,T; L^2) \cap L^2(0,T; W^{1,2})$, $\theta \in L^\infty(0,T; L^1) \cap L^{5/3}(Q_T)$, y $F \in L^\infty(0,T; L^2) \cap L^4(Q_T)$.

5. Significado

Este trabajo representa un avance significativo en la teoría matemática de fluidos complejos. Anteriormente, los resultados de existencia global para fluidos viscoelásticos conductores del calor a menudo dependían de:

• Modelos simplificados (por ejemplo, Oldroyd-B con difusión de esfuerzos).
• Suposiciones de datos pequeños.
• Configuraciones bidimensionales.
• Suposiciones antinaturales sobre el crecimiento de la conductividad térmica.

Al eliminar estas restricciones, Bulíček, Los y Woźnicki proporcionan una base rigurosa para modelar flujos realistas, a gran escala y tridimensionales de materiales complejos (como polímeros y fluidos biológicos) donde los efectos térmicos y los esfuerzos elásticos están fuertemente acoplados. La metodología desarrollada, particularmente la combinación de compacidad basada en entropía y el tratamiento específico de la no linealidad Giesekus sin difusión artificial, establece un nuevo estándar para el análisis de fluidos viscoelásticos de tipo tasa.