On a Local Existence Theorem for the Evolution Equation of Viscous Gaseous Stars in a Physical Vacuum

Resumen Técnico: Teorema de Existencia Local para Estrellas Gaseosas Viscosas en un Vacío Físico

Formulación del Problema
Este artículo aborda la existencia local de una solución bien planteada para el problema de frontera libre (VFBP) de las ecuaciones de Navier–Stokes–Poisson compresibles (CNSP) en tres dimensiones que describen estrellas gaseosas autogravitantes viscosas. El sistema se caracteriza por viscosidades degeneradas, donde los coeficientes de viscosidad de corte y de volumen dependen de la densidad mediante leyes de potencia, $\mu(\rho) = a_1\rho^\delta$ y $\lambda(\rho) = a_2\rho^\delta$, con $\delta \in (0, 1)$.

El entorno físico involucra una frontera de "vacío físico", definida por la condición de que la velocidad del sonido al cuadrado, $c^2$, tiene una derivada normal negativa en la frontera ($-\infty < \partial_n c_0^2 < 0$). Esta condición implica que la densidad inicial $\rho_0$ se anula en la frontera como una función de distancia ($\rho_0^{\gamma-1} \sim \text{dist}(x, \partial\Omega)$). Los autores se centran en el movimiento barotrópico con simetría esférica, reduciendo el problema de 3-D a un problema de valor inicial-frontera (IBVP) de 1-D en un intervalo móvil $I(t) = [0, R(t))$. El principal desafío matemático radica en la fuerte degeneración de las ecuaciones cerca de la frontera del vacío y en la singularidad de coordenadas en el origen, lo que complica la propagación de la regularidad del campo de velocidades.

Metodología
Los autores emplean una reformulación lagrangiana del problema para fijar el dominio espacial a $I = [0, 1)$. La metodología central consiste en un marco analítico de múltiples pasos:

1. Linealización y Aproximación: Se linealiza el problema no lineal alrededor de un mapa de flujo $\bar{\eta}$ dado. Para resolver el problema linealizado resultante, los autores utilizan un esquema de Galerkin modificado. Construyen una base de Hilbert $\{\xi_j\}$ derivada de un problema de Sturm–Liouville específico, adaptado a la simetría esférica y a las condiciones de contorno de Neumann. Esto permite la construcción de soluciones de aproximación globales en el tiempo.
2. Estimaciones de Energía Uniformes: Una parte significativa del trabajo se dedica a establecer estimaciones de energía uniformes para las soluciones lineales. Los autores definen un funcional de energía total $E(t)$ y un funcional de disipación $D(t)$, los cuales se dividen en componentes interiores y exteriores.
   • Estimaciones Interiores: Cerca del origen, se utilizan normas de Sobolev estándar, aprovechando la suavidad de la solución en el interior.
   • Estimaciones Exteriores: Cerca de la frontera del vacío, los autores introducen espacios de Sobolev ponderados con pesos que involucran potencias de la densidad inicial $\rho_0$. Estos pesos se eligen cuidadosamente para manejar la degeneración y asegurar que la solución permanezca suave hasta la frontera móvil. Los pesos dependen del exponente adiabático $\gamma$ y del parámetro de viscosidad $\delta$.
3. Bootstrap y Regularidad: A través de un argumento de bootstrapping que involucra estimaciones elípticas e inserciones de derivadas cruzadas (detalladas en el Apéndice D), los autores derivan una regularidad de orden superior para la velocidad y el mapa de flujo. Abordan específicamente los términos de "derivada cruzada" que surgen en la ecuación de cantidad de movimiento, demostrando que las derivadas de orden superior pueden controlarse en los espacios ponderados.
4. Iteración de Picard: Finalmente, la existencia local de soluciones para el problema no lineal se establece mediante una iteración de Picard clásica. Las estimaciones uniformes derivadas para el problema linealizado se utilizan para demostrar la convergencia fuerte de la secuencia iterativa hacia una solución clásica única.

Contribuciones Clave y Resultados
El artículo establece los siguientes resultados principales:

• Existencia Local (Well-Posedness): Los autores demuestran la existencia y unicidad de soluciones clásicas para las ecuaciones CNSP con simetría esférica, viscosidades degeneradas y fronteras de vacío físico. Esto es válido para todos los exponentes adiabáticos $\gamma > 4/3$ y parámetros de viscosidad $\delta \in (0, 1)$.
• Regularidad hasta la Frontera: A diferencia de resultados previos para flujos viscosos que a menudo producían soluciones fuertes que no eran suaves hasta la frontera móvil, este trabajo demuestra que las soluciones obtenidas son suaves en todo el camino hasta la frontera de vacío móvil para $t > 0$. Las soluciones capturan el comportamiento específico del vacío físico de la configuración de la estrella de Lane-Emden.
• Datos Iniciales Generales: Los resultados se mantienen sin supuestos de pequeñez sobre los datos iniciales, siempre que la densidad inicial cumpla la condición de vacío físico ($\rho_0^{\gamma-1} \sim 1-r$) y pertenezca a los espacios de Sobolev ponderados apropiados ($H^3$ para la densidad, $H^6$ ponderado para la velocidad).
• Condiciones de Contorno: El artículo deriva rigurosamente las condiciones de contorno que satisface la solución, mostrando que la condición de ausencia de tensión (stress-free) se cumple automáticamente y que la velocidad satisface una relación específica en la frontera del vacío que involucra los coefores de viscosidad.

Significancia
La importancia de este trabajo, tal como la presentan los autores, reside en superar las dificultades intrínsecas de la fuerte degeneración de las ecuaciones cerca del vacío físico. Estudios previos sobre estrellas gaseosas viscosas a menudo resultaban en soluciones que carecían de suavidad en la frontera o requerían supuestos restrictivos sobre los datos iniciales (por ejemplo, tasas de decaimiento de la densidad más altas cerca del vacío). Al desarrollar un marco de energía a medida con normas ponderadas específicas y utilizar un esquema de Galerkin modificado, este artículo extiende la teoría de existencia local al caso físicamente relevante de las estrellas de Lane-Emden con viscosidades degeneradas. Los resultados cierran la brecha entre el caso de Euler-Poisson no viscoso, bien comprendido, y el caso viscoso más complejo, proporcionando una base rigurosa para el estudio de la evolución de las estrellas gaseosas viscosas en un vacío físico.