Phase-space averaging for stellar convection I. Liouvillian dynamics

Resumen Técnico: Promediado en el espacio de fases para la convección estelar I. Dinámica de Liouville

Planteamiento del Problema
La convección térmica es un mecanismo crítico que vincula la dinámica de fluidos local con la estructura global de las estrellas, gobernando la redistribución de energía, la mezcla química y la excitación de ondas. Sin embargo, un desafío persistente en la modelización de la evolución estelar es cerrar la brecha entre la hidrodinámica multidimensional y las ecuaciones de campo medio unidimensionales (1D) utilizadas en los códigos estándar. Actualmente, estos códigos dependen de la Teoría de la Longitud de Mezcla (MLT) y sus variantes, que tratan la convección de forma fenomenológica mediante una longitud de mezcla prescrita. Este enfoque carece de un vínculo estadístico o dinámico riguroso con las ecuaciones hidrodinámicas subyacentes, particularmente en regiones de transición complejas como las capas cercanas a la superficie, los límites convectivos y las zonas de sobreimpulso (overshooting). Las cantidades medias esenciales (por ejemplo, velocidades convectivas, entropías, flujos) a menudo carecen de una interpretación única, y los operadores de promediado utilizados son típicamente impuestos a priori en lugar de ser derivados de la dinámica del flujo.

Metodología
Los autores proponen un formalismo de Promediado en el Espacio de Fases (PSA, por sus siglas en inglés) que deriva descripciones de campo medio directamente de la moción hidrodinámica de los elementos de fluido convectivo. La metodología procede a través de los siguientes pasos:

1. Descripción Mesoscópica en el Espacio de Fases: El flujo se modela no como un continuo de puntos infinitesimales, sino como una distribución de "partículas de fluido mesoscópicas" en un espacio de posición-velocidad de seis dimensiones $(r, u)$. Estas partículas son lo suficientemente grandes como para estar en equilibrio termodinámico local (LTE), pero pequeñas en comparación con las escalas de la estructura estelar.
2. Ley de Conservación: Se define una distribución de densidad de masa $\rho(t, r, u)$ en este espacio de fases. Su evolución está gobernada por una ecuación de continuidad en el espacio $(r, u)$, que sirve como la ley de conservación fundamental. Esta ecuación es análoga a la ecuación de Liouville en la teoría cinética, pero tiene en cuenta la compresibilidad del flujo en el espacio de fases debido a procesos disipativos.
3. Derivación de las Ecuaciones de Campo Medio: Las cantidades de campo medio de Reynolds y Favre se definen como momentos del espacio de velocidades de la distribución $\rho$. Al integrar la ecuación de continuidad del espacio de fases sobre el espacio de velocidades, los autores derivan las ecuaciones estándar de conservación de masa, momento y energía de Reynolds-Favre. Crucialmente, se demuestra que la conmutación del operador de promediado con las derivadas espaciales y temporales se cumple exactamente dentro de este marco, a diferencia de las manipulaciones clásicas de campo medio donde suele ser una aproximación.
4. Formulación de Liouville: Bajo condiciones estándar de interior (alto número de Reynolds, bajo número de Mach y ajuste de presión inmediato), la dinámica de las partículas de fluido se reestructura en una forma de Liouville. Se introduce un generador $\psi$, análogo a un Hamiltoniano pero adaptado para medios disipativos. Esto permite que las ecuaciones de movimiento se escriban como ecuaciones generalizadas de Hamilton.
5. Diagnóstico de Estabilidad: Los autores analizan la divergencia del flujo en el espacio de fases, $\text{div}(\dot{\phi})$, la cual se identifica como el corchete de Poisson $\{s, T\}$ entre la entropía específica ($s$) y la temperatura ($T$).

Contribuciones Clave y Resultados

• Derivación Unificada: El artículo proporciona una derivación dinámicamente fundamentada de las ecuaciones de campo medio de Reynolds-Favre. El operador de promediado y la jerarquía de momentos emergen naturalmente como consecuencias de la ley de continuidad del espacio de fases, en lugar de ser ingredientes fenomenológicos separados.
• Generador de Liouville: Los autores demuestran que, incluso en medios estratificados y disipativos, la dinámica de las partículas de fluido puede ser generada por un potencial $\psi$ en el espacio de fases. Esta estructura separa claramente las contribuciones conservativas de los efectos térmicos no conservativos (producción de entropía).
• Divergencia del Espacio de Fases como Criterio de Estabilidad: Un resultado central es la identificación de la divergencia del espacio de fases $\{s, T\}$ como una medida lagrangiana local de la contracción o expansión de las trayectorias de partículas de fluido cercanas.
  • Límite Cuasi-Adiabático: En el límite donde la localidad y la cuasi-adiabaticidad se mantienen, el signo de $\{s, T\}$ se reduce exactamente al criterio de estabilidad de Schwarzschild clásico ($\text{sign}\{s, T\} = \text{sign}(-N^2)$, donde $N$ es la frecuencia de Brunt-Väisälä).
  • Caso General: Lejos del límite cuasi-adiabático (por ejemplo, en capas superficiales o regiones de penetración), el diagnóstico permanece resuelto en velocidad. Puede distinguir entre diferentes partes de la población convectiva dentro de la misma capa espacial, identificando regiones donde las trayectorias pueden estar convergiendo (estables) o divergiendo (inestables) basándose en sus velocidades específicas.
• Extensiones: Se muestra que el marco es extensible. Los apéndices esbozan cómo la rotación, los campos magnéticos y los cambios de composición pueden incorporarse mediante la modificación de la estructura del espacio de fases y los corchetes de Poisson asociados.

Significado y Reivindicaciones
El artículo afirma establecer una "ruta de Liouville" para la modelización de la convección estelar. Su principal significancia radica en proporcionar una base teórica donde el procedimiento de promediado, la jerarquía de momentos y el diagnóstico de estabilidad emanan de una única imagen de espacio de fases dinámicamente consistente.

Los autores enfatizan que este enfoque:

1. Elimina la necesidad de una longitud de mezcla prescrita externamente al fundamentar la descripción de campo medio en la acción real de los elementos de fluido.
2. Proporciona un diagnóstico de estabilidad local basado en trayectorias que es más matizado que el clásico criterio de Schwarzschild por capas, particularmente en regiones de transición donde este último puede fallar.
3. Sirve como el esqueleto geométrico y dinámico necesario para los modelos de cierre desarrollados en el artículo complementario (Artículo II), donde se construyen representaciones de máxima entropía de la distribución de fase $\rho$ para cerrar las ecuaciones de evolución estelar 1D.

El trabajo es modesto en cuanto a su alcance inmediato, señalando que la construcción de Liouville actual asume números de Mach bajos y un ajuste de presión inmediato. Excluye explícitamente los regímenes transónicos donde las fluctuaciones de presión dinámica rompen la correspondencia simple entre entropía y densidad, reconociendo que tales regímenes requieren hidrodinámica de radiación totalmente compresible o formalismos diferentes. Sin embargo, para los interiores profundos y las envolturas subsónicas de las estrellas, los autores argumentan que este marco ofrece un camino riguroso desde la hidrodinámica hacia las ecuaciones de la convección estelar de campo medio.