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

    • Thermodynamic fluctuations in solar photospheric three-dimensional convection simulations and observations

      Context. Numerical three-dimensional (3D) radiative (magneto-)hydrodynamical [(M)HD] simulations of solar convection are nowadays used to understand the physical properties of the solar photosphere and convective envelope, and, in particular, to determine the Sun's photospheric chemical abundances. To validate this approach, it is important to

      Beck, C. et al.

      Fecha de publicación:

      9
      2013
      Número de citas
      17
    • Physical properties of a sunspot chromosphere with umbral flashes

      We present new high-resolution spectro-polarimetric Ca IIλ8542 observations of umbral flashes in sunspots. At nearly 0.18 arcsec, and spanning about one hour of continuous observation, this is the most detailed dataset published thus far. Our study involves both LTE and non-LTE inversions (but includes also a weak field analysis as a sanity check)

      de la Cruz Rodríguez, J. et al.

      Fecha de publicación:

      8
      2013
      Número de citas
      102
    • Improved Search of Principal Component Analysis Databases for Spectro-polarimetric Inversion

      We describe a simple technique for the acceleration of spectro-polarimetric inversions based on principal component analysis (PCA) of Stokes profiles. This technique involves the indexing of the database models based on the sign of the projections (PCA coefficients) of the first few relevant orders of principal components of the four Stokes

      Casini, R. et al.

      Fecha de publicación:

      8
      2013
      Número de citas
      8
    • Target Selection for the Apache Point Observatory Galactic Evolution Experiment (APOGEE)

      The Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution infrared spectroscopic survey spanning all Galactic environments (i.e., bulge, disk, and halo), with the principal goal of constraining dynamical and chemical evolution models of the Milky Way. APOGEE takes advantage of the reduced effects of extinction at

      Zasowski, G. et al.

      Fecha de publicación:

      10
      2013
      Número de citas
      329
    • Plasma Jets and Eruptions in Solar Coronal Holes: A Three-dimensional Flux Emergence Experiment

      A three-dimensional (3D) numerical experiment of the launching of a hot and fast coronal jet followed by several violent eruptions is analyzed in detail. These events are initiated through the emergence of a magnetic flux rope from the solar interior into a coronal hole. We explore the evolution of the emerging magnetically dominated plasma dome

      Moreno-Insertis, F. et al.

      Fecha de publicación:

      7
      2013
      Número de citas
      140
    • Atmosphere Dynamics of the Active Region NOAA 11024

      We present results of the study of chromospheric and photospheric line-of-sight velocity fields in the young active region NOAA 11024. Multi-layer, multi-wavelength observational data were used for the analysis of the emerging flux in this active region. Spectropolarimetric observations were carried out with the telescope THEMIS on Tenerife (Canary

      Khomenko, E. V. et al.

      Fecha de publicación:

      6
      2013
      Número de citas
      16
    • Is Magnetic Reconnection the Cause of Supersonic Upflows in Granular Cells?

      In a previous work, we reported on the discovery of supersonic magnetic upflows on granular cells in data from the SUNRISE/IMaX instrument. In the present work, we investigate the physical origin of these events employing data from the same instrument but with higher spectral sampling. By means of the inversion of Stokes profiles we are able to

      Borrero, J. M. et al.

      Fecha de publicación:

      5
      2013
      Número de citas
      21
    • Bayesian Analysis of Multiple Harmonic Oscillations in the Solar Corona

      The detection of multiple mode harmonic kink oscillations in coronal loops enables us to obtain information on coronal density stratification and magnetic field expansion using seismology inversion techniques. The inference is based on the measurement of the period ratio between the fundamental mode and the first overtone and theoretical results

      Díaz, A. J. et al.

      Fecha de publicación:

      3
      2013
      Número de citas
      20
    • Constraining clumpy dusty torus models using optimized filter sets

      Recent success in explaining several properties of the dusty torus around the central engine of active galactic nuclei has been gathered with the assumption of clumpiness. The properties of such clumpy dusty tori can be inferred by analysing spectral energy distributions (SEDs), sometimes with scarce sampling given that large aperture telescopes

      Ramos-Almeida, C. et al.

      Fecha de publicación:

      1
      2013
      Número de citas
      11
    • Solar Fe abundance and magnetic fields. Towards a consistent reference metallicity

      Aims: We investigate the impact on Fe abundance determination of including magnetic flux in series of 3D radiation-magnetohydrodynamics (MHD) simulations of solar convection, which we used to synthesize spectral intensity profiles corresponding to disc centre. Methods: A differential approach is used to quantify the changes in theoretical

      Fabbian, D. et al.

      Fecha de publicación:

      12
      2012
      Número de citas
      42
    • Analytic Approximate Seismology of Propagating Magnetohydrodynamic Waves in the Solar Corona

      Observations show that propagating magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. The technique of MHD seismology uses the wave observations combined with MHD wave theory to indirectly infer physical parameters of the solar atmospheric plasma and magnetic field. Here, we present an analytical seismological inversion scheme

      Goossens, M. et al.

      Fecha de publicación:

      12
      2012
      Número de citas
      19
    • Shear Photospheric Forcing and the Origin of Turbulence in Coronal Loops

      We present a series of numerical simulations aimed at understanding the nature and origin of turbulence in coronal loops in the framework of the Parker model for coronal heating. A coronal loop is studied via reduced magnetohydrodynamic (MHD) simulations in Cartesian geometry. A uniform and strong magnetic field threads the volume between the two

      Rappazzo, A. F. et al.

      Fecha de publicación:

      10
      2010
      Número de citas
      50
    • Simulation of a flux emergence event and comparison with observations by Hinode

      Aims: We study the observational signature of flux emergence in the photosphere using synthetic data from a 3D MHD simulation of the emergence of a twisted flux tube. Methods: Several stages in the emergence process are considered. At every stage we compute synthetic Stokes spectra of the two iron lines Fe I 6301.5 Å and Fe I 6302.5 Å and degrade

      Yelles Chaouche, L. et al.

      Fecha de publicación:

      12
      2009
      Número de citas
      15
    • Nonlinear force-free modelling: influence of inaccuracies in the measured magnetic vector

      Context. Solar magnetic fields are regularly extrapolated into the corona starting from photospheric magnetic measurements that can be affected by significant uncertainty. Aims: We study how inaccuracies introduced into the maps of the photospheric magnetic vector by the inversion of ideal and noisy Stokes parameters influence the extrapolation of

      Wiegelmann, T. et al.

      Fecha de publicación:

      2
      2010
      Número de citas
      18
    • Comparison of the thin flux tube approximation with 3D MHD simulations

      Context: The structure and dynamics of small vertical photospheric magnetic flux concentrations has been often treated in the framework of an approximation based upon a low-order truncation of the Taylor expansions of all quantities in the horizontal direction, together with the assumption of instantaneous total pressure balance at the boundary to

      Yelles Chaouche, L. et al.

      Fecha de publicación:

      9
      2009
      Número de citas
      30
    • The emergence of toroidal flux tubes from beneath the solar photosphere

      Context: Models of flux emergence frequently use a twisted cylindrical loop as the initial starting configuration and ignore the coupling between the radiation field and plasma. In these models, the axis of the original tube never emerges through the photosphere. Without the axis emerging, it is very difficult to form a realistic sunspot. Aims: The

      Hood, A. W. et al.

      Fecha de publicación:

      9
      2009
      Número de citas
      64
    • On the Origin of the Type II Spicules: Dynamic Three-dimensional MHD Simulations

      Recent high temporal and spatial resolution observations of the chromosphere have forced the definition of a new type of spicule, "type II's," that are characterized by rising rapidly, having short lives, and by fading away at the end of their lifetimes. Here, we report on features found in realistic three-dimensional simulations of the outer solar

      Martínez-Sykora, Juan et al.

      Fecha de publicación:

      7
      2011
      Número de citas
      78
    • Magnetic flux emergence into the solar photosphere and chromosphere

      Aims: We model the emergence of magnetized plasma across granular convection cells and the low atmosphere, including layers up to the mid-chromosphere. Methods: Three-dimensional numerical experiments are carried out in which the equations of MHD and radiative transfer are solved self-consistently. We use the MURaM code, which assumes local

      Tortosa-Andreu, A. et al.

      Fecha de publicación:

      11
      2009
      Número de citas
      65
    • Jets in Coronal Holes: Hinode Observations and Three-dimensional Computer Modeling

      Recent observations of coronal hole areas with the XRT and EIS instruments on board the Hinode satellite have shown with unprecedented detail the launching of fast, hot jets away from the solar surface. In some cases these events coincide with episodes of flux emergence from beneath the photosphere. In this Letter we show results of a three

      Moreno-Insertis, F. et al.

      Fecha de publicación:

      2
      2008
      Número de citas
      202
    • The Three-dimensional Structure of an Active Region Filament as Extrapolated from Photospheric and Chromospheric Observations

      The three-dimensional structure of an active region filament is studied using nonlinear force-free field extrapolations based on simultaneous observations at a photospheric and a chromospheric height. To that end, we used the Si I 10827 Å line and the He I 10830 Å triplet obtained with the Tenerife Infrared Polarimeter at the Vacuum Tower Telescope

      Yelles-Chaouche, L. et al.

      Fecha de publicación:

      3
      2012
      Número de citas
      30

    Charlas relacionadas

    No se han encontrado charlas relacionadas.

    Congresos relacionados

    No se han encontrado congresos relacionados.