Simulación Numérica de Procesos Astrofísicos

Año de inicio
2003
Unidad organizativa

Subvenciones relacionadas:

    General
    Descripción

    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.

    Investigador principal

    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).

    Publicaciones relacionadas

    • Deciphering Solar Coronal Heating: Energizing Small-scale Loops through Surface Convection
      The solar atmosphere is filled with clusters of hot small-scale loops commonly known as coronal bright points (CBPs). These ubiquitous structures stand out in the Sun by their strong X-ray and/or extreme-ultraviolet (EUV) emission for hours to days, which makes them a crucial piece when solving the solar coronal heating puzzle. In addition, they
      Nóbrega-Siverio, D. et al.

      Fecha de publicación:

      12
      2023
      Número de citas
      6
    • Bayesian evidence for two slow-wave damping models in hot coronal loops
      We computed the evidence in favour of two models, one based on field-aligned thermal conduction alone and another that includes thermal misbalance as well, to explain the damping of slow magneto-acoustic waves in hot coronal loops. Our analysis is based on the computation of the marginal likelihood and the Bayes factor for the two damping models
      Arregui, I. et al.

      Fecha de publicación:

      8
      2023
      Número de citas
      4
    • Opacity for realistic 3D MHD simulations of cool stellar atmospheres
      Context. Realistic three-dimensional time-dependent simulations of stellar near-surface convection employ the opacity binning method for the efficient and accurate computation of the radiative energy exchange. The method provides several orders of magnitude of speedup, but its implementation includes a number of free parameters. Aims: Our aim is to
      Perdomo García, A. et al.

      Fecha de publicación:

      7
      2023
      Número de citas
      4
    • A comparative study of resistivity models for simulations of magnetic reconnection in the solar atmosphere
      Context. Magnetic reconnection is a fundamental mechanism in astrophysics. A common challenge in mimicking this process numerically in particular for the Sun is that the solar electrical resistivity is small compared to the diffusive effects caused by the discrete nature of codes. Aims: We aim to study different anomalous resistivity models and
      Færder, Ø. H. et al.

      Fecha de publicación:

      7
      2023
      Número de citas
      4
    • A study of the capabilities for inferring atmospheric information from high-spatial-resolution simulations
      In this work, we study the accuracy that can be achieved when inferring the atmospheric information from realistic numerical magneto-hydrodynamic simulations that reproduce the spatial resolution we will obtain with future observations made by the 4m class telescopes DKIST and EST. We first study multiple inversion configurations using the SIR code
      Quintero Noda, C. et al.

      Fecha de publicación:

      7
      2023
      Número de citas
      7
    • Interaction of solar jets with filaments: Triggering of large-amplitude filament oscillations
      Context. Large-amplitude oscillations (LAOs) are often detected in filaments. Using multi-wavelength observations, their origin can be traced back to the interaction with eruptions and jets. Aims: We present two different case studies as observational evidence in support of 2.5D numerical magnetohydrodynamics (MHD) experiments that show that the
      Joshi, Reetika et al.

      Fecha de publicación:

      4
      2023
      Número de citas
      4
    • The Chromosphere Underneath a Coronal Bright Point
      Coronal bright points (CBPs) are sets of small-scale coronal loops, connecting opposite magnetic polarities, primarily characterized by their enhanced extreme-ultraviolet (EUV) and X-ray emission. Being ubiquitous, they are thought to play an important role in heating the solar corona. We aim at characterizing the barely explored chromosphere
      Bose, Souvik et al.

      Fecha de publicación:

      2
      2023
      Número de citas
      11
    • The role of cooling induced by mixing in the mass and energy cycles of the solar atmosphere
      In many astrophysical systems, mixing between cool and hot temperature gas/plasma through Kelvin-Helmholtz-instability-driven turbulence leads to the formation of an intermediate temperature phase with increased radiative losses that drive efficient cooling. The solar atmosphere is a potential site for this process to occur with interaction between
      Hillier, Andrew et al.

      Fecha de publicación:

      4
      2023
      Número de citas
      6
    • Large Ion-neutral Drift Velocities and Plasma Heating in Partially Ionized Coronal Rain Blobs
      In this paper we present a numerical study of the dynamics of partially ionized coronal rain blobs. We use a two-fluid model to perform a high-resolution 2D simulation that takes into account the collisional interaction between the charged and neutral particles contained in the plasma. We follow the evolution of a cold plasma condensation as it
      Martínez-Gómez, David et al.

      Fecha de publicación:

      12
      2022
      Número de citas
      7
    • Observational and numerical characterization of a recurrent arc-shaped front propagating along a coronal fan
      Context. Recurrent, arc-shaped intensity disturbances were detected by extreme-ultraviolet channels in an active region. The fronts were observed to propagate along a coronal loop bundle rooted in a small area within a sunspot umbra. Previous works have linked these intensity disturbances to slow magnetoacoustic waves that propagate from the lower
      Sieyra, M. V. et al.

      Fecha de publicación:

      11
      2022
      Número de citas
      2
    • Ambipolar Diffusion in the Lower Solar Atmosphere: Magnetohydrodynamic Simulations of a Sunspot
      Magnetohydrodynamic (MHD) simulations of the solar atmosphere are often performed under the assumption that the plasma is fully ionized. However, in the lower solar atmosphere a reduced temperature often results in only the partial ionization of the plasma. The interaction between the decoupled neutral and ionized components of such a partially
      MacBride, Conor D. et al.

      Fecha de publicación:

      10
      2022
      Número de citas
      5
    • The European Solar Telescope
      The European Solar Telescope (EST) is a project aimed at studying the magnetic connectivity of the solar atmosphere, from the deep photosphere to the upper chromosphere. Its design combines the knowledge and expertise gathered by the European solar physics community during the construction and operation of state-of-the-art solar telescopes
      Quintero Noda, C. et al.

      Fecha de publicación:

      10
      2022
      Número de citas
      44
    • A 2D Model for Coronal Bright Points: Association with Spicules, UV Bursts, Surges, and EUV Coronal Jets
      Coronal bright points (CBPs) are ubiquitous structures in the solar atmosphere composed of hot small-scale loops observed in extreme-ultraviolet (EUV) or X-rays in the quiet Sun and coronal holes. They are key elements to understanding the heating of the corona; nonetheless, basic questions regarding their heating mechanisms, the chromosphere
      Nóbrega-Siverio, D. et al.

      Fecha de publicación:

      8
      2022
      Número de citas
      16
    • Doppler-velocity Drifts Detected in a Solar Prominence
      We analyzed multiline observations of a quiescent prominence from the slit spectrograph located at the Ondřejov Observatory. Dopplergrams and integrated intensity maps of the whole prominence were obtained from observations in six spectral lines: Ca II H, Hϵ, Hβ, He I D3, Hα, and Ca II IR. By combining integrated intensity maps with non-LTE
      Zapiór, Maciej et al.

      Fecha de publicación:

      7
      2022
      Número de citas
      12
    • Modeling the thermal conduction in the solar atmosphere with the code MANCHA3D
      Context. Thermal conductivity is one of the important mechanisms of heat transfer in the solar corona. In the limit of strongly magnetized plasma, it is typically modeled by Spitzer's expression where the heat flux is aligned with the magnetic field. Aims: This paper describes the implementation of the heat conduction into the code MANCHA3D with an
      Navarro, A. et al.

      Fecha de publicación:

      7
      2022
      Número de citas
      10
    • Generalized Fluid Models of the Braginskii Type
      Several generalizations of the well-known fluid model of Braginskii (1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multispecies plasmas with arbitrary
      Hunana, P. et al.

      Fecha de publicación:

      6
      2022
      Número de citas
      13
    • Ambipolar diffusion: Self-similar solutions and MHD code testing. Cylindrical symmetry
      Context. Ambipolar diffusion is a process occurring in partially ionised astrophysical systems that imparts a complicated mathematical and physical nature to Ohm's law. The numerical codes that solve the magnetohydrodynamic (MHD) equations have to be able to deal with the singularities that are naturally created in the system by the ambipolar
      Moreno-Insertis, F. et al.

      Fecha de publicación:

      6
      2022
      Número de citas
      3
    • Construction of coronal hole and active region magnetohydrostatic solutions in two dimensions: Force and energy balance
      Coronal holes and active regions are typical magnetic structures found in the solar atmosphere. We propose several magnetohydrostatic equilibrium solutions that are representative of these structures in two dimensions. Our models include the effect of a finite plasma-β and gravity, but the distinctive feature is that we incorporate a thermal
      Terradas, J. et al.

      Fecha de publicación:

      4
      2022
      Número de citas
      5
    • Recent Applications of Bayesian Methods to the Solar Corona
      Solar coronal seismology is based on the remote diagnostics of physical conditions in the corona of the Sun by comparison between model predictions and observations of magnetohydrodynamic wave activity. Our lack of direct access to the physical systems of interest makes information incomplete and uncertain so our conclusions are at best
      Arregui, I.

      Fecha de publicación:

      3
      2022
    • Methodology for Predicting the Probability Distribution of the Amplitude of Solar Cycle 25
      A number of precursor-type methods for solar-cycle prediction are based on the use of regression models and confidence-level estimates. A drawback of these methods is that they do not permit one to make probability statements, nor do they offer straightforward ways to propagate the uncertainty from observations to the quantities of interest. We
      Arregui, Iñigo

      Fecha de publicación:

      2
      2022
      Número de citas
      2

    Charlas relacionadas

    No se han encontrado charlas relacionadas.

    Congresos relacionados

    No se han encontrado congresos relacionados.