Radiation magnetohydrodynamics modeling of an impulsively driven chromospheric jet in the solar atmosphere

Resumen Técnico: Modelado Magnetohidrodinámico con Radiación de un Chorro Cromosférico Impulsado de Forma Impulsiva

Enunciado del Problema
Los chorros solares, particularmente las espículas y macroespículas, son eyecciones dinámicas de plasma observadas en múltiples longitudes de onda. Aunque la reconexión magnética y las perturbaciones impulsivas (como pulsos de presión o velocidad) son ampliamente aceptadas como mecanismos desencadenantes, el papel de los procesos radiativos en la evolución de estos chorros sigue siendo objeto de investigación. Estudios numéricos anteriores han modelado chorros impulsados de forma impulsiva utilizando magnetohidrodinámica (MHD) ideal o no ideal, a menudo acoplada con funciones de enfriamiento radiativo empíricas. Sin embargo, existe una brecha en la simulación de estos chorros mediante ecuaciones de MHD no ideal acopladas completamente con un sistema de ecuaciones de transporte de radiación promediadas en frecuencia y ángulo. Específicamente, hay una falta de simulaciones que evolucionen independientemente tanto la densidad de energía radiativa como los componentes del flujo radiativo para determinar su impacto específico en la morfología y la disipación del chorro en la corona solar ópticamente delgada.

Metodología
El estudio emplea una simulación numérica bidimensional utilizando el código PLUTO (v4.4-patch3) para modelar un chorro cromosférico impulsado de forma impulsiva. El marco físico consiste en:

• Ecuaciones Gobernantes: El sistema resuelve ecuaciones de MHD no ideal acopladas con un esquema de hidrodinámica radiativa de dos momentos bajo la aproximación de gris. Esto incluye la evolución de la densidad del plasma, el momento, la energía total y el campo magnético, junto con la densidad de energía radiativa ($E_r$) y el flujo radiativo ($F_r$).
• Transporte de Radiación: Se utiliza el cierre M1 para cerrar las ecuaciones de momentos de radiación, permitiendo la descripción de ambos límites de difusión y de flujo libre. Los términos de interacción radiación-materia se derivan mediante transformaciones de Lorentz hasta el primer orden en $\beta = v/c$, teniendo en cuenta las opacidades de absorción y dispersión.
• Modelo de la Atmósfera Solar: El estado inicial es un equilibrio hidrostático basado en el modelo semiempírico C7, extendido a la corona. La composición del plasma asume un peso molecular medio de $\mu = 0.6$ (hidrógeno y helio completamente ionizados).
• Parámetros Físicos: La simulación utiliza una opacidad de absorción promedio de Planck baja ($\kappa = 1 \times 10^{-6}$ cm$^2$ g$^{-1}$) y una opacidad de dispersión promedio de Rosseland ($\sigma = 0.34$ cm$^2$ g$^{-1}$), características de la corona solar ópticamente delgada. La conducción térmica se modela con un enfoque limitado por flujo, y una función de calentamiento coronal empírica mantiene la atmósfera de fondo.
• Perturbación: La formación del chorro se inicia mediante un pulso gaussiano localizado en el componente de velocidad vertical ($v_y$) en $y = 1.75$ Mm con una amplitud de 100 km s$^{-1}$.
• Esquema Numérico: Las ecuaciones se resuelven utilizando un integrador temporal Runge-Kutta TVD de segundo orden. Las ecuaciones de MHD utilizan el solucionador de Riemann HLLD, mientras que el transporte de radiación utiliza el solucionador HLLC. Los esquemas Runge-Kutta Implícito-Explícito (IMEX) manejan los términos rígidos de interacción radiación-materia, y se emplea la aproximación de velocidad reducida de la luz (RSLA) con $\hat{c} = 0.1c$ para garantizar la viabilidad computacional.

Resultados Clave
La simulación reproduce con éxito la formación y evolución de una estructura de plasma colimada que se asemeja a una macroespícula:

• Morfología y Dinámica del Chorro: El pulso impulsivo genera una onda de choque que propulsa el plasma hacia arriba, formando un chorro que alcanza una altura de aproximadamente 8 Mm. El chorro exhibe una velocidad vertical máxima de $\sim$60 km s$^{-1}$ en el frente de choque y de $\sim$30 km s$^{-1}$ dentro del cuerpo del chorro. La estructura muestra una morfología de doble hilo en etapas posteriores, una característica comúnmente asociada con las espículas solares.
• Evolución Térmica: La cabeza del chorro alcanza temperaturas de $\sim$2 $\times$ 10$^5$ K, mientras que la mayor parte del chorro se estabiliza alrededor de 2 $\times$ 10$^5$ K. La región de choque alcanza inicialmente $\sim$9 $\times$ 10$^5$ K antes de enfriarse.
• Efectos Radiativos: La inclusión del transporte de radiación autoconsistente revela que la radiación actúa como un mecanismo de enfriamiento eficiente en la corona ópticamente delgada. La densidad de energía radiativa aumenta linealmente, mientras que el flujo radiativo disminuye, lo cual es consistente con las pérdidas radiativas ópticamente delgadas.
• Mecanismo de Disipación: Crucialmente, los términos de transporte de radiación contribuyen a la disipación del chorro. El acoplamiento entre el plasma y el campo de radiación induce velocidades de flujo negativas (movimiento descendente) en la corona, causando que el chorro se difunda y disipe de nuevo hacia la atmósfera inferior. Las variables radiativas no dictan la forma o morfología inicial del chorro, pero gobiernan significativamente su evolución y vida útil.
• Estabilidad del Fondo: La atmósfera solar de fondo permanece globalmente estable frente a los términos fuente (conducción térmica, calentamiento y radiación) en ausencia de la perturbación del chorro, observándose solo variaciones menores de temperatura.

Significado y Afirmaciones
El artículo afirma que este trabajo representa la primera simulación numérica de un chorro cromosférico impulsado de forma impulsiva donde la densidad de energía radiativa y el flujo se evolucionan de manera autoconsistente mediante ecuaciones de transporte de dos momentos, en lugar de depender de funciones de enfriamiento empíricas. El significado principal radica en demostrar que:

1. El Enfriamiento Radiativo es Crítico: Incluso en el régimen ópticamente delgado, donde las opacidades son bajas, el tratamiento explícito del transporte de radiación actúa como un mecanismo de enfriamiento dominante que limita la altura, la duración y el contenido energético de eventos similares a macroespículas.
2. Dinámica de Disipación: Los términos radiativos inducen flujos descendentes que son responsables del desvanecimiento rápido y la disipación del chorro, un proceso consistente con estudios anteriores que utilizan funciones de pérdida empíricas, pero ahora derivado de primeros principios dentro del marco del transporte de radiación.
3. Avance Metodológico: El estudio valida el uso del cierre M1 acoplado con MHD no ideal para modelar chorros solares, mostrando que, aunque la radiación no altera la morfología inicial del chorro, es esencial para capturar la evolución termodinámica correcta y las escalas de tiempo de disipación.

Los autores concluyen que, aunque el chorro simulado reproduce características dinámicas clave de las macroespículas observadas (como la formación de choques y estructuras de doble hilo), sus parámetros específicos (altura, velocidad, temperatura) caen dentro de los límites inferiores de los rangos observados, lo que sugiere que el evento modelado es un suceso de energía relativamente baja. El trabajo enfatiza que un tratamiento autoconsistente del transporte de radiación es necesario para modelar con precisión la disipación de tales fenómenos atmosféricos solares.