Expanding Solutions to Free Boundary 3D Spherically Symmetric Compressible Navier-Stokes-Poisson Equations near the Lane-Emden Stars

Resumen Técnico: Expansión de Soluciones para las Ecuaciones de Navier-Stokes-Poisson Compresibles con Frontera Libre en 3D con Simetría Esférica cerca de las Estrellas de Lane-Emden

Formulación del Problema
Este artículo investiga la existencia global y la dinámica a largo plazo de soluciones con simetría esférica de las ecuaciones de Navier-Stokes-Poisson (NSP) compresibles con una frontera libre en tres dimensiones. El sistema modela estrellas gaseosas politrópicas viscosas con una ecuación de estado $P(\rho) = K\rho^\gamma$, donde $\gamma \in (6/5, 4/3]$. Las ecuaciones describen la evolución de la densidad $\rho$ y la velocidad $u$ bajo la influencia de la autogravedad y el esfuerzo viscoso, con un dominio que se expande o contrae según una frontera libre $a(t)$ donde la densidad se desvanece continuamente hacia el vacío.

El estudio se centra en dos regímenes específicos:

1. Caso crítico de masa ($\gamma = 4/3$): Se requiere que la masa inicial sea estrictamente menor que la masa de la estrella de Lane-Emden (masa de Chandrasekhar).
2. Caso subcrítico ($\gamma \in (6/5, 4/3)$): Los datos iniciales deben pertenecer a un conjunto invariante específico $I$ definido por restricciones variacionales, el cual incluye datos cercanos a las estrellas de Lane-Emden pero permite configuraciones no cubiertas por la literatura previa.

Metodología
Los autores emplean una combinación de métodos variacionales, estimaciones de energía a priori y argumentos de compacidad para construir soluciones débiles globales.

1. Control Variacional de la Energía Gravitatoria:
   Un desafío central en el rango $\gamma \in (6/5, 4/3]$ es controlar la energía gravitatoria negativa, que puede dominar la energía interna. Basándose en los métodos variacionales establecidos en [7], los autores utilizan las propiedades de las soluciones de Lane-Emden para derivar cotas inferiores para el funcional $S_\mu(\rho)$.

   • Para $\gamma = 4/3$, la condición $M(\rho_0) < M_{ch}$ asegura que la energía gravitatoria esté controlada por la energía interna.
   • Para $\gamma \in (6/5, 4/3)$, los autores definen un conjunto invariante $I$ donde el funcional $Q(\rho)$ permanece positivo. Demuestran que este conjunto incluye datos iniciales cerca de las estrellas de Lane-Emden que fueron excluidos previamente en trabajos como [25], específicamente aquellos con pequeña energía cinética pero perfiles de densidad específicos.

2. Aproximación y Regularización:
   Para manejar la frontera del vacío y la singularidad en el origen, los autores consideran un problema aproximado en un dominio $[\xi, a(t)]$ con $\xi > 0$. Establecen la existencia de soluciones fuertes locales para este problema aproximado y luego las extienden globalmente.

3. Estimaciones A Priori:
   La prueba se basa en la derivación de estimaciones uniformes independientes del parámetro de regularización $\xi$:

   • Estimación de Energía Básica: Establece la conservación de la energía total menos la disipación debida a la viscosidad.
   • Estimaciones de Orden Superior: Los autores derivan una mayor integrabilidad espacio-temporal para la densidad ($\rho^{2\gamma}$) y estimaciones de $L^2$ ponderadas para las derivadas de la velocidad y la densidad en coordenadas lagrangianas.
   • Separación de Frontera y Origen: El dominio espacial se divide en una región cerca del origen y una región cerca de la frontera libre.
     • Cerca de la frontera libre, la mayor regularidad permite una convergencia fuerte mediante el teorema de Arzela-Ascoli.
     • Cerca del origen, donde la regularidad es menor, los autores adaptan técnicas de [22] (originalmente para Navier-Stokes en $\mathbb{R}^n$) para probar la convergencia fuerte de la densidad. Esto implica probar que el límite débil del término de presión $\rho^\gamma$ es igual a la presión del límite débil, utilizando la convexidad de la función de presión e desigualdades integrales específicas.

4. Compacidad y Paso al Límite:
   Al combinar la convergencia fuerte cerca de la frontera y del origen, los autores extraen una subsucesión que converge a una solución débil global que satisface la formulación distributiva de las ecuaciones NSP.

Resultados Clave

1. Existencia Global de Soluciones Débiles (Teoremas 2.1 y 2.2):
   El artículo demuestra la existencia de soluciones débiles globales para $\gamma \in (6/5, 4/3]$ bajo las condiciones especificadas de masa o del conjunto invariante.

   • Teorema 2.1: Establece la existencia de una solución débil global $(\rho, u, a)$ con densidad no negativa y energía acotada.
   • Teorema 2.2: Proporciona resultados de mayor regularidad. Específicamente, establece la regularidad interior ($H^1$ para la densidad, $H^2$ para la velocidad) lejos del origen y de la frontera, así como una regularidad específica cerca de la frontera libre, incluyendo cotas sobre el flujo viscoso efectivo $\sigma$.

2. Expansión de la Frontera Libre (Teorema 2.3):
   Asumiendo la existencia de soluciones fuertes, el artículo demuestra que el soporte de la solución se expande al infinito a una tasa algebraica. Esto indica la fuerte inestabilidad de las soluciones de Lane-Emden en este régimen.

   • Para $\gamma \in (6/5, 4/3)$, el radio $a(t)$ satisface $a(t) \approx (t / (\eta + 4/3\varepsilon))^{1/3}$.
   • Para $\gamma = 4/3$, la tasa de expansión está acotada por $(1+t)^{1/4} \lesssim a(t) \lesssim (t / (\eta + 4/3\varepsilon))^{1/3}$.
   • La prueba utiliza un argumento de tipo Virial involucrando el funcional $H(t) = \int \rho r^4 dr$ y sus derivadas temporales.

3. Viscosidad Dependiente de la Densidad (Observación 2.3):
   Los autores extienden el análisis de la tasa de expansión a casos con viscosidad dependiente de la densidad $\varepsilon(\rho), \eta(\rho) \propto \rho^\alpha$. Muestran que la tasa de expansión depende del exponente $\alpha$:

   • Si $0 \le \alpha < 2/3$, la tasa es $a(t) \gtrsim (t/\eta)^{1/(3(1-\alpha))}$.
   • Si $2/3 \le \alpha \le \gamma$, la tasa es lineal, $a(t) \gtrsim t$.
   • Esto resalta que el desvanecimiento de la viscosidad volumétrica altera significativamente la dinámica de expansión cuando $\alpha$ es pequeño.

Significado y Reivindicaciones
El artículo afirma resolver la existencia global de soluciones débiles para las ecuaciones NSP de frontera libre en los rangos de exponente adiabático crítico y subcrítico ($\gamma \in (6/5, 4/3]$) con viscosidad constante. Una contribución primaria es el refinamiento del conjunto de datos iniciales para $\gamma \in (6/5, 4/3)$; los autores demuestran que su conjunto invariante $I$ es más amplio que las restricciones previas (por ejemplo, [25]), permitiendo datos iniciales cerca de las estrellas de Lane-Emden que anteriormente se consideraban inalcanzables.

Además, el artículo establece la expansión algebraica de la frontera libre, proporcionando una confirmación matemática rigurosa de la inestabilidad de las estrellas de Lane-Emden en este régimen. Los autores señalan que, si bien la tasa de expansión se demuestra para soluciones fuertes, la existencia global de tales soluciones fuertes sigue siendo un problema abierto, lo que hace que el resultado de inestabilidad sea condicional. El trabajo cierra la brecha entre los resultados de estabilidad para $\gamma > 4/3$ y los resultados de inestabilidad para $\gamma < 4/3$ en el contexto de estrellas gaseosas viscosas.