-
Mancha3D code: Multi-purpose Advanced Non-ideal MHD Code for High resolution simulations in Astrophysics
Authors:
M. Modestov,
E. Khomenko,
N. Vitas,
A. de Vicente,
A. Navarro,
P. A. Gonzalez-Morales,
M. Collados,
T. Felipe,
D. Martinez-Gomez,
P. Hunana,
M. Luna,
M. Koll Pistarini,
B. Popescu Braileanu,
A. Perdomo Garcia,
V. Liakh,
I. Santamaria,
M. M. Gomez Miguez
Abstract:
The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic processes in solar/stellar atmospheres. The code includes non-ideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high quality realistic simulations of magneto-convection, as well as idealized simulations of particular processes, suc…
▽ More
The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic processes in solar/stellar atmospheres. The code includes non-ideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high quality realistic simulations of magneto-convection, as well as idealized simulations of particular processes, such as wave propagation, instabilities or energetic events. The paper summarizes the equations and methods used in the Mancha3D code. It also describes its numerical stability and parallel performance and efficiency. The code is based on a finite difference discretization and memory-saving Runge-Kutta (RK) scheme. It handles non-ideal effects through super-time stepping and Hall diffusion schemes, and takes into account thermal conduction by solving an additional hyperbolic equation for the heat flux. The code is easily configurable to perform different kinds of simulations. Several examples of the code usage are given. It is demonstrated that splitting variables into equilibrium and perturbation parts is essential for simulations of wave propagation in a static background. A perfectly matched layer (PML) boundary condition built into the code greatly facilitates a non-reflective open boundary implementation. Spatial filtering is an important numerical remedy to eliminate grid-size perturbations enhancing the code stability. Parallel performance analysis reveals that the code is strongly memory bound, which is a natural consequence of the numerical techniques used, such as split variables and PML boundary conditions. Both strong and weak scalings show adequate performance up till several thousands of CPUs.
△ Less
Submitted 7 December, 2023;
originally announced December 2023.
-
The polarization of the Lyman-$α$ lines of H I and He II as a tool for exploring the solar corona
Authors:
Supriya Hebbur Dayananda,
Javier Trujillo Bueno,
Ángel de Vicente,
Tanausú del Pino Alemán
Abstract:
The near-Earth space weather is driven by the quick release of magnetic free energy in the solar corona. Probing this extremely hot and rarified region of the extended solar atmosphere requires modeling the polarization of forbidden and permitted coronal lines. To this end, it is important to develop efficient codes to calculate the Stokes profiles that emerge from given three-dimensional (3D) cor…
▽ More
The near-Earth space weather is driven by the quick release of magnetic free energy in the solar corona. Probing this extremely hot and rarified region of the extended solar atmosphere requires modeling the polarization of forbidden and permitted coronal lines. To this end, it is important to develop efficient codes to calculate the Stokes profiles that emerge from given three-dimensional (3D) coronal models, and this should be done taking into account the symmetry breaking produced by the presence of magnetic fields and non-radial solar wind velocities. We have developed such a tool with the aim of theoretically predicting and interpreting spectropolarimetric observations of the solar corona in permitted and forbidden lines. In this paper we show the results of a theoretical investigation of the linear polarization signals produced by scattering processes in the H I Ly-$α$ line at 1216 A and in the He II Ly-$α$ line at 304 A using 3D coronal models by Predictive Science Inc. These spectral lines have very different critical magnetic fields for the onset of the Hanle effect (53 G and 850 G, respectively), as well as different sensitivities to the Doppler effect caused by the solar wind velocities. We study under which circumstances simultaneous observations of the scattering polarization in these Ly-$α$ lines can facilitate the determination of magnetic fields and macroscopic velocities in the solar corona.
△ Less
Submitted 2 July, 2021;
originally announced July 2021.
-
Long characteristics vs. short characteristics in 3D radiative transfer simulations of polarized radiation
Authors:
A. de Vicente,
T. del Pino Alemán,
J. Trujillo Bueno
Abstract:
We compare maps of scattering polarization signals obtained from three-dimensional (3D) radiation transfer calculations in a magneto-convection model of the solar atmosphere using formal solvers based on the "short characteristics" (SC) and the "long characteristics" (LC) methods. The SC method requires less computational work, but it is known to introduce spatial blurring in the emergent radiatio…
▽ More
We compare maps of scattering polarization signals obtained from three-dimensional (3D) radiation transfer calculations in a magneto-convection model of the solar atmosphere using formal solvers based on the "short characteristics" (SC) and the "long characteristics" (LC) methods. The SC method requires less computational work, but it is known to introduce spatial blurring in the emergent radiation for inclined lines of sight. For polarized radiation this effect is generally more severe due to it being a signed quantity and to the sensitivity of the scattering polarization to the model's inhomogeneities. We study the differences in the polarization signals of the emergent spectral line radiation calculated with such formal solvers. We take as a case study already published results of the scattering polarization in the Sr I $4607~\unicode{xC5}$ line obtained with the SC method, demonstrating that in high-resolution grids it is accurate enough for that type of study. In general, the LC method is the preferred one for accurate calculations of the emergent radiation, reason why it is now one of the options in the public version of the 3D radiative transfer code PORTA.
△ Less
Submitted 6 March, 2021;
originally announced March 2021.
-
Two-fluid simulations of Rayleigh-Taylor instability in a magnetized solar prominence thread II. Effects of collisionality
Authors:
B. Popescu Braileanu,
V. S. Lukin,
E. Khomenko,
A. de Vicente
Abstract:
In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid m…
▽ More
In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid model takes into account neutral viscosity, thermal conductivity, and collisional interaction between neutrals and charges: ionization/recombination, energy and momentum transfer, and frictional heating. In this paper II, the sensitivity of the RTI dynamics to collisional effects for different magnetic field configurations supporting the prominence thread is explored. This is done by artificially varying, or eliminating, effects of both elastic and inelastic collisions by modifying the model equations. We find that ionization and recombination reactions between ionized and neutral fluids, if in equilibrium prior to the onset of the instability, do not substantially impact the development of the primary RTI. However, such reactions can impact development of secondary structures during mixing of the cold prominence and hotter surrounding coronal material. We find that collisionality within and between ionized and neutral particle populations play an important role in both linear and non-linear development of RTI, with ion-neutral collision frequency as the primary determining factor in development or damping of small scale structures. We also observe that degree and signatures of flow decoupling between ion and neutral fluids can depend both on the inter-particle collisionality and the magnetic field configuration of the prominence thread.
△ Less
Submitted 8 April, 2021; v1 submitted 29 January, 2021;
originally announced January 2021.
-
Two-fluid simulations of Rayleigh-Taylor instability in a magnetized solar prominence thread. I. Effects of prominence magnetization and mass loading
Authors:
B. Popescu Braileanu,
V. S. Lukin,
E. Khomenko,
A. de Vicente
Abstract:
Solar prominences are formed by partially ionized plasma with inter-particle collision frequencies generally warranting magnetohydrodynamic treatment. In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. We perform 2.5D two-fluid (charges - neutrals) sim…
▽ More
Solar prominences are formed by partially ionized plasma with inter-particle collision frequencies generally warranting magnetohydrodynamic treatment. In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. We perform 2.5D two-fluid (charges - neutrals) simulations of the Rayleigh-Taylor instability (RTI) at a smoothly changing interface between a solar prominence thread and the corona. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid model takes into account viscosity, thermal conductivity, and collisional interaction between neutrals and charges: ionization/recombination, energy and momentum transfer, and frictional heating. In this paper I, the sensitivity of the RTI dynamics to the prominence equilibrium configuration, including the impact of the magnetic field strength and shear supporting the prominence thread, and the amount of prominence mass-loading is explored. We show that, at small scales, a realistically smooth prominence-corona interface leads to qualitatively different linear RTI evolution than that expected for a discontinuous interface, while magnetic field shear has the stabilizing effect of reducing the growth rate or eliminating the instability. In the non-linear phase, we observe that in the presence of field shear the development of the instability leads to formation of coherent and interacting 2.5D magnetic structures, which, in turn, can lead to substantial plasma flow across magnetic field lines and associated decoupling of the fluid velocities of charges and neutrals.
△ Less
Submitted 15 December, 2020; v1 submitted 31 July, 2020;
originally announced July 2020.
-
Two-fluid simulations of waves in the solar chromosphere II. Propagation and damping of fast magneto-acoustic waves and shocks
Authors:
B. Popescu Braileanu,
V. S. Lukin,
E. Khomenko,
A. de Vicente
Abstract:
Waves and shocks traveling through the solar chromospheric plasma are influenced by its partial ionization and weak collisional coupling, and may become susceptible to multi-fluid effects, similar to interstellar shock waves. In this study, we consider fast magneto-acoustic shock wave formation and propagation in a stratified medium, that is permeated by a horizontal magnetic field, with propertie…
▽ More
Waves and shocks traveling through the solar chromospheric plasma are influenced by its partial ionization and weak collisional coupling, and may become susceptible to multi-fluid effects, similar to interstellar shock waves. In this study, we consider fast magneto-acoustic shock wave formation and propagation in a stratified medium, that is permeated by a horizontal magnetic field, with properties similar to that of the solar chromosphere. The evolution of plasma and neutrals is modeled using a two-fluid code that evolves a set of coupled equations for two separate fluids. We observed that waves in neutrals and plasma, initially coupled at the upper photosphere, become uncoupled at higher heights in the chromosphere. This decoupling can be a consequence of either the characteristic spatial scale at the shock front, that becomes similar to the collisional scale, or the change in the relation between the wave frequency, ion cyclotron frequency, and the collisional frequency with height. The decoupling height is a sensitive function of the wave frequency, wave amplitude, and the magnetic field strength. We observed that decoupling causes damping of waves and an increase in the background temperature due to the frictional heating. The comparison between analytical and numerical results allows us to separate the role of the nonlinear effects from the linear ones on the decoupling and damping of waves.
△ Less
Submitted 14 August, 2019;
originally announced August 2019.
-
Two-fluid simulations of waves in the solar chromosphere I: numerical code verification
Authors:
B. Popescu Braileanu,
V. S. Lukin,
E. Khomenko,
A. de Vicente
Abstract:
Solar chromosphere consists of a partially ionized plasma, which makes modeling the solar chromosphere a particularly challenging numerical task. Here we numerically model chromospheric waves using a two-fluid approach with a newly developed numerical code. The code solves two-fluid equations of conservation of mass, momentum and energy, together with the induction equation, for the case of the pu…
▽ More
Solar chromosphere consists of a partially ionized plasma, which makes modeling the solar chromosphere a particularly challenging numerical task. Here we numerically model chromospheric waves using a two-fluid approach with a newly developed numerical code. The code solves two-fluid equations of conservation of mass, momentum and energy, together with the induction equation, for the case of the purely hydrogen plasma with collisional coupling between the charged and neutral fluid components. The implementation of a semi-implicit algorithm allows us to overcome the numerical stability constraints due to the stiff collisional terms. We test the code against analytical solutions of acoustic and Alfvén wave propagation in uniform medium in several regimes of collisional coupling.The results of our simulations are consistent with the analytical estimates, and with other results described in the literature. In the limit of a large collisional frequency, the waves propagate with a common speed of a single fluid. In the other limit of a vanishingly small collisional frequency, the Alfvén waves propagate with an Alfvén speed of the charged fluid only, while the perturbation in neutral fluid is very small. The acoustic waves in these limits propagate with the sound speed corresponding to either the charges or the neutrals, while the perturbation in the other fluid component is very small. Otherwise, when the collision frequency is similar to the real part of the wave frequency, the interaction between charges and neutrals through momentum transfer collisions cause alterations of the waves frequencies and damping of the wave amplitudes.
△ Less
Submitted 9 May, 2019;
originally announced May 2019.
-
Three-dimensional simulations of solar magneto-convection including effects of partial ionization
Authors:
E. Khomenko,
N. Vitas,
M. Collados,
A. de Vicente
Abstract:
Over the last decades, realistic 3D radiative-MHD simulations have become the dominant theoretical tool for understanding the complex interactions between the plasma and the magnetic field on the Sun. Most of such simulations are based on approximations of magnetohydrodynamics, without directly considering the consequences of the very low degree of ionization of the solar plasma in the photosphere…
▽ More
Over the last decades, realistic 3D radiative-MHD simulations have become the dominant theoretical tool for understanding the complex interactions between the plasma and the magnetic field on the Sun. Most of such simulations are based on approximations of magnetohydrodynamics, without directly considering the consequences of the very low degree of ionization of the solar plasma in the photosphere and bottom chromosphere. The presence of large amount of neutrals leads to a partial decoupling of the plasma and the magnetic field. As a consequence of that, a series of non-ideal effects (ambipolar diffusion, Hall effect and battery effect) arises. The ambipolar effect is the dominant one in the solar chromosphere. Here we report on the first three-dimensional realistic simulations of magneto-convection including ambipolar diffusion and battery effects. The simulations are done using the newly developed Mancha3D code. Our results reveal that ambipolar diffusion causes measurable effects on the amplitudes of waves excited by convection in the simulations, on the absorption of Poynting flux and heating and on the formation of chromospheric structures. We provide a low limit on the chromospheric temperature increase due to the ambipolar effect using the simulations with battery-excited dynamo fields.
△ Less
Submitted 3 July, 2018;
originally announced July 2018.
-
MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma
Authors:
P. A. González-Morales,
E. Khomenko,
T. P. Downes,
A. de Vicente
Abstract:
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 extra terms in the generalised induction equation…
▽ More
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 extra terms in the generalised induction equation and in the energy conservation equation. These effects are: Ohmic diffusion, ambipolar diffusion, the Hall effect, and the Biermann battery effect. From the point of view of the numerical solution of the single-fluid equations, when ambipolar diffusion or Hall effects dominate can introduce severe restrictions on the integration time step and can compromise the stability of the numerical scheme. In this paper we introduce two numerical schemes to overcome those limitations. The first of them is known as Super Time-Stepping (STS) and it is designed to overcome the limitations imposed when the ambipolar diffusion term is dominant. The second scheme is called the Hall Diffusion Scheme (HDS) and it is used when the Hall term becomes dominant. These two numerical techniques can be used together by applying Strang operator splitting. This paper describes the implementation of the STS and HDS schemes in the single-fluid code Mancha3D. The validation for each of these schemes is provided by comparing the analytical solution with the numerical one for a suite of numerical tests.
△ Less
Submitted 13 March, 2018;
originally announced March 2018.
-
Numerical simulations of quiet Sun magnetic fields seeded by Biermann battery
Authors:
E. Khomenko,
N. Vitas,
M. Collados,
A. de Vicente
Abstract:
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 due to the action of local dynamo in the upper convection zone of the Sun. Existing simulations of the local solar dynamo require an initial seed field, and sufficien…
▽ More
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 due to the action of local dynamo in the upper convection zone of the Sun. Existing simulations of the local solar dynamo require an initial seed field, and sufficiently high spatial resolution, in order to achieve the amplification of the seed field to the observed values in the quiet Sun. Here we report an alternative model of seeding based on the action of the Bierman battery effect. This effect generates a magnetic field due to the local imbalances in electron pressure in the partially ionized solar plasma. We show that the battery effect self-consistently creates from zero an initial seed field of a strength of the order of micro G, and together with dynamo amplification, allows the generation of quiet Sun magnetic fields of a similar strength to those from solar observations.
△ Less
Submitted 19 June, 2017;
originally announced June 2017.
-
High frequency waves in the corona due to null points
Authors:
I. C. Santamaria,
E. Khomenko,
M. Collados,
A. de Vicente
Abstract:
This work aims to understand the behavior of non-linear waves in the vicinity of a coronal null point. In previous works we have showed that high frequency waves are generated in such magnetic configuration. This paper studies those waves in detail in order to provide a plausible explanation of their generation. We demonstrate that slow magneto-acoustic shock waves generated in the chromosphere pr…
▽ More
This work aims to understand the behavior of non-linear waves in the vicinity of a coronal null point. In previous works we have showed that high frequency waves are generated in such magnetic configuration. This paper studies those waves in detail in order to provide a plausible explanation of their generation. We demonstrate that slow magneto-acoustic shock waves generated in the chromosphere propagate through the null point and produce a train of secondary shocks that escape along the field lines. A particular combination of the shock wave speeds generates waves at a frequency of 80 mHz. We speculate that this frequency may be sensitive to the atmospheric parameters in the corona and therefore can be used to probe the structure of this solar layer.
△ Less
Submitted 21 April, 2017;
originally announced April 2017.
-
Simulated interaction of MHD shock waves with a complex network-like region
Authors:
Irantzu C. Santamaria,
Elena Khomenko,
Manuel Collados,
Angel de Vicente
Abstract:
We provide estimates of the wave energy reaching the solar chromosphere and corona in a network-like magnetic field topology, including a coronal null point. The waves are excited by an instantaneous strong subphotospheric source and propagate through the subphotosphere, photosphere, chromosphere, transition region, and corona with the plasma beta and other atmospheric parameters varying by severa…
▽ More
We provide estimates of the wave energy reaching the solar chromosphere and corona in a network-like magnetic field topology, including a coronal null point. The waves are excited by an instantaneous strong subphotospheric source and propagate through the subphotosphere, photosphere, chromosphere, transition region, and corona with the plasma beta and other atmospheric parameters varying by several orders of magnitude. We compare two regimes of the wave propagation: a linear and nonlinear regime. While the amount of energy reaching the corona is similar in both regimes, this energy is transmitted at different frequencies. In both cases the dominant periods of waves at each height strongly depend on the local magnetic field topology, but this distribution is only in accordance with observations in the nonlinear case.
△ Less
Submitted 29 April, 2016;
originally announced April 2016.
-
Heating of the partially ionized solar chromosphere by waves in magnetic structures
Authors:
S. Shelyag,
E. Khomenko,
A. de Vicente,
D. Przybylski
Abstract:
In this paper, we show a "proof of concept" of the heating mechanism of the solar chromosphere due to wave dissipation caused by the effects of partial ionization. Numerical modeling of non-linear wave propagation in a magnetic flux tube, embedded in the solar atmosphere, is performed by solving a system of single-fluid quasi-MHD equations, which take into account the ambipolar term from the gener…
▽ More
In this paper, we show a "proof of concept" of the heating mechanism of the solar chromosphere due to wave dissipation caused by the effects of partial ionization. Numerical modeling of non-linear wave propagation in a magnetic flux tube, embedded in the solar atmosphere, is performed by solving a system of single-fluid quasi-MHD equations, which take into account the ambipolar term from the generalized Ohm's law. It is shown that perturbations caused by magnetic waves can be effectively dissipated due to ambipolar diffusion. The energy input by this mechanism is continuous and shown to be more efficient than dissipation of static currents, ultimately leading to chromospheric temperature increase in magnetic structures.
△ Less
Submitted 10 February, 2016;
originally announced February 2016.
-
On the robustness of the pendulum model for large-amplitude longitudinal oscillations in prominences
Authors:
M. Luna,
J. Terradas,
E. Khomenko,
M. Collados,
A. de Vicente
Abstract:
Large-amplitude longitudinal oscillations (LALOs) in prominences are spectacular manifestations of the solar activity. In such events nearby energetic disturbances induce periodic motions on filaments with displacements comparable to the size of the filaments themselves and with velocities larger than 20 km/s. The pendulum model, in which the gravity projected along a rigid magnetic field is the r…
▽ More
Large-amplitude longitudinal oscillations (LALOs) in prominences are spectacular manifestations of the solar activity. In such events nearby energetic disturbances induce periodic motions on filaments with displacements comparable to the size of the filaments themselves and with velocities larger than 20 km/s. The pendulum model, in which the gravity projected along a rigid magnetic field is the restoring force, was proposed to explain these events. However, it can be objected that in a realistic situation where the magnetic field reacts to the mass motion of the heavy prominence, the simplified pendulum model could be no longer valid. We have performed non-linear time-dependent numerical simulations of LALOs considering a dipped magnetic field line structure. In this work we demonstrate that for even relatively weak magnetic fields the pendulum model works very well. We therefore validate the pendulum model and show its robustness, with important implications for prominence seismology purposes. With this model it is possible to infer the geometry of the dipped field lines that support the prominence.
△ Less
Submitted 16 December, 2015;
originally announced December 2015.
-
Rayleigh-Taylor instability in prominences from numerical simulations including partial ionization effects
Authors:
E. Khomenko,
A. Diaz,
A. de Vicente,
M. Collados,
M. Luna
Abstract:
We study the Rayleigh-Taylor instability (RTI) at a prominence-corona transition region in a non-linear regime. Our aim is to understand how the presence of neutral atoms in the prominence plasma influences the instability growth rate, and the evolution of velocity, magnetic field vector and thermodynamic parameters of turbulent drops. We perform 2.5D numerical simulations of the instability initi…
▽ More
We study the Rayleigh-Taylor instability (RTI) at a prominence-corona transition region in a non-linear regime. Our aim is to understand how the presence of neutral atoms in the prominence plasma influences the instability growth rate, and the evolution of velocity, magnetic field vector and thermodynamic parameters of turbulent drops. We perform 2.5D numerical simulations of the instability initiated by a multi-mode perturbation at the corona-prominence interface using a single-fluid MHD approach including a generalized Ohm's law. The initial equilibrium configuration is purely hydrostatic and contains a homogeneous horizontal magnetic field forming an angle with the direction in which the plasma is perturbed. We analyze simulations with two different orientations of the magnetic field. For each field orientation we compare two simulations, one for the pure MHD case, and one including the ambipolar diffusion in the Ohm's law (AD case). Other than that, both simulations for each field orientation are identical. The numerical results in the initial stage of the instability are compared with the analytical linear calculations. We find that the configuration is always unstable in the AD case. The growth rate of the small-scale modes in the non-linear regime is up to 50% larger in the AD case than in the purely MHD case and the average velocities of flows are a few percent larger. Significant drift momenta are found at the interface between the coronal and the prominence material at all stages of the instability, produced by the faster downward motion of the neutral component with respect to the ionized component. The differences in temperature of the bubbles between the ideal and non-ideal case are also significant, reaching 30%. There is an asymmetry between large rising bubbles and small-scale down flowing fingers, favoring the detection of upward velocities in observations.
△ Less
Submitted 18 March, 2014;
originally announced March 2014.
-
Rayleigh-Taylor instability in partially ionized prominence plasma
Authors:
E. Khomenko,
A. Diaz,
A. de Vicente,
M. Collados,
M. Luna
Abstract:
We study Rayleigh-Taylor instability (RTI) at the coronal-prominence boundary by means of 2.5D numerical simulations in a single-fluid MHD approach including a generalized Ohm's law. The initial configuration includes a homogeneous magnetic field forming an angle with the direction in which the plasma is perturbed. For each field inclination we compare two simulations, one for the pure MHD case, a…
▽ More
We study Rayleigh-Taylor instability (RTI) at the coronal-prominence boundary by means of 2.5D numerical simulations in a single-fluid MHD approach including a generalized Ohm's law. The initial configuration includes a homogeneous magnetic field forming an angle with the direction in which the plasma is perturbed. For each field inclination we compare two simulations, one for the pure MHD case, and one including the ambipolar diffusion in the Ohm's law, otherwise identical. We find that the configuration containing neutral atoms is always unstable. The growth rate of the small-scale modes in the non-linear regime is larger than in the purely MHD case.
△ Less
Submitted 28 October, 2013;
originally announced October 2013.
-
New ATLAS9 And MARCS Model Atmosphere Grids for the Apache Point Observatory Galactic Evolution Experiment (APOGEE)
Authors:
Sz. Mészáros,
C. Allende Prieto,
B. Edvardsson,
F. Castelli,
A. E. García Pérez,
B. Gustafsson,
S. R. Majewski,
B. Plez,
R. Schiavon,
M. Shetrone,
A. de Vicente
Abstract:
We present a new grid of model photospheres for the SDSS-III/APOGEE survey of stellar populations of the Galaxy, calculated using the ATLAS9 and MARCS codes. New opacity distribution functions were generated to calculate ATLAS9 model photospheres. MARCS models were calculated based on opacity sampling techniques. The metallicity ([M/H]) spans from -5 to 1.5 for ATLAS and -2.5 to 0.5 for MARCS mode…
▽ More
We present a new grid of model photospheres for the SDSS-III/APOGEE survey of stellar populations of the Galaxy, calculated using the ATLAS9 and MARCS codes. New opacity distribution functions were generated to calculate ATLAS9 model photospheres. MARCS models were calculated based on opacity sampling techniques. The metallicity ([M/H]) spans from -5 to 1.5 for ATLAS and -2.5 to 0.5 for MARCS models. There are three main differences with respect to previous ATLAS9 model grids: a new corrected H2O linelist, a wide range of carbon ([C/M]) and alpha element [alpha/M] variations, and solar reference abundances from Asplund et al. 2005. The added range of varying carbon and alpha element abundances also extends the previously calculated MARCS model grids. Altogether 1980 chemical compositions were used for the ATLAS9 grid, and 175 for the MARCS grid. Over 808 thousand ATLAS9 models were computed spanning temperatures from 3500K to 30000K and log g from 0 to 5, where larger temperatures only have high gravities. The MARCS models span from 3500K to 5500K, and log g from 0 to 5. All model atmospheres are publically available online.
△ Less
Submitted 9 August, 2012;
originally announced August 2012.
-
Automatic unsupervised classification of all SDSS/DR7 galaxy spectra
Authors:
J. Sanchez Almeida,
J. A. L. Aguerri,
C. Munoz-Tunon,
A. de Vicente
Abstract:
Using the 'k-means' cluster analysis algorithm, we carry out an unsupervised classification of all galaxy spectra in the seventh and final Sloan Digital Sky Survey data release (SDSS/DR7). Except for the shift to restframe wavelengths, and the normalization to the g-band flux, no manipulation is applied to the original spectra. The algorithm guarantees that galaxies with similar spectra belong to…
▽ More
Using the 'k-means' cluster analysis algorithm, we carry out an unsupervised classification of all galaxy spectra in the seventh and final Sloan Digital Sky Survey data release (SDSS/DR7). Except for the shift to restframe wavelengths, and the normalization to the g-band flux, no manipulation is applied to the original spectra. The algorithm guarantees that galaxies with similar spectra belong to the same class. We find that 99 % of the galaxies can be assigned to only 17 major classes, with 11 additional minor classes including the remaining 1%. The classification is not unique since many galaxies appear in between classes, however, our rendering of the algorithm overcomes this weakness with a tool to identify borderline galaxies. Each class is characterized by a template spectrum, which is the average of all the spectra of the galaxies in the class. These low noise template spectra vary smoothly and continuously along a sequence labeled from 0 to 27, from the reddest class to the bluest class. Our Automatic Spectroscopic K-means-based (ASK) classification separates galaxies in colors, with classes characteristic of the red sequence, the blue cloud, as well as the green valley. When red sequence galaxies and green valley galaxies present emission lines, they are characteristic of AGN activity. Blue galaxy classes have emission lines corresponding to star formation regions. We find the expected correlation between spectroscopic class and Hubble type, but this relationship exhibits a high intrinsic scatter. Several potential uses of the ASK classification are identified and sketched, including fast determination of physical properties by interpolation, classes as templates in redshift determinations, and target selection in follow-up works (we find classes of Seyfert galaxies, green valley galaxies, as well as a significant number of outliers). The ASK classification is publicly accessible through various websites.
△ Less
Submitted 16 March, 2010;
originally announced March 2010.
-
Characterization of horizontal flows around solar pores from high-resolution time series of images
Authors:
S. Vargas Dominguez,
A. de Vicente,
J. A. Bonet,
V. Martinez Pillet
Abstract:
Though there is increasing evidence linking the moat flow and the Evershed flow along the penumbral filaments, there is not a clear consensus regarding the existence of a moat flow around umbral cores and pores, and the debate is still open. Solar pores appear to be a suitable scenario to test the moat-penumbra relation as evidencing the direct interaction between the umbra and the convective plas…
▽ More
Though there is increasing evidence linking the moat flow and the Evershed flow along the penumbral filaments, there is not a clear consensus regarding the existence of a moat flow around umbral cores and pores, and the debate is still open. Solar pores appear to be a suitable scenario to test the moat-penumbra relation as evidencing the direct interaction between the umbra and the convective plasma in the surrounding photosphere, without any intermediate structure in between. The present work studies solar pores based on high resolution ground-based and satellite observations. Local correlation tracking techniques have been applied to different-duration time series to analyze the horizontal flows around several solar pores. Our results establish that the flows calculated from different solar pore observations are coherent among each other and show the determinant and overall influence of exploding events in the granulation around the pores. We do not find any sign of moat-like flows surrounding solar pores but a clearly defined region of inflows surrounding them. The connection between moat flows and flows associated to penumbral filaments is hereby reinforced by this work.
△ Less
Submitted 10 March, 2010;
originally announced March 2010.
-
Big Science with a Small Budget: Non-Embarrassingly Parallel Applications in a Non-Dedicated Network of Workstations
Authors:
Angel de Vicente,
Nayra Rodriguez
Abstract:
Many astronomers and astrophysicists require large computing resources for their research, which are usually obtained via dedicated (and expensive) parallel machines. Depending on the type of the problem to be solved, an alternative solution can be provided by creating dynamically a computer cluster out of non-dedicated workstations using the Condor High Throughput Computing System and the Maste…
▽ More
Many astronomers and astrophysicists require large computing resources for their research, which are usually obtained via dedicated (and expensive) parallel machines. Depending on the type of the problem to be solved, an alternative solution can be provided by creating dynamically a computer cluster out of non-dedicated workstations using the Condor High Throughput Computing System and the Master-Worker (MW) framework. As an example of this we show in this paper how a radiative transfer application previously coded with MPI is solved using this solution without the need for dedicated machines.
△ Less
Submitted 31 October, 2005;
originally announced October 2005.