Optimal existence of weak solutions for the generalised Navier-Stokes-Voigt equations

Resumen Técnico: Existencia Óptima de Soluciones Débiles para las Ecuaciones Generalizadas de Navier-Stokes-Voigt

1. Planteamiento del Problema

El artículo investiga las ecuaciones de Navier-Stokes-Voigt (NSV) generalizadas para fluidos incompresibles en un dominio acotado $\Omega \subset \mathbb{R}^d$ ($d \ge 2$) durante un intervalo de tiempo $[0, T]$. El modelo describe fluidos viscoelásticos no newtonianos que siguen una ley de potencia.

El sistema de ecuaciones gobernante es:
$$\begin{cases} \partial_t(v - \kappa\Delta v) + \nabla \cdot (v \otimes v) + \nabla \pi - \nu\nabla \cdot (|D(v)|^{p-2}D(v)) = f, & \text{en } Q_T, \\ \nabla \cdot v = 0, & \text{en } Q_T, \\ v = v_0, & \text{en } \Omega \times \{0\}, \\ v = 0, & \text{en } \partial\Omega \times (0, T), \end{cases}$$
donde:

• $v$ es el campo de velocidad y $\pi$ la presión.
• $\kappa > 0$ es el tiempo de relajación (término de Voigt) y $\nu > 0$ la viscosidad cinemática.
• $D(v) = \frac{1}{2}(\nabla v + (\nabla v)^\top)$ es la parte simétrica del gradiente de velocidad.
• $p \in (1, \infty)$ es el exponente de la ley de potencia que caracteriza el comportamiento del fluido (shear-thinning si $1<p<2$, shear-thickening si $p>2$).
• $f$ es la fuerza externa.

El objetivo principal es establecer la existencia y unicidad de soluciones débiles para rangos de $p$ que sean óptimos, es decir, los más bajos posibles permitidos por la teoría matemática actual.

2. Metodología y Enfoque Analítico

Los autores emplean una combinación de técnicas avanzadas de análisis funcional y teoría de ecuaciones en derivadas parciales (EDP):

• Aproximación de Galerkin: Se construyen soluciones aproximadas en subespacios de dimensión finita utilizando una base de eigenfunciones.
• Estimaciones Uniformes: Se derivan estimaciones de energía a priori para controlar las normas de las soluciones aproximadas en espacios de Sobolev y Lebesgue.
• Teoría de Operadores Monótonos: Para ciertos rangos de $p$, se utiliza la teoría estándar de operadores monótonos (técnica de Minty-Browder) para identificar el límite no lineal.
• Descomposición de Presión (Técnica Clave): Para los casos más difíciles (donde la inmersión de Sobolev estándar falla), los autores utilizan una descomposición de la presión inspirada en trabajos de Wolf y Diening. La presión se divide en componentes ($\pi_1, \pi_2, \pi_h$) asociadas a términos inerciales, viscosos y convectivos, permitiendo manejar la falta de regularidad en el campo de velocidad.
• Lema de Compacidad de Aubin-Dubinski˘ı: Se aplica para extraer convergencia fuerte de las soluciones aproximadas, esencial para pasar al límite en los términos no lineales.
• Evitación de la Proyección de Helmholtz-Hodge: Una innovación metodológica notable es que los autores no utilizan la proyección de Helmholtz-Hodge. Esta proyección suele restringir el análisis a $p \ge 2$. Al evitarla, pueden tratar directamente el rango $p > 1$ para dimensiones bajas, superando obstáculos técnicos asociados al control del término convectivo en $p < 2$.

3. Contribuciones Clave y Novedades

1. Rangos Óptimos de Existencia:
   • Para dimensiones $2 \le d \le 3$, se establece la existencia y unicidad para $p > 1$. Este es el límite inferior óptimo, ya que $p>1$ es necesario para la coercividad del operador $p$-Laplaciano.
   • Para dimensión $d = 4$, se establece la existencia para $p > \frac{2d}{d+2}$ (específicamente $p > 4/3$). Este es un resultado de frontera, ya que para $d > 3$ y $p$ muy pequeño, la falta de un triple de Gelfand adecuado impide la aplicación directa de lemas de compacidad estándar.
2. Tratamiento de Dimensiones Superiores ($d=4$): El artículo aborda analíticamente el caso $d=4$, que es crítico debido al equilibrio delicado entre la dimensión espacial y el crecimiento de la no linealidad. Se demuestra que el marco analítico es robusto incluso en estas dimensiones "límite".
3. Unicidad de Soluciones: Se prueba la unicidad de la solución débil bajo las mismas condiciones de existencia, utilizando estimaciones de energía y la monotonía del operador no lineal.
4. Novedad en la Técnica de Presión: La aplicación refinada de la descomposición de presión en el contexto de NSV generalizados permite compensar la pérdida de compacidad en dimensiones altas y rangos bajos de $p$.

4. Resultados Principales

• Teorema de Existencia (Teoremas 3.2 y 3.3):
  • Si $2 \le d \le 3$ y $p > 1$, existe al menos una solución débil $v \in C([0, T]; V) \cap L^p(0, T; V_p)$.
  • Si $d = 4$ y $\frac{2d}{d+2} < p < \frac{d}{2}$ (es decir, $4/3 < p < 2$), existe al menos una solución débil bajo condiciones de regularidad del dominio ($C^2$) y datos iniciales adecuados.
• Teorema de Unicidad (Teorema 3.4): Bajo las mismas hipótesis, la solución débil es única.
• Convergencia: Se demuestra que las soluciones aproximadas convergen fuertemente en espacios adecuados, lo que permite identificar el límite del término no lineal $|D(v)|^{p-2}D(v)$.

5. Significado e Impacto

Este trabajo representa un avance significativo en la teoría matemática de fluidos viscoelásticos no newtonianos:

• Optimalidad: Cierra la brecha en la literatura sobre los rangos de $p$ para los cuales se garantiza la existencia de soluciones en el modelo NSV generalizado, alcanzando el límite inferior teórico ($p>1$) en dimensiones físicas relevantes ($d=2,3$).
• Generalización: Extiende la teoría más allá de los fluidos newtonianos ($p=2$) y los modelos de ley de potencia estándar, incorporando efectos elásticos (término de Voigt) de manera rigurosa.
• Aplicabilidad Teórica: Aunque las aplicaciones físicas suelen limitarse a $d \le 3$, el análisis en $d=4$ valida la solidez estructural del método matemático, ofreciendo herramientas para futuros estudios en dimensiones superiores o en problemas de regularidad crítica.
• Metodología: La demostración de que se puede evitar la proyección de Helmholtz-Hodge para tratar $p < 1$ (en realidad $p>1$ pero $<2$) abre nuevas vías para el análisis de sistemas de fluidos con no linealidades complejas donde la proyección estándar falla o es difícil de controlar.

En resumen, el artículo proporciona un marco matemático completo y óptimo para la existencia y unicidad de soluciones en el modelo de Navier-Stokes-Voigt con viscosidad de ley de potencia, resolviendo problemas abiertos en dimensiones críticas y rangos de exponentes bajos.