Global existence and optimal decay for a three-dimensional penalized Navier--Stokes system with biharmonic damping

Resumen Técnico: Existencia Global y Decaimiento Óptimo para un Sistema de Navier–Stokes Penalizado Tridimensional con Amortiguamiento Biharmónico

Planteamiento del Problema
Este artículo investiga un sistema parabólico tridimensional que sirve como una aproximación hiperviscosa y penalizada de las ecuaciones de Navier–Stokes incompresibles. El modelo aborda los desafíos matemáticos de las ecuaciones clásicas de Navier–Stokes (1.1), particularmente el problema abierto de la regularidad global y la unicidad para datos iniciales arbitrarios en tres dimensiones. El sistema propuesto (1.3) combina tres mecanismos:

1. Difusión viscosa clásica ($-\nu\Delta u$).
2. Amortiguamiento biharmónico (hiperviscoso) ($\beta\Delta^2 u$), que proporciona una disipación de orden superior.
3. Penalización de la divergencia ($-\frac{1}{\varepsilon}\nabla \text{div} u$), que relaja la restricción estricta de incompresibilidad ($\text{div} u = 0$) al imponer una casi-incompresibilidad conforme el parámetro $\varepsilon \to 0$.

Crucialmente, el término de convección no lineal se modifica utilizando una corrección de tipo Temam asimétrica respecto al signo (skew-symmetric) ($N(u) = (u \cdot \nabla)u + \frac{1}{2}u \text{div} u$). Esta corrección compensa la débil compresibilidad introducida por la penalización, restaurando la propiedad de cancelación $L^2$ del término de transporte, la cual se pierde típicamente en formulaciones compresibles.

Metodología
El análisis se basa en una combinación de métodos de energía, teoría de semigrupos y técnicas de división de Fourier, con un enfoque central en la derivación de estimaciones que son uniformes con respecto al parámetro de penalización $\varepsilon > 0$.

• Soluciones Débiles (Teorema 1.3): La existencia global de soluciones débiles para cualquier dato inicial $L^2$ se establece mediante una aproximación de tipo Galerkin de Friedrichs. Los autores proyectan la ecuación sobre bolas de Fourier utilizando operadores de corte $J_n$ en lugar del proyector de Leray. Esta distinción es crítica; el uso del proyector de Leray anularía el término de penalización, alterando el problema límite. Se derivan estimaciones a priori uniformes a partir de la identidad de energía, y la compacidad se obtiene mediante el lema de Aubin–Lions–Simon.
• Soluciones Fuertes (Teorema 1.4): Para datos iniciales pequeños en $H^2(\mathbb{R}^3)$, se prueban soluciones fuertes globales utilizando un argumento de punto fijo sobre la formulación débil asociada al semigrupo lineal $S_\varepsilon(t) = e^{tL_\varepsilon}$. El operador lineal $L_\varepsilon = \nu\Delta - \beta\Delta^2 + \frac{1}{\varepsilon}\nabla\text{div}$ genera un semigrupo con propiedades de suavizado. La prueba depende de una estimación de energía $H^2$ uniforme donde los términos no lineales son absorbidos por la disipación viscosa y biharmónica.
• Estimaciones de Decaimiento (Teorema 1.6): Se derivan tasas de decaimiento de tiempo largo óptimas utilizando el método de división de Fourier (Schonbek). El análisis separa las frecuencias bajas y altas. El comportamiento de baja frecuencia se controla mediante la fórmula de Duhamel y la descomposición específica de la no linealidad, mientras que las frecuencias altas se gestionan mediante desigualdades de energía diferenciadas e interpolación de Gagliardo–Nirenberg.

Contribuciones Clave y Resultados

1. Existencia Global de Soluciones Débiles: El artículo demuestra la existencia global de soluciones débiles para cualquier dato inicial $u_0 \in L^2(\mathbb{R}^3)$. Las soluciones satisfacen clases de regularidad específicas, incluyendo $L^\infty(0, T; L^2)$, $L^2(0, T; H^2)$ y $L^4(0, T; L^4)$.
2. Soluciones Fuertes Globales con Pequeñez Uniforme: Para datos iniciales suficientemente pequeños en $H^2(\mathbb{R}^3)$, el sistema admite una solución fuerte global única. Una contribución significativa es que el umbral de pequeñez $\delta_0$ y las estimaciones de energía asociadas son independientes de $\varepsilon$. Esta uniformidad asegura que los resultados se apliquen a toda la familia de sistemas penalizados, no solo a una instancia fija.
3. Tasas de Decaimiento Óptimas: Para datos iniciales en $L^1(\mathbb{R}^3) \cap H^2(\mathbb{R}^3)$, los autores derivan estimaciones de decaimiento de tiempo largo óptimas:
   $$\|u(t)\|_{L^2} \leq C(1+t)^{-3/4}$$
   $$\|\nabla^k u(t)\|_{L^2} \leq C(1+t)^{-3/4 - k/2}, \quad k=1, 2$$
   Estas tasas coinciden con el decaimiento asintótico de la ecuación del calor clásica en tres dimensiones. El análisis confirma que, si bien el término biharmónico domina la disipación de alta frecuencia, el decaimiento a largo plazo está gobernado por el comportamiento de flujo de calor de las bajas frecuencias.
4. Uniformidad en $\varepsilon$: Una característica definitoria del análisis es que todas las constantes y estimaciones permanecen uniformes cuando $\varepsilon \to 0$. El término de penalización $\frac{1}{\varepsilon}\|\text{div} u\|^2$ se trata como un término de disipación no negativo que es o bien retenido o bien descartado, nunca acotado superiormente por una cantidad que involucre $1/\varepsilon$.

Significancia y Reivindicaciones
El artículo afirma proporcionar una "base analítica estable y rigurosa" para la aproximación penalizada de flujos incompresibles. Al incorporar el amortiguamiento biharmónico y la corrección de Temam, el modelo logra resultados de bienestar mejorados en comparación con las ecuaciones de Navier–Stokes clásicas, específicamente en lo que respecta a la disponibilidad de estimaciones de energía $L^2$ incondicionales y suavizado de orden superior.

Los autores señalan que sus estimaciones uniformes sugieren que, a medida que $\varepsilon \to 0$, las soluciones fuertes convergen a una solución de un sistema de Navier–Stokes hiperviscoso incompresible. Sin embargo, declaran explícitamente que la identificación del sistema límite y la demostración de las tasas de convergencia en $\varepsilon$ son direcciones dejadas para trabajos futuros. Además, aunque la regularidad global para datos arbitrarios es conocida para el sistema de Navier–Stokes hiperviscoso incompresible (debido a Lions), los autores no reclaman tal resultado para el sistema penalizado con datos grandes, señalando que la estructura compresible asimétrica de Temam modifica las estimaciones de energía de orden superior de una manera que impide una extensión directa del resultado de datos grandes.

El trabajo se distingue de estudios previos sobre ecuaciones de Navier–Stokes amortiguadas (por ejemplo, con términos de absorción) al centrarse en un mecanismo de amortiguamiento de orden superior y un método de penalización que elimina el término de presión, simplificando así la formulación débil al eliminar la restricción de divergencia nula sobre las funciones de prueba.