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

    • The C/O ratio at low metallicity: constraints on early chemical evolution from observations of Galactic halo stars

      Aims: We present new measurements of the abundances of carbon and oxygen derived from high-excitation C i and O i absorption lines in metal-poor halo stars, with the aim of clarifying the main sources of these two elements in the early stages of the chemical enrichment of the Galaxy. Methods: We target 15 new stars compared to our previous study

      Fabbian, D. et al.

      Fecha de publicación:

      6
      2009
      Número de citas
      102
    • Stellar Kinematics in Double-Barred Galaxies: The σ-Hollows

      We present SAURON integral-field stellar velocity and velocity dispersion maps for four double-barred early-type galaxies: NGC 2859, NGC 3941, NGC 4725, and NGC 5850. The presence of the inner bar does not produce major changes in the line-of-sight velocity, but it appears to have an important effect in the stellar velocity dispersion maps: we find

      de Lorenzo-Cáceres, A. et al.

      Fecha de publicación:

      9
      2008
      Número de citas
      34
    • Solar Abundance Corrections Derived Through Three-dimensional Magnetoconvection Simulations

      We explore the effect of the magnetic field when using realistic three-dimensional convection experiments to determine solar element abundances. By carrying out magnetoconvection simulations with a radiation-hydro code (the Copenhagen stagger code) and through a posteriori spectral synthesis of three Fe I lines, we obtain evidence that moderate

      Fabbian, D. et al.

      Fecha de publicación:

      12
      2010
      Número de citas
      58
    • Neutral oxygen spectral line formation revisited with new collisional data: large departures from LTE at low metallicity

      Aims: A detailed study is presented, including estimates of the impact on elemental abundance analysis, of the non-local thermodynamic equilibrium (non-LTE) formation of the high-excitation neutral oxygen 777 nm triplet in model atmospheres representative of stars with spectral types F to K. Methods: We have applied the statistical equilibrium code

      Fabbian, D. et al.

      Fecha de publicación:

      6
      2009
      Número de citas
      84
    • The Electrical Current Density Vector in the Inner Penumbra of a Sunspot

      We determine the entire electrical current density vector in a geometrical three-dimensional volume of the inner penumbra of a sunspot from an inversion of spectropolarimetric data obtained with Hinode/SP. Significant currents are seen to wrap around the hotter, more elevated regions with lower and more horizontal magnetic fields that harbor strong

      Puschmann, K. G. et al.

      Fecha de publicación:

      9
      2010
      Número de citas
      21
    • Magnetoacoustic Wave Energy from Numerical Simulations of an Observed Sunspot Umbra

      We aim at reproducing the height dependence of sunspot wave signatures obtained from spectropolarimetric observations through three-dimensional MHD numerical simulations. A magnetostatic sunspot model based on the properties of the observed sunspot is constructed and perturbed at the photosphere, introducing the fluctuations measured with the Si I

      Felipe, T. et al.

      Fecha de publicación:

      7
      2011
      Número de citas
      34
    • A Geometrical Height Scale for Sunspot Penumbrae

      Inversions of spectropolarimetric observations of penumbral filaments deliver the stratification of different physical quantities in an optical depth scale. However, without establishing a geometrical height scale, their three-dimensional geometrical structure cannot be derived. This is crucial in understanding the correct spatial variation of

      Puschmann, K. G. et al.

      Fecha de publicación:

      9
      2010
      Número de citas
      43
    • Three-dimensional Radiative Transfer Modeling of the Polarization of the Sun's Continuous Spectrum

      Polarized light provides the most reliable source of information at our disposal for diagnosing the physical properties of astrophysical plasmas, including the three-dimensional (3D) structure of the solar atmosphere. Here we formulate and solve the 3D radiative transfer problem of the linear polarization of the solar continuous radiation, which is

      Trujillo Bueno, Javier et al.

      Fecha de publicación:

      4
      2009
      Número de citas
      42
    • Non-LTE Determination of the Silicon Abundance Using a Three-dimensional Hydrodynamical Model of the Solar Photosphere

      Confrontations of spectroscopic observations with local thermodynamic equilibrium (LTE) spectral syntheses in a three-dimensional (3D) hydrodynamical model of the solar photosphere led to a downward revision of the photospheric and meteoritic silicon abundances. Here we derive the photospheric silicon abundance taking into account non-LTE (NLTE)

      Shchukina, N. et al.

      Fecha de publicación:

      8
      2012
      Número de citas
      14
    • Gauss-Seidel and Successive Overrelaxation Methods for Radiative Transfer with Partial Frequency Redistribution

      The linearly polarized solar limb spectrum that is produced by scattering processes contains a wealth of information on the physical conditions and magnetic fields of the solar outer atmosphere, but the modeling of many of its strongest spectral lines requires solving an involved non-local thermodynamic equilibrium radiative transfer problem

      Sampoorna, M. et al.

      Fecha de publicación:

      4
      2010
      Número de citas
      14
    • Channeling 5 Minute Photospheric Oscillations into the Solar Outer Atmosphere through Small-Scale Vertical Magnetic Flux Tubes

      We report two-dimensional MHD simulations which demonstrate that photospheric 5 minute oscillations can leak into the chromosphere inside small-scale vertical magnetic flux tubes. The results of our numerical experiments are compatible with those inferred from simultaneous spectropolarimetric observations of the photosphere and chromosphere

      Khomenko, E. et al.

      Fecha de publicación:

      3
      2008
      Número de citas
      75
    • The Sunrise Mission

      The first science flight of the balloon-borne Sunrise telescope took place in June 2009 from ESRANGE (near Kiruna/Sweden) to Somerset Island in northern Canada. We describe the scientific aims and mission concept of the project and give an overview and a description of the various hardware components: the 1-m main telescope with its postfocus

      Barthol, P. et al.

      Fecha de publicación:

      1
      2011
      Número de citas
      196
    • The Frontier between Small-scale Bipoles and Ephemeral Regions in the Solar Photosphere: Emergence and Decay of an Intermediate-scale Bipole Observed with SUNRISE/IMaX

      We report on the photospheric evolution of an intermediate-scale (≈4 Mm footpoint separation) magnetic bipole, from emergence to decay, observed in the quiet Sun at high spatial (0farcs3) and temporal (33 s) resolution. The observations were acquired by the Imaging Magnetograph Experiment imaging magnetograph during the first science flight of the

      Guglielmino, S. L. et al.

      Fecha de publicación:

      2
      2012
      Número de citas
      40
    • The Filter Imager SuFI and the Image Stabilization and Light Distribution System ISLiD of the Sunrise Balloon-Borne Observatory: Instrument Description

      We describe the design of the Sunrise Filter Imager (SuFI) and the Image Stabilization and Light Distribution (ISLiD) unit onboard the Sunrise balloon borne solar observatory. This contribution provides the necessary information which is relevant to understand the instruments' working principles, the relevant technical data, and the necessary

      Gandorfer, A. et al.

      Fecha de publicación:

      1
      2011
      Número de citas
      94
    • Mesogranulation and the Solar Surface Magnetic Field Distribution

      The relation of the solar surface magnetic field with mesogranular cells is studied using high spatial (≈100 km) and temporal (≈30 s) resolution data obtained with the IMaX instrument on board SUNRISE. First, mesogranular cells are identified using Lagrange tracers (corks) based on horizontal velocity fields obtained through local correlation

      Yelles-Chaouche, L. et al.

      Fecha de publicación:

      2
      2011
      Número de citas
      43
    • Magnetic field emergence in mesogranular-sized exploding granules observed with sunrise/IMaX data

      We report on magnetic field emergences covering significant areas of exploding granules. The balloon-borne mission Sunrise provided high spatial and temporal resolution images of the solar photosphere. Continuum images, longitudinal and transverse magnetic field maps and Dopplergrams obtained by IMaX onboard Sunrise are analyzed by local

      Palacios, J. et al.

      Fecha de publicación:

      1
      2012
      Número de citas
      24
    • Center-to-limb variation of the area covered by magnetic bright points in the quiet Sun

      Context. The quiet Sun magnetic fields produce ubiquitous bright points (BPs) that cover a significant fraction of the solar surface. Their contribution to the total solar irradiance (TSI) is so-far unknown. Aims: We aim at measuring the center-to-limb variation (CLV) of the fraction of solar surface covered by quiet Sun magnetic bright points. The

      Bonet, J. A. et al.

      Fecha de publicación:

      3
      2012
      Número de citas
      15
    • Pattern Speeds of Bars and Spiral Arms from Hα Velocity Fields

      We have applied the Tremaine-Weinberg method to 10 late-type barred spiral galaxies using data cubes, in Hα emission, from the FaNTOmM and GHAFAS Fabry-Perot spectrometers. We have combined the derived bar (and/or spiral) pattern speeds with angular frequency plots to measure the corotation radii for the bars in these galaxies. We base our results

      Fathi, K. et al.

      Fecha de publicación:

      10
      2009
      Número de citas
      59
    • Stray-light contamination and spatial deconvolution of slit-spectrograph observations

      Context. Stray light caused by scattering on optical surfaces and in the Earth's atmosphere degrades the spatial resolution of observations. Whereas post-facto reconstruction techniques are common for 2D imaging and spectroscopy, similar options for slit-spectrograph data are rarely applied. Aims: We study the contribution of stray light to the two

      Beck, C. et al.

      Fecha de publicación:

      11
      2011
      Número de citas
      34
    • Magneto-acoustic waves in sunspots from observations and numerical simulations

      We study the propagation of waves from the photosphere to the chromosphere of sunspots. From time series of cospatial Ca II H (including its line blends) intensity spectra and polarimetric spectra of Si I λ 1082.7 nm and He I λ 1083.0 nm we retrieve the line-of-sight velocity at several heights. The analysis of the phase difference and

      Felipe, T. et al.

      Fecha de publicación:

      1
      2011
      Número de citas
      0

    Charlas relacionadas

    No se han encontrado charlas relacionadas.

    Congresos relacionados

    No se han encontrado congresos relacionados.