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

    • Self-similar Approach for Rotating Magnetohydrodynamic Solar and Astrophysical Structures
      Rotating magnetic structures are common in astrophysics, from vortex tubes and tornadoes in the Sun all the way to jets in different astrophysical systems. The physics of these objects often combine inertial, magnetic, gas pressure, and gravitational terms. Also, they often show approximate symmetries that help simplify the otherwise rather
      Luna, M. et al.

      Fecha de publicación:

      8
      2018
      Número de citas
      2
    • MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma
      The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere and low chromosphere. Under the single-fluid formalism the main non-ideal effects are described by a series of
      González-Morales, P. A. et al.

      Fecha de publicación:

      7
      2018
      Número de citas
      30
    • Inversions of synthetic umbral flashes: Effects of scanning time on the inferred atmospheres
      Context. The use of instruments that record narrowband images at selected wavelengths is a common approach in solar observations. They allow scanning of a spectral line by sampling the Stokes profiles with two-dimensional images at each line position, but require a compromise between spectral resolution and temporal cadence. The interpretation and
      Felipe, T. et al.

      Fecha de publicación:

      6
      2018
      Número de citas
      20
    • Small-scale Magnetic Flux Emergence in the Quiet Sun
      Small bipolar magnetic features are observed to appear in the interior of individual granules in the quiet Sun, signaling the emergence of tiny magnetic loops from the solar interior. We study the origin of those features as part of the magnetoconvection process in the top layers of the convection zone. Two quiet-Sun magnetoconvection models
      Moreno-Insertis, F. et al.

      Fecha de publicación:

      6
      2018
      Número de citas
      26
    • GONG Catalog of Solar Filament Oscillations Near Solar Maximum
      We have cataloged 196 filament oscillations from the Global Oscillation Network Group Hα network data during several months near the maximum of solar cycle 24 (2014 January–June). Selected examples from the catalog are described in detail, along with our statistical analyses of all events. Oscillations were classified according to their velocity
      Luna, M. et al.

      Fecha de publicación:

      6
      2018
      Número de citas
      38
    • On the Importance of the Nonequilibrium Ionization of Si IV and O IV and the Line of Sight in Solar Surges
      Surges are ubiquitous cool ejections in the solar atmosphere that often appear associated with transient phenomena like UV bursts or coronal jets. Recent observations from the Interface Region Imaging Spectrograph show that surges, although traditionally related to chromospheric lines, can exhibit enhanced emission in Si IV with brighter spectral
      Nóbrega-Siverio, D. et al.

      Fecha de publicación:

      5
      2018
      Número de citas
      20
    • Prominence oscillations
      Prominences are intriguing, but poorly understood, magnetic structures of the solar corona. The dynamics of solar prominences has been the subject of a large number of studies, and of particular interest is the study of prominence oscillations. Ground- and space-based observations have confirmed the presence of oscillatory motions in prominences
      Arregui, I. et al.

      Fecha de publicación:

      4
      2018
      Número de citas
      48
    • Rayleigh-Taylor instabilities with sheared magnetic fields in partially ionised plasmas
      Aims: In the present study we investigate the nature of the magnetic Rayleigh-Taylor instability appearing at a tangential discontinuity in a partially ionised plasma when the effect of magnetic shear is taken into account. Methods: The partially ionised character of the plasma is described by the ambipolar diffusion in the induction equation. The
      Ruderman, M. S. et al.

      Fecha de publicación:

      1
      2018
      Número de citas
      11
    • Large-amplitude Longitudinal Oscillations Triggered by the Merging of Two Solar Filaments: Observations and Magnetic Field Analysis
      We follow the eruption of two related intermediate filaments observed in Hα (from GONG) and EUV (from Solar Dynamics Observatory SDO/Atmospheric Imaging assembly AIA) and the resulting large-amplitude longitudinal oscillations of the plasma in the filament channels. The events occurred in and around the decayed active region AR12486 on 2016 January
      Luna, M. et al.

      Fecha de publicación:

      12
      2017
      Número de citas
      22
    • Surges and Si IV Bursts in the Solar Atmosphere: Understanding IRIS and SST Observations through RMHD Experiments
      Surges often appear as a result of the emergence of magnetized plasma from the solar interior. Traditionally, they are observed in chromospheric lines such as Hα 6563 \mathringA and Ca II 8542 \mathringA . However, whether there is a response to the surge appearance and evolution in the Si IV lines or, in fact, in many other transition region lines
      Nóbrega-Siverio, D. et al.

      Fecha de publicación:

      12
      2017
      Número de citas
      42
    • Intermittent Reconnection and Plasmoids in UV Bursts in the Low Solar Atmosphere
      Magnetic reconnection is thought to drive a wide variety of dynamic phenomena in the solar atmosphere. Yet, the detailed physical mechanisms driving reconnection are difficult to discern in the remote sensing observations that are used to study the solar atmosphere. In this Letter, we exploit the high-resolution instruments Interface Region Imaging
      Rouppe van der Voort, L. et al.

      Fecha de publicación:

      12
      2017
      Número de citas
      66
    • Magnetic topological analysis of coronal bright points
      Context. We report on the first of a series of studies on coronal bright points which investigate the physical mechanism that generates these phenomena. Aims: The aim of this paper is to understand the magnetic-field structure that hosts the bright points. Methods: We use longitudinal magnetograms taken by the Solar Optical Telescope with the
      Galsgaard, K. et al.

      Fecha de publicación:

      10
      2017
      Número de citas
      17
    • Comparison of Solar Fine Structure Observed Simultaneously in Lyα and Mg ii h
      The Chromospheric Lyman Alpha Spectropolarimeter (CLASP) observed the Sun in H i Lyα during a suborbital rocket flight on 2015 September 3. The Interface Region Imaging Telescope (IRIS) coordinated with the CLASP observations and recorded nearly simultaneous and co-spatial observations in the Mg ii h and k lines. The Mg ii h and Lyα lines are
      Schmit, D. et al.

      Fecha de publicación:

      10
      2017
      Número de citas
      7
    • Dissecting the long-term emission behaviour of the BL Lac object Mrk 421
      We report on long-term multiwavelength monitoring of blazar Mrk 421 by the GLAST-AGILE Support Program of the Whole Earth Blazar Telescope (GASP-WEBT) collaboration and Steward Observatory, and by the Swift and Fermi satellites. We study the source behaviour in the period 2007-2015, characterized by several extreme flares. The ratio between the
      Carnerero, M. I. et al.

      Fecha de publicación:

      12
      2017
      Número de citas
      33
    • Two-dimensional Radiative Magnetohydrodynamic Simulations of Partial Ionization in the Chromosphere. II. Dynamics and Energetics of the Low Solar Atmosphere
      We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion–neutral interaction effects using the generalized Ohm’s law, i.e., we include the Hall term and the ambipolar diffusion
      Martínez-Sykora, J. et al.

      Fecha de publicación:

      9
      2017
      Número de citas
      62
    • A Si I atomic model for NLTE spectropolarimetric diagnostics of the 10 827 Å line
      Aims: The Si i 10 827 Å line is commonly used for spectropolarimetric diagnostics of the solar atmosphere. First, we aim at quantifying the sensitivity of the Stokes profiles of this line to non-local thermodynamic equilibrium (NLTE) effects. Second, we aim at facilitating NLTE diagnostics of the Si i 10 827 Å line. To this end, we propose the use
      Shchukina, N. G. et al.

      Fecha de publicación:

      7
      2017
      Número de citas
      18
    • Numerical simulations of quiet Sun magnetic fields seeded by the Biermann battery
      The magnetic fields of the quiet Sun cover at any time more than 90% of its surface and their magnetic energy budget is crucial to explain the thermal structure of the solar atmosphere. One of the possible origins of these fields is the action of the local dynamo in the upper convection zone of the Sun. Existing simulations of the local solar
      Khomenko, E. et al.

      Fecha de publicación:

      8
      2017
      Número de citas
      38
    • Helioseismic holography of simulated sunspots: dependence of the travel time on magnetic field strength and Wilson depression
      Improving methods for determining the subsurface structure of sunspots from their seismic signature requires a better understanding of the interaction of waves with magnetic field concentrations. We aim to quantify the impact of changes in the internal structure of sunspots on local helioseismic signals. We have numerically simulated the
      Felipe, T. et al.

      Fecha de publicación:

      8
      2017
      Número de citas
      9
    • High-frequency waves in the corona due to null points
      This work aims to understand the behavior of non-linear waves in the vicinity of a coronal null point. In previous works we have shown that high-frequency waves are generated in such a magnetic configuration. This paper studies those waves in detail in order to provide a plausible explanation of their generation. We demonstrate that slow magneto
      Santamaria, I. C. et al.

      Fecha de publicación:

      6
      2017
      Número de citas
      11
    • Dependence of sunspot photospheric waves on the depth of the source of solar p-modes
      Photospheric waves in sunspots moving radially outward at speeds faster than the characteristic wave velocities have been recently detected. It has been suggested that they are the visual pattern of p-modes excited around 5 Mm beneath the sunspot's surface. Using numerical simulations, we performed a parametric study of the waves observed at the
      Felipe, T. et al.

      Fecha de publicación:

      2
      2017
      Número de citas
      12

    Charlas relacionadas

    No se han encontrado charlas relacionadas.

    Congresos relacionados

    No se han encontrado congresos relacionados.