Subvenciones relacionadas:
General
La simulación numérica mediante códigos complejos de ordenador es una herramienta fundamental en la investigación física y en la técnica desde hace décadas. El crecimiento vertiginoso de las capacidades informáticas junto con el avance notable de la matemática numérica ha hecho accesible a los centros de investigación de tamaño medio esta rama de la investigación, a caballo entre la física teórica y la física experimental. La astrofísica no es excepción a lo anterior, habiéndose desarrollado desde finales de los 70 una especialidad de la misma, la astrofísica computacional, que ha permitido llegar a comprender gran variedad de fenómenos inaccesibles a la investigación teórica pura y dar cuenta de observaciones hasta entonces inexplicadas. Su mayor campo de aplicación en las décadas pasadas han sido los fenómenos (magneto) hidrodinámicos y de dinámica de gases en multiplicidad de entornos cósmicos, por ejemplo los interiores y atmósferas estelares y planetarios y el medio interestelar, incluyendo magnetoconvección y dínamo, discos de acreción, evolución de nebulosas planetarias, explosiones y restos de supernova, etc. La incorporación a las simulaciones numéricas de las ecuaciones del transporte radiativo, ocurrida ya en décadas pasadas, ha permitido dotar de mayor realismo a los estudios de procesos hidrodinámicos en fotosferas y cromosferas estelares.
El presente Proyecto quiere apoyar el desarrollo en el IAC de la investigación astrofísica basada en el uso de grandes códigos numéricos que requieren el uso de ordenadores masivamente paralelos y su enlace con los resultados de observación. Objetivo general de este Proyecto es la realización de cálculos de física de fluidos cósmicos y de transporte radiativo. La temática de dichos cálculos se centrará en
- fenómenos de dinámica de gases magnetizados en interiores y atmósferas estelares
- transporte de radiación y señales de polarización en líneas espectrales en base a modelos atómicos y moleculares realistas y los efectos Hanle y Zeeman
- comparación de resultados teórico/numéricos con datos de observación
Este Proyecto es especialmente relevante a la vista de la involucración, cada vez mayor, del IAC en las redes de supercomputación nacionales y europeas y, en general, en grandes iniciativas de instalación de superordenadores.
Miembros
Resultados
A continuación, destacamos los resultados de nuestro resumen anual de 2022.
A lo largo del año 2022, los efectos de ionización parcial, los efectos de ionización fuera de equilibrio y los fluidos múltiples han sido uno de los principales bloques de desarrollo tanto desde la perspectiva teórica como numérica. Por ejemplo, se ha logrado una generalización de las ecuaciones de Braginskii de 1965 para plasmas multi-especies generales con masas y temperaturas arbitrarias, donde todas las viscosidades y flujos de calor en el modelo se describen mediante sus propias ecuaciones de evolución. Este nuevo enfoque tiene una ventaja crucial en la que las componentes paralelas a lo largo de las líneas de campo magnético no se vuelven ilimitadas (infinitamente grandes) en regímenes de baja colisionalidad de interés para este grupo, como, por ejemplo, la corona solar (Hunana et al. 2022). En este bloque temático, también se han realizado simulaciones 2D y 3D utilizando un modelo de dos fluidos que trata las especies neutras e ionizadas como dos componentes separados, para analizar el efecto que tiene la interacción colisional entre ambas componentes en la dinámica de la lluvia coronal, la evolución de la inestabilidad de Kelvin-Helmholtz, la propagación de ondas magnetoacústicas a través de la cromosfera solar o el calentamiento del plasma (Martínez-Gómez et al. 2022a). Otro ejemplo de desarrollo teórico con posibles aplicaciones numéricas ha sido la búsqueda de los efectos de la difusión ambipolar en la cromosfera desde una perspectiva más fundamental mediante soluciones analíticas. Las soluciones obtenidas para casos con simetría cilíndrica se demuestran como una prueba exigente, pero no obstante viable, para los códigos magnetohidrodinámicos (MHD) que incorporan difusión ambipolar. Además, se han realizado ejecuciones detalladas de tablas de las soluciones disponibles públicamente para la comunidad (Moreno-Insertis et al. 2022). Por último, se han comenzado a estudiar los efectos de la ionización fuera de equilibrio del átomo de hidrógeno junto con el estudio de los efectos Lyman α en configuraciones simples para aplicarlos más tarde en simulaciones realistas que incluyan la cromosfera.
La mejora y prueba de las capacidades de los códigos MHD disponibles en el grupo solar ha sido otro de los principales desarrollos clave realizados en 2022. Por ejemplo, los resultados obtenidos por Moreno-Insertis et al. 2022 se utilizaron para verificar que el código MHD Bifrost es capaz de reproducir las soluciones teóricas con la suficiente precisión hasta tiempos de difusión muy avanzados, así como para explorar las propiedades asintóticas de estas soluciones teóricas. Además de eso, se han realizado varios cambios en el código MANCHA, cuyo objetivo era aumentar la eficiencia y agregar nuevas características que permitieran a los investigadores realizar experimentos más realistas y explorar nuevas áreas de investigación. Por ejemplo, el código MANCHA se ha extendido para poder simular simulaciones solares hasta la corona, agregando un nuevo módulo que calcula de manera eficiente uno de los ingredientes clave en la corona: la conducción térmica (Navarro et al. 2022). La preparación del código MANCHA para su extensión multifluida con radiación también ha sido otra rama de trabajo relacionada con el desarrollo numérico en 2022. Además, se han desarrollado nuevas rutinas de ecuaciones de estado y opacidad que permiten separar las contribuciones de fondo en equilibrio de las tratadas fuera del equilibrio. Además de enfrentar diferentes desafíos en la física solar, el gran desarrollo generado en MANCHA es útil para estudiar estrellas frías de secuencia principal (G, K, M), lo que contribuye a una mejor comprensión de la física estelar. Para llevar a cabo todas estas tareas, fue necesario no solo realizar numerosas pruebas de escalado y experimentos numéricos en máquinas locales en el IAC, sino también en supercomputadoras como LaPalma, PICASSO, PizDaint y MareNostrum4; así como trabajar en colaboración con colaboradores externos.
Durante 2022, en este proyecto también se ha centrado en diferentes fenómenos de la atmósfera solar y la correspondiente comparación con observaciones. Como ejemplo ilustrativo, se han modelado por primera vez Puntos Brillantes Coronales (CBPs) con la suficiente realismo para desentrañar los mecanismos que los generan y proporcionarles energía, siendo capaces también de explicar diferentes características observadas desde satélites espaciales. La comparación con observaciones se realiza a través de imágenes sintéticas de SDO/AIA, Solar Orbiter EUI-HRI e IRIS que se han calculado a partir del experimento numérico realizado con el código Bifrost (Nóbrega-Siverio y Moreno-Insertis, 2022). Otro ejemplo es la combinación de experimentos numéricos en 3D con el código MoLMH y modelado directo utilizando la línea Hα para estudiar oscilaciones transversales de hilo prominencial. Los resultados contienen implicaciones relevantes para el campo de la sismología de prominencias, mostrando que la emisión Hα se puede utilizar para detectar el modo fundamental de las oscilaciones (Martínez-Gómez et al. 2022b). Además, se han analizado observaciones de alta resolución en tierra de fenómenos ejectivos como chorros en la atmósfera solar, encontrando similitudes sorprendentes con resultados obtenidos de experimentos numéricos. Además, ha habido contribuciones significativas de los miembros de este proyecto al avance de las observaciones y la construcción de nuevos telescopios (Quintero et al. 2022) y satélites (De Pontieu et al. 2022, Cheung et al. 2022), utilizando el conocimiento adquirido de los experimentos teórico-numéricos. Finalmente, se realizó un primer intento exploratorio para comprender la física de agujeros coronales y regiones activas desde un punto de vista global a través de soluciones magnetohidrostáticas en 2D (Terradas et al. 2022), lo que requerirá un mayor desarrollo en los próximos años para su comparación con observaciones.
Por último, pero no menos importante, se han aplicado herramientas de vanguardia como las proporcionadas por el Aprendizaje Automático (Machine Learning) y la estadística Bayesiana a problemas de la atmósfera solar. En este sentido, se lanzó un proyecto para caracterizar los límites de los métodos de k-means y su aplicación a observaciones solares. Además, se han iniciado nuevos desarrollos en códigos de transferencia radiativa para utilizarlos en un estudio preliminar de un enfoque de aprendizaje automático para el cálculo de términos radiativos. El desarrollo de la aplicación de técnicas bayesianas a la comparación de modelos en la sismología de la atmósfera solar continuó en 2022, con la publicación de un artículo de revisión que recoge los principales resultados obtenidos en la última década (Arregui 2022a). Además, el formalismo bayesiano se ha aplicado con éxito a la predicción de la amplitud del ciclo de actividad solar, proponiendo una nueva metodología para cuantificar la bondad tanto de la predicción como del modelo subyacente (Arregui 2022b).
Actividad científica
Publicaciones relacionadas
-
Probing the Physics of the Solar Atmosphere with the Multi-slit Solar Explorer (MUSE). I. Coronal HeatingThe Multi-slit Solar Explorer (MUSE) is a proposed mission composed of a multislit extreme ultraviolet (EUV) spectrograph (in three spectral bands around 171 Å, 284 Å, and 108 Å) and an EUV context imager (in two passbands around 195 Å and 304 Å). MUSE will provide unprecedented spectral and imaging diagnostics of the solar corona at high spatial (De Pontieu, Bart et al.
Fecha de publicación:
22022 -
Probing the Physics of the Solar Atmosphere with the Multi-slit Solar Explorer (MUSE). II. Flares and EruptionsCurrent state-of-the-art spectrographs cannot resolve the fundamental spatial (subarcseconds) and temporal (less than a few tens of seconds) scales of the coronal dynamics of solar flares and eruptive phenomena. The highest-resolution coronal data to date are based on imaging, which is blind to many of the processes that drive coronal energeticsCheung, Mark C. M. et al.
Fecha de publicación:
22022 -
Transverse kink oscillations of inhomogeneous prominence threads: Numerical analysis and Hα forward modellingContext. Prominence threads are very long and thin flux tubes that are partially filled with cold plasma. Observations have shown that transverse oscillations are frequent in these solar structures. The observations are usually interpreted as the fundamental kink mode, while the detection of the first harmonic remains elusive. Aims: The propertiesMartínez-Gómez, David et al.
Fecha de publicación:
22022 -
-
Solar surges related to UV bursts. Characterization through k-means, inversions, and density diagnosticsContext. Surges are cool and dense ejections typically observed in chromospheric lines and closely related to other solar phenomena such as UV bursts or coronal jets. Even though surges have been observed for decades now, questions regarding their fundamental physical properties such as temperature and density, as well as their impact on upperNóbrega-Siverio, D. et al.
Fecha de publicación:
112021 -
Large-amplitude longitudinal oscillations in solar prominences simulated with different resolutionsContext. Large-amplitude longitudinal oscillations (LALOs) in solar prominences have been widely studied in recent decades. However, their damping and amplification mechanisms are not well understood. Aims: In this study, we investigate the attenuation and amplification of LALOs using high-resolution numerical simulations with progressivelyLiakh, V. et al.
Fecha de publicación:
102021 -
Effect of momentum and heat losses on the hydrodynamic instability of a premixed equidiffusive flame in a Hele-Shaw cellThe linear stage of hydrodynamic instability of a laminar premixed flame propagating in a Hele-Shaw cell is investigated. Our theoretical model takes into account momentum and heat losses, temperature-dependent transport coefficients, and the continuous internal structure of the flame front. The dispersion relation is obtained numerically as aHan, Yifan et al.
Fecha de publicación:
102021 -
Evidence of the multi-thermal nature of spicular downflows. Impact on solar atmospheric heatingContext. Spectroscopic observations of the emission lines formed in the solar transition region commonly show persistent downflows on the order of 10−15 km s −1. The cause of such downflows, however, is still not fully clear and has remained a matter of debate. Aims: We aim to understand the cause of such downflows by studying the coronal andBose, Souvik et al.
Fecha de publicación:
102021 -
Newly formed downflow lanes in exploding granules in the solar photosphereContext. Exploding granules have drawn renewed interest because of their interaction with the magnetic field (either emerging or already present). Especially the newly forming downflow lanes developing in their centre seem to be eligible candidates for the intensification of magnetic fields. We analyse spectroscopic data from two differentEllwarth, M. et al.
Fecha de publicación:
92021 -
Two-fluid simulations of Rayleigh-Taylor instability in a magnetized solar prominence thread. II. Effects of collisionalitySolar prominences are formed by partially ionized plasma with inter-particle collision frequencies generally warranting magnetohydrodynamic treatment. In this work we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluidsPopescu Braileanu, B. et al.
Fecha de publicación:
62021 -
Bayesian Evidence for a Nonlinear Damping Model for Coronal Loop OscillationsRecent observational and theoretical studies indicate that the damping of solar coronal loop oscillations depends on the oscillation amplitude. We consider two mechanisms: linear resonant absorption and a nonlinear damping model. We confront theoretical predictions from these models with observed data in the plane of observables defined by theArregui, Iñigo
Fecha de publicación:
72021 -
Chromospheric Heating by Magnetohydrodynamic Waves and InstabilitiesThe importance of the chromosphere in the mass and energy transport within the solar atmosphere is now widely recognized. This review discusses the physics of magnetohydrodynamic waves and instabilities in large-scale chromospheric structures as well as in magnetic flux tubes. We highlight a number of key observational aspects that have helped ourSrivastava, A. K. et al.
Fecha de publicación:
62021 -
Simulations of the Biermann battery mechanism in two-fluid partially ionised plasmasContext. In the absence of an initial seed, the Biermann battery term of a non-ideal induction equation acts as a source that generates weak magnetic fields. These fields are then amplified via a dynamo mechanism. The Kelvin-Helmholtz instability is a fluid phenomenon that takes place in many astrophysical scenarios and can trigger the action ofMartínez-Gómez, D. et al.
Fecha de publicación:
62021 -
A New View of the Solar Interface Region from the Interface Region Imaging Spectrograph (IRIS)The Interface Region Imaging Spectrograph (IRIS) has been obtaining near- and far-ultraviolet images and spectra of the solar atmosphere since July 2013. IRIS is the highest resolution observatory to provide seamless coverage of spectra and images from the photosphere into the low corona. The unique combination of near- and far-ultraviolet spectraDe Pontieu, Bart et al.
Fecha de publicación:
52021 -
Large-amplitude Prominence Oscillations following Impact by a Coronal JetObservational evidence shows that coronal jets can hit prominences and set them in motion. The impact leads to large-amplitude oscillations (LAOs) in the prominence. In this paper, we attempt to understand this process via 2.5D MHD numerical experiments. In our model, the jets are generated in a sheared magnetic arcade above a parasitic bipolarLuna, Manuel et al.
Fecha de publicación:
52021 -
Critical Science Plan for the Daniel K. Inouye Solar Telescope (DKIST)The National Science Foundation's Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand, and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinaryRast, Mark P. et al.
Fecha de publicación:
42021 -
Acoustic-gravity wave propagation characteristics in three-dimensional radiation hydrodynamic simulations of the solar atmosphereThere has been tremendous progress in the degree of realism of three-dimensional radiation magneto-hydrodynamic simulations of the solar atmosphere in the past decades. Four of the most frequently used numerical codes are Bifrost, CO5BOLD, MANCHA3D and MURaM. Here we test and compare the wave propagation characteristics in model runs from theseFleck, B. et al.
Fecha de publicación:
22021 -
Influence of ambipolar and Hall effects on vorticity in three-dimensional simulations of magneto-convectionThis paper presents the results of the analysis of three-dimensional simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation runs are compared: one including both ambipolar diffusion and the Hall effect; one including only ambipolar diffusion and one without any of theseKhomenko, E. et al.
Fecha de publicación:
22021 -
The chromospheric component of coronal bright points. Coronal and chromospheric responses to magnetic-flux emergenceContext. We investigate the chromospheric counterpart of small-scale coronal loops constituting a coronal bright point (CBP) and its response to a photospheric magnetic-flux increase accompanied by co-temporal CBP heating. Aims: The aim of this study is to simultaneously investigate the chromospheric and coronal layers associated with a CBP, and inMadjarska, Maria S. et al.
Fecha de publicación:
22021 -
Stability Analysis for an Interface with a Continuous Internal StructureA general method for solving a linear stability problem of an interface with a continuous internal structure is described. Such interfaces or fronts are commonly found in various branches of physics, such as combustion and plasma physics. It extends simplified analysis of an infinitely thin discontinuous front by means of numerical integrationModestov, Mikhail
Fecha de publicación:
12021