Stability of purely convective steady-states of fractional Boussinesq equations in an exterior domain

Resumen Técnico: Estabilidad de Estados Estacionarios Puramente Convectivos de Ecuaciones de Boussinesq Fraccionarias en un Dominio Exterior

Formulación del Problema
El artículo investiga la estabilidad global y el comportamiento asintótico de flujos de convección térmica en un dominio no acotado tridimensional $\Omega = \{x \in \mathbb{R}^3; |x| > 1\}$, exterior a una esfera unitaria. El fluido se calienta desde la frontera de la esfera, creando una fuerza de flotación impulsada por la gravedad. El estudio modifica el sistema clásico de Boussinesq reemplazando la fuerza viscosa estándar (determinada por el operador de Stokes $-\Delta$) con una potencia fraccionaria del operador de Stokes, denotada como $A^\kappa$ (donde $0 < \kappa \leq 1$).

El sistema rector para las variables de perturbación $(u, \theta)$, que representan las desviaciones de velocidad y temperatura de un estado estacionario puramente conductivo $(u_0, T_0)$, viene dado por:
$$\begin{cases} \partial_t u + A^\kappa u + P(\alpha g \theta) = -P(u \cdot \nabla u), \\ \partial_t \theta - \Delta \theta + u \cdot \nabla T_0 = -u \cdot \nabla \theta, \\ u|_{\partial \Omega} = 0, \quad \theta|_{\partial \Omega} = 0, \end{cases}$$
donde $P$ es la proyección de Helmholtz, $g$ es el vector de aceleración gravitacional y $\alpha$ es una constante relacionada con la temperatura de la frontera. El objetivo principal es establecer la estabilidad global en $L^2$ del estado estacionario y derivar tasas de decaimiento algebraico para las perturbaciones cuando $t \to \infty$.

Metodología
El análisis procede a través de cuatro etapas principales:

1. Teoría de Operadores y Generación de Semigrupos: Los autores definen un operador lineal $L$ asociado con el sistema linealizado. Utilizando el teorema de Kato–Rellich, demuestran que $L$ es autoadjunto y positivo bajo la condición $0 < \alpha < 2^{-\kappa}$. Esto permite la construcción de un semigrupo analítico fuertemente continuo $e^{-tL}$, que gobierna la evolución lineal del sistema.
2. Existencia de Soluciones Débiles: Para manejar los términos no lineales, los autores emplean una técnica de regularización mediante operadores $J_n$ e $I_n$. Construyen soluciones aproximadas a la forma integral de la ecuación del sistema. Al establecer estimaciones de energía uniformes y utilizar argumentos de compacidad (específicamente resultados del tipo lema de Aubin-Lions), demuestran la existencia de una solución débil global en el espacio $L^\infty(0, \tau; L^2) \cap L^2(0, \tau; D(A^{\kappa/2}) \times W^{1,2}_0)$.
3. Estimaciones de Decaimiento en $L^2$: El decaimiento de la solución se analiza dividiendo el término integral no lineal en la formulación de la solución suave. Los autores utilizan la representación espectral del operador $L$ y derivan tasas de decaimiento específicas para el semigrupo $e^{-tL}$ actuando sobre los datos iniciales. Crucialmente, emplean desigualdades de interpolación (Lema 2.2 y Lema 2.3) y embebimientos de Sobolev para acotar los términos convectivos no lineales frente a las tasas de decaimiento lineales.
4. Iteración de Bootstrap: Para establecer la tasa de decaimiento algebraico óptima, el artículo utiliza un esquema de iteración de bootstrap. Comenzando desde una tasa de decaimiento base, los autores mejoran iterativamente la estimación alimentando el límite de decaimiento anterior de nuevo en las estimaciones integrales no lineales, convergiendo eventualmente a la tasa de decaimiento objetivo determinada por el problema linealizado.

Contribuciones y Resultados Clave
El artículo establece los siguientes resultados principales bajo las condiciones $3/4 < \kappa \leq 1$, $0 < \beta < 1/2$ y $0 < \alpha < 2^{-\kappa}$:

• Existencia Global: Para cualquier dato inicial $(a, b) \in L^2_\sigma(\Omega)^3 \times L^2(\Omega)$, el sistema de Boussinesq fraccionario admite una solución débil global.
• Estabilidad Global: La solución $(u(t), \theta(t))$ converge a cero en la norma $L^2$ cuando $t \to \infty$, confirmando la estabilidad global del estado estacionario puramente conductivo.
• Decaimiento Algebraico: Si el sistema linealizado exhibe una tasa de decaimiento algebraico de $t^{-\beta}$, el sistema no lineal hereda esta tasa. Específicamente, la solución satisface:
  $$\|(u(t), \theta(t))\|_2 \leq C t^{-\beta}, \quad t > 0.$$
  La tasa de decaimiento está restringida por la potencia fraccionaria $\kappa$ y las propiedades de los datos iniciales.

Significado y Alcance
El artículo se posiciona como el primer intento de comprender la estabilidad de flujos de Boussinesq exteriores fraccionarios. Extiende trabajos previos sobre las ecuaciones de Navier-Stokes clásicas (donde $\kappa=1$) y las ecuaciones estándar de Boussinesq al caso fraccionario en un dominio exterior.

Los autores señalan explícitamente que el resultado de estabilidad depende de la suposición de que el parámetro de fuerza de flotación $\alpha$ es suficientemente pequeño. Reconocen que si la diferencia de temperatura es grande, pueden surgir inestabilidades (como las observadas en la convección de Rayleigh-Bénard), haciendo inviable la estabilidad. El estudio no propone nuevas aplicaciones o configuraciones experimentales, sino que proporciona una base matemática rigurosa para el comportamiento asintótico de estos modelos específicos de fluidos fraccionarios, llenando un vacío en la literatura respecto a flujos exteriores fraccionarios.