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

    • Synchrotron emission from the blazar PG 1553+113. An analysis of its flux and polarization variability

      In 2015 July 29-September 1, the satellite XMM-Newton pointed at the BL Lac object PG 1553+133 six times, collecting data for 218 h. During one of these epochs, simultaneous observations by the Swift satellite were requested to compare the results of the X-ray and optical-UV instruments. Optical, near-infrared and radio monitoring was carried out

      Raiteri, C. M. et al.

      Fecha de publicación:

      4
      2017
      Número de citas
      22
    • On the effects of ion-neutral interactions in solar plasmas

      Solar photosphere and chromosphere are composed of weakly ionized plasma for which collisional coupling decreases with height. This implies a breakdown of some hypotheses underlying magnetohydrodynamics at low altitudes and gives rise to non-ideal MHD effects such as ambipolar diffusion, Hall effect, etc. Recently, there has been progress in

      Khomenko, E.

      Fecha de publicación:

      1
      2017
      Número de citas
      26
    • Multi-parametric Study of Rising 3D Buoyant Flux Tubes in an Adiabatic Stratification Using AMR

      We study the buoyant rise of magnetic flux tubes embedded in an adiabatic stratification using two-and three-dimensional, magnetohydrodynamic simulations. We analyze the dependence of the tube evolution on the field line twist and on the curvature of the tube axis in different diffusion regimes. To be able to achieve a comparatively high spatial

      Martínez-Sykora, J. et al.

      Fecha de publicación:

      11
      2015
      Número de citas
      24
    • Model Comparison for the Density Structure across Solar Coronal Waveguides

      The spatial variation of physical quantities, such as the mass density, across solar atmospheric waveguides governs the timescales and spatial scales for wave damping and energy dissipation. The direct measurement of the spatial distribution of density, however, is difficult, and indirect seismology inversion methods have been suggested as an

      Arregui, I. et al.

      Fecha de publicación:

      10
      2015
      Número de citas
      23
    • On the nature of transverse coronal waves revealed by wavefront dislocations

      Context. Coronal waves are an important aspect of the dynamics of the plasma in the corona. Wavefront dislocations are topological features of most waves in nature and also of magnetohydrodynamic waves. Are there dislocations in coronal waves? Aims: The finding and explanation of dislocations may shed light on the nature and characteristics of the

      López Ariste, A. et al.

      Fecha de publicación:

      7
      2015
      Número de citas
      4
    • Are Tornado-like Magnetic Structures Able to Support Solar Prominence Plasma?

      Recent high-resolution and high-cadence observations have surprisingly suggested that prominence barbs exhibit apparent rotating motions suggestive of a tornado-like structure. Additional evidence has been provided by Doppler measurements. The observations reveal opposite velocities for both hot and cool plasma on the two sides of a prominence barb

      Luna, M. et al.

      Fecha de publicación:

      7
      2015
      Número de citas
      21
    • The impact of surface dynamo magnetic fields on the solar iron abundance

      Most chemical abundance determinations ignore that the solar photosphere is significantly magnetized by the ubiquitous presence of a small-scale magnetic field. A recent investigation has suggested that there should be a significant impact on the derived iron abundance, owing to the magnetically induced changes on the photospheric temperature and

      Trujillo Bueno, J. et al.

      Fecha de publicación:

      7
      2015
      Número de citas
      14
    • Evidence for Evaporation-incomplete Condensation Cycles in Warm Solar Coronal Loops

      Quasi-constant heating at the footpoints of loops leads to evaporation and condensation cycles of the plasma: thermal non-equilibrium (TNE). This phenomenon is believed to play a role in the formation of prominences and coronal rain. However, it is often discounted as being involved in the heating of warm loops because the models do not reproduce

      Solomon, J. et al.

      Fecha de publicación:

      7
      2015
      Número de citas
      57
    • Model comparison for the density structure along solar prominence threads

      Context. Quiescent solar prominence fine structures are typically modelled as density enhancements, called threads, which occupy a fraction of a longer magnetic flux tube. This is justified from the spatial distribution of the imaged plasma emission or absorption of prominences at small spatial scales. The profile of the mass density along the

      Arregui, I. et al.

      Fecha de publicación:

      6
      2015
      Número de citas
      9
    • Wave heating of the solar atmosphere

      Magnetic waves are a relevant component in the dynamics of the solar atmosphere. Their significance has increased because of their potential as a remote diagnostic tool and their presumed contribution to plasma heating processes. We discuss our current understanding of coronal heating by magnetic waves, based on recent observational evidence and

      Arregui, I.

      Fecha de publicación:

      4
      2015
    • Magnetohydrodynamic wave propagation from the subphotosphere to the corona in an arcade-shaped magnetic field with a null point

      Aims: The aim of this work is to study the energy transport by means of Magnetohydrodynamic (MHD) waves propagating in quiet-Sun magnetic topology from layers below the surface to the corona. Upwardly propagating waves find obstacles, such as the equipartition layer with plasma β = 1, the transition region, and null points, and they get transmitted

      Santamaria, I. C. et al.

      Fecha de publicación:

      5
      2015
      Número de citas
      36
    • Three-dimensional Radiative Transfer Simulations of the Scattering Polarization of the Hydrogen Lyα Line in a Magnetohydrodynamic Model of the Chromosphere–Corona Transition Region

      Probing the magnetism of the upper solar chromosphere requires measuring and modeling the scattering polarization produced by anisotropic radiation pumping in UV spectral lines. Here we apply PORTA (a novel radiative transfer code) to investigate the hydrogen Lyα line in a three-dimensional model of the solar atmosphere resulting from a state of

      Štěpán, J. et al.

      Fecha de publicación:

      4
      2015
      Número de citas
      47
    • Height Variation of the Vector Magnetic Field in Solar Spicules

      Proving the magnetic configuration of solar spicules has hitherto been difficult due to the lack of spatial resolution and image stability during off-limb ground-based observations. We report spectropolarimetric observations of spicules taken in the He i 1083 nm spectral region with the Tenerife Infrared Polarimeter II at the German Vacuum Tower

      Orozco Suárez, D. et al.

      Fecha de publicación:

      4
      2015
      Número de citas
      28
    • Continuum Intensity and [O I] Spectral Line Profiles in Solar 3D Photospheric Models: The Effect of Magnetic Fields

      The importance of magnetic fields in three-dimensional (3D) magnetoconvection models of the Sun’s photosphere is investigated in terms of their influence on the continuum intensity at different viewing inclination angles and on the intensity profile of two [O i] spectral lines. We use the RH numerical radiative transfer code to perform a posteriori

      Fabbian, D. et al.

      Fecha de publicación:

      4
      2015
      Número de citas
      18
    • Synthetic Observations of Wave Propagation in a Sunspot Umbra

      Spectropolarimetric temporal series from Fe I λ6301.5 Å and Ca II infrared triplet lines are obtained by applying the Stokes synthesis code NICOLE to a numerical simulation of wave propagation in a sunspot umbra from MANCHA code. The analysis of the phase difference between Doppler velocity and intensity core oscillations of the Fe I λ6301.5 Å line

      Felipe, T. et al.

      Fecha de publicación:

      11
      2014
      Número de citas
      26
    • Determination of the cross-field density structuring in coronal waveguides using the damping of transverse waves

      Context. Time and spatial damping of transverse magnetohydrodynamic (MHD) kink oscillations is a source of information on the cross-field variation of the plasma density in coronal waveguides. Aims: We show that a probabilistic approach to the problem of determining the density structuring from the observed damping of transverse oscillations

      Arregui, I. et al.

      Fecha de publicación:

      5
      2014
      Número de citas
      17
    • The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

      The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the Tenth Public Data Release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon

      Anders, Friedrich et al.

      Fecha de publicación:

      4
      2014
      Número de citas
      899
    • The power spectrum of solar convection flows from high-resolution observations and 3D simulations

      Context. Understanding solar surface magnetoconvection requires the study of the Fourier spectra of the velocity fields. Nowadays, observations are available that resolve very small spatial scales, well into the subgranular range, almost reaching the scales routinely resolved in numerical magnetoconvection simulations. Comparison of numerical and

      Yelles Chaouche, L. et al.

      Fecha de publicación:

      3
      2014
      Número de citas
      8
    • Twisting solar coronal jet launched at the boundary of an active region

      Aims: A broad jet was observed in a weak magnetic field area at the edge of active region NOAA 11106 that also produced other nearby recurring and narrow jets. The peculiar shape and magnetic environment of the broad jet raised the question of whether it was created by the same physical processes of previously studied jets with reconnection

      Schmieder, B. et al.

      Fecha de publicación:

      11
      2013
      Número de citas
      88
    • PORTA: A three-dimensional multilevel radiative transfer code for modeling the intensity and polarization of spectral lines with massively parallel computers

      The interpretation of the intensity and polarization of the spectral line radiation produced in the atmosphere of the Sun and of other stars requires solving a radiative transfer problem that can be very complex, especially when the main interest lies in modeling the spectral line polarization produced by scattering processes and the Hanle and

      Štěpán, J. et al.

      Fecha de publicación:

      9
      2013
      Número de citas
      70

    Charlas relacionadas

    No se han encontrado charlas relacionadas.

    Congresos relacionados

    No se han encontrado congresos relacionados.