Numerical Simulation of Astrophysical Processes

Start year
2003
Organizational Unit

Grants related:

    General
    Description

    Numerical simulation through complex computer codes has been a fundamental tool in physics and technology research for decades. The rapid growth of computing capabilities, coupled with significant advances in numerical mathematics, has made this branch of research accessible to medium-sized research centers, bridging the gap between theoretical and experimental physics. Astrophysics is no exception to this trend. Since the late 1970s, a specialized field known as computational astrophysics has emerged, allowing us to understand a wide range of phenomena that were previously inaccessible to pure theoretical research and to account for previously unexplained observations.

    In recent decades, its primary areas of application have included (magneto)hydrodynamic phenomena and gas dynamics in various cosmic environments. For example, this includes the interiors and atmospheres of stars and planets, the interstellar medium, including magnetohydrodynamics and dynamos, accretion disks, the evolution of planetary nebulae, supernova explosions and remnants, and more. The incorporation of radiative transport equations into numerical simulations, which occurred in past decades, has added greater realism to the study of hydrodynamic processes in stellar photospheres and chromospheres.

    The current project aims to support the development of astrophysical research at the IAC based on the use of large numerical codes that require massively parallel computers and their connection with observational results. The general objective of this project is to perform calculations related to cosmic fluid dynamics and radiative transport. The topics of these calculations will focus on:

    1. Magnetized gas dynamics in the interiors and atmospheres of stars.
    2. Radiation transport and polarization signals in spectral lines based on realistic atomic and molecular models, including Hanle and Zeeman effects.
    3. Comparing theoretical/numerical results with observational data.

    This project is particularly relevant given the increasing involvement of the IAC in national and European supercomputing networks and, more generally, in large-scale supercomputer installation initiatives.

    Principal investigator

    In the following, we highlight the results of our annual year-end summary.

    Throughout the year 2022, partial ionization effects, nonequilibrium ionization effects, and multi-fluids have been one the main blocks of development both from the theoretical and numerical perspective. For instance, a generalization of the the Braginskii 1965 equations has been achieved for a general multi-species plasmas with arbitrary masses and temperatures, and where all of the viscosities and heat fluxes in the model are described by their own evolution equations. This new approach has a crucial advantage that the parallel components along the magnetic field lines do not become unbounded (infinitely large) in regimes of low-collisionallity of interest for this group as, for example, the solar corona (Hunana et al. 2022).  In this thematic block, 2D and 3D simulations, using a two-fluid model that treats the neutral and ionized species as two separate components, have also been performed to analyze the effect that the collisional interaction between both components has on the dynamics of coronal rain, the evolution of the instability of Kelvin-Helmholtz, the propagation of magneto-acoustic waves through the solar chromosphere or the heating of the plasma (Martínez-Gómez et al. 2022a). Another example of the theoretical development with potential numerical applications has been the pursuit of the effects of the ambipolar diffusion in the chromosphere from a more fundamental perspective by means of analytical solutions. The obtained solutions for cases with cylindrical symmetry are shown to constitute a demanding, but nonetheless viable, test for magnetohydrodynamic (MHD) codes that incorporate ambipolar diffusion. In addition, detailed tabulated runs of the solutions have been made available public for the community (Moreno-Insertis et al. 2022).  Lastly, nonequilibrium ionization effects of the Hydrogen atom together with the study of the Lyman α effects have been started to study in simple configurations to be applied later in realistic simulations that include the chromosphere.

     

    Improving and testing the capabilities of the available MHD codes in the solar group has been another of the major key developments carried out in 2022. For example, the results obtained by Moreno-Insertis et al. 2022 were used to check that the MHD Bifrost code is able to reproduce the theoretical solutions with sufficient accuracy up to very advanced diffusive times, as well as to explore the asymptotic properties of these theoretical solutions. In addition to that, several changes have been performed in the MANCHA code whose aim was to increase the efficiency and to add new features that will allow the researchers to perform more realistic experiments as well as exploring new research areas. For instance, MANCHA code has been extended to be able to simulate solar simulations up to the corona, adding a new module that efficiently calculates one of the key ingredients in the corona: the thermal conduction (Navarro et al. 2022). The preparation of the MANCHA code for its multi-fluid extension with radiation has also been another working branch concerning the numerical development in 2022. In addition, new equation-of-state and opacity routines have been developed that allow separating the equilibrium background contributions from those treated out of the equilibrium. Besides facing different challenges in solar physics, the huge development brought about in MANCHA is useful to study main sequence cool stars (G,K,M), which contributes to the better understanding of the stellar physics. To accomplish all these tasks, it was necessary not only to carry out numerous scaling tests and numerical experiments in local machines at the IAC, as well as on Supercomputers such as LaPalma, PICASSO, PizDaint, and MareNostrum4; but also to work together with external collaborators.

     

    During 2022, in this project there has also been a focus on different solar atmosphere phenomena and the corresponding comparison with observations. As an illustrative example, Coronal Bright Points (CBPs) have been modeled for the first time with enough realism to unravel the mechanisms that generate them and provide them with energy, being also able to explain different characteristics observed from space satellites. The comparison with observations is through synthetic SDO/AIA, Solar Orbiter EUI-HRI, and IRIS images that have been computed from the numerical experiment performed with the Bifrost code (Nóbrega-Siverio and Moreno-Insertis, 2022). Another example is the combination of 3D numerical experiments with the MoLMH code and forward modelling using Hα line to study transverse kink oscillations in prominence threads. The results contain relevant implications for the field of prominence seismology, showing that the Hα emission can be used to detect the fundamental mode of the oscillations (Martínez-Gómez et al. 2022b). In addition, ground high-resolution observations of ejective phenomena such as surges in the solar atmosphere have been analyzed, finding striking similarities with results obtained from numerical experiments. On top of that, there have also been significant contributions from the members of this project to the further advance of the observations and construction of new telescopes (Quintero et al. 2022) and satellites (De Pontieu et al. 2022, Cheung et al. 2022), using the earned knowledge from the theoretical-numerical experiments. Finally, an exploratory first attempt at understanding the physics of coronal holes and active regions from a global point of view through 2D magnetohydrostatic solutions was performed (Terradas et al. 2022), which will need of further development in the incoming years for comparisons with observations.

     

    Last but not least, state-of-the-art tools such as the ones provided by Machine Learning and Bayesian statistics have been applied to solar atmosphere problems. In this vein, a project to characterize the limits of the k-means methods and its application to solar observations was launched. In addition, new development in radiative transfer codes have being started to use in a preliminary study of machine learning approach to the computations of radiative terms. Development of the application of Bayesian techniques to the comparison of models in seismology of the solar atmosphere continued in 2022, with a review article published that accounts for the main results obtained in the last decade (Arregui 2022a). Moreover, the Bayesian formalism has been successfully applied to the prediction of the amplitude of the solar activity cycle, proposing a new methodology to quantify the goodness of both the prediction and the underlying model (Arregui 2022b).

    Related publications

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

      Advertised on:

      4
      2017
      Citations
      23
    • 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.

      Advertised on:

      1
      2017
      Citations
      29
    • 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.

      Advertised on:

      11
      2015
      Citations
      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.

      Advertised on:

      10
      2015
      Citations
      24
    • 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.

      Advertised on:

      7
      2015
      Citations
      5
    • 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.

      Advertised on:

      7
      2015
      Citations
      22
    • 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.

      Advertised on:

      7
      2015
      Citations
      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.

      Advertised on:

      7
      2015
      Citations
      60
    • 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.

      Advertised on:

      6
      2015
      Citations
      11
    • 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.

      Advertised on:

      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.

      Advertised on:

      5
      2015
      Citations
      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.

      Advertised on:

      4
      2015
      Citations
      48
    • 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.

      Advertised on:

      4
      2015
      Citations
      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.

      Advertised on:

      4
      2015
      Citations
      19
    • 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.

      Advertised on:

      11
      2014
      Citations
      27
    • 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.

      Advertised on:

      5
      2014
      Citations
      19
    • 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.

      Advertised on:

      4
      2014
      Citations
      905
    • 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.

      Advertised on:

      3
      2014
      Citations
      9
    • 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.

      Advertised on:

      11
      2013
      Citations
      90
    • 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.

      Advertised on:

      9
      2013
      Citations
      74

    Related talks

    No related talks were found.

    Related conferences

    No related conferences were found.