-
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.
-
Influence of ambipolar and Hall effects on vorticity in 3D simulations of magneto-convection
Authors:
E. Khomenko,
M. Collados,
N. Vitas,
P. A. Gonzalez-Morales
Abstract:
This paper presents the results of the analysis of 3D simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation-runs are compared: one including both ambipolar diffusion and Hall effect; one including only ambipolar diffusion; and one without any of these two effects. The magnetic field is amplified from initial field to…
▽ More
This paper presents the results of the analysis of 3D simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation-runs are compared: one including both ambipolar diffusion and Hall effect; one including only ambipolar diffusion; and one without any of these two effects. The magnetic field is amplified from initial field to saturation level by the action of turbulent local dynamo. In each of these cases, we study 2 hours of simulated solar time after the local dynamo reaches the saturation regime. We analyze the power spectra of vorticity, of magnetic field fluctuations and of the different components of the magnetic Poynting flux responsible for the transport of vertical or horizontal perturbations. Our preliminary results show that the ambipolar diffusion produces a strong reduction of vorticity in the upper chromospheric layers and that it dissipates the vortical perturbations converting them into thermal energy. The Hall effect acts in the opposite way, strongly enhancing the vorticity. When the Hall effect is included, the magnetic field in the simulations becomes, on average, more vertical and long-lived flux tube-like structures are produced. We trace a single magnetic structure to study its evolution pattern and the magnetic field intensification, and their possible relation to the Hall effect.
△ Less
Submitted 21 September, 2020;
originally announced September 2020.
-
Joint action of Hall and ambipolar effects in 3D magneto-convection simulations of the quiet Sun. I. Dissipation and generation of waves
Authors:
P. A. Gonzalez-Morales,
E. Khomenko,
N. Vitas,
M. Collados
Abstract:
The partial ionization of the solar plasma causes several nonideal effects such as the ambipolar diffusion, the Hall effect, and the Biermann battery effect. Here we report on the first three-dimensional realistic simulations of solar local dynamo where all three effects were taken into account. The simulations started with a snapshot of already saturated battery-seeded dynamo, where two new serie…
▽ More
The partial ionization of the solar plasma causes several nonideal effects such as the ambipolar diffusion, the Hall effect, and the Biermann battery effect. Here we report on the first three-dimensional realistic simulations of solar local dynamo where all three effects were taken into account. The simulations started with a snapshot of already saturated battery-seeded dynamo, where two new series were developed: one with solely ambipolar diffusion and another one also taking into account the Hall term in the generalized Ohm's law. The simulations were then run for about 4 hours of solar time to reach the stationary regime and improve the statistics. In parallel, a purely MHD dynamo simulation was also run for the same amount of time. The simulations are compared in a statistical way. The results show that, with the inclusion of the ambipolar diffusion, the amplitudes of the incompressible perturbations related to Alfven waves are reduced, and the Poynting flux is absorbed, with a frequency dependence. The Hall effect causes the opposite action: significant excess of incompressible perturbations is generated and an excess of the Poynting flux is observed in the chromospheric layers. The model with ambipolar diffusion shows, on average, sharper current sheets and slightly more abundant fast magneto-acoustic shocks in the chromosphere. The model with the Hall effect has higher temperatures at the lower chromosphere and stronger and more vertical magnetic field concentrations all over the chromosphere. The study of high-frequency waves reveals that significant power of incompressible perturbations is associated with areas with intense and more vertical magnetic fields and larger temperatures. We find a positive correlation between the magnitude of the ambipolar heating and the temperature increase at the same location after a characteristic time of 10^2 sec.
△ Less
Submitted 24 August, 2020;
originally announced August 2020.
-
Fast-to-Alfvén mode conversion mediated by Hall current. II Application to the solar atmosphere
Authors:
P. A. González-Morales,
E. Khomenko,
P. S. Cally
Abstract:
Coupling between fast magneto-acoustic and Alfvén waves can be observe in fully ionized plasmas mediated by stratification and 3D geometrical effects. In Paper I, Cally & Khomenko (2015) have shown that in a weakly ionized plasma, such as the solar photosphere and chromosphere, the Hall current introduces a new coupling mechanism. The present study extends the results from Paper I to the case of w…
▽ More
Coupling between fast magneto-acoustic and Alfvén waves can be observe in fully ionized plasmas mediated by stratification and 3D geometrical effects. In Paper I, Cally & Khomenko (2015) have shown that in a weakly ionized plasma, such as the solar photosphere and chromosphere, the Hall current introduces a new coupling mechanism. The present study extends the results from Paper I to the case of warm plasma. We report on numerical experiments where mode transformation is studied using quasi-realistic stratification in thermodynamic parameters resembling the solar atmosphere. This redresses the limitation of the cold plasma approximation assumed in Paper I, in particular allowing the complete process of coupling between fast and slow magneto-acoustic modes and subsequent coupling of the fast mode to the Alfvén mode through the Hall current. Our results confirm the efficacy of the mechanism proposed in Paper I for the solar case. We observe that the efficiency of the transformation is a sensitive function of the angle between the wave propagation direction and the magnetic field, and of the wave frequency. The efficiency increases when the field direction and the wave direction are aligned for increasing wave frequencies. After scaling our results to typical solar values, the maximum amplitude of the transformed Alfvén waves, for a frequency of 1 Hz, corresponds to an energy flux (measured above the height of peak Hall coupling) of $\sim10^3$ $\rm W\,m^{-2}$, based on an amplitude of 500 $\rm m\,s^{-1}$ at $β=1$, which is sufficient to play a major role in both quiet and active region coronal heating.
△ Less
Submitted 15 November, 2018;
originally announced November 2018.
-
The extreme HBL behaviour of Markarian 501 during 2012
Authors:
M. L. Ahnen,
S. Ansoldi,
L. A. Antonelli,
C. Arcaro,
A. Babić,
B. Banerjee,
P. Bangale,
U. Barres de Almeida,
J. A. Barrio,
J. Becerra González,
W. Bednarek,
E. Bernardini,
A. Berti,
W. Bhattacharyya,
O. Blanch,
G. Bonnoli,
R. Carosi,
A. Carosi,
A. Chatterjee,
S. M. Colak,
P. Colin,
E. Colombo,
J. L. Contreras,
J. Cortina,
S. Covino
, et al. (254 additional authors not shown)
Abstract:
A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration.
Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0…
▽ More
A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration.
Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands.
The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.
△ Less
Submitted 14 August, 2018; v1 submitted 13 August, 2018;
originally announced August 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.
-
Exceptional outburst of the blazar CTA 102 in 2012: The GASP-WEBT campaign and its extension
Authors:
V. M. Larionov,
M. Villata,
C. M. Raiteri,
S. G. Jorstad,
A. P. Marscher,
I. Agudo,
P. S. Smith,
J. A. Acosta-Pulido,
M. J. Arévalo,
A. A. Arkharov,
R. Bachev,
D. A. Blinov,
G. Borisov,
G. A. Borman,
V. Bozhilov,
A. Bueno,
M. I. Carnerero,
D. Carosati,
C. Casadio,
W. P. Chen,
D. P. Clemens,
A. Di Paola,
Sh. A. Ehgamberdiev,
J. L. Gómez,
P. A González-Morales
, et al. (27 additional authors not shown)
Abstract:
After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished…
▽ More
After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004-2005. We find remarkable similarity between the optical and gamma-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ~1 hour, indicative of co-spatiality of the optical and gamma-ray emission regions. The relation between the gamma-ray and optical fluxes is consistent with the SSC mechanism, with a quadratic dependence of the SSC gamma-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the gamma-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in gamma-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV-NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet.
△ Less
Submitted 21 June, 2016;
originally announced June 2016.
-
Infrared properties of blazars: putting the GASP-WEBT sources into context
Authors:
C. M. Raiteri,
M. Villata,
M. I. Carnerero,
J. A. Acosta-Pulido,
V. M. Larionov,
F. D'Ammando,
M. J. Arévalo,
A. A. Arkharov,
A. Bueno Bueno,
A. Di Paola,
N. V. Efimova,
P. A. González-Morales,
D. L. Gorshanov,
A. B. Grinon-Marin,
C. Lázaro,
A. Manilla-Robles,
A. Pastor Yabar,
I. Puerto Giménez,
S. Velasco
Abstract:
The infrared properties of blazars can be studied from the statistical point of view with the help of sky surveys, like that provided by the Wide-field Infrared Survey Explorer (WISE) and the Two Micron All Sky Survey (2MASS). However, these sources are known for their strong and unpredictable variability, which can be monitored for a handful of objects only. In this paper we consider the 28 blaza…
▽ More
The infrared properties of blazars can be studied from the statistical point of view with the help of sky surveys, like that provided by the Wide-field Infrared Survey Explorer (WISE) and the Two Micron All Sky Survey (2MASS). However, these sources are known for their strong and unpredictable variability, which can be monitored for a handful of objects only. In this paper we consider the 28 blazars (14 BL Lac objects and 14 flat-spectrum radio quasars, FSRQs) that are regularly monitored by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) since 2007. They show a variety of infrared colours, redshifts, and infrared-optical spectral energy distributions (SEDs), and thus represent an interesting mini-sample of bright blazars that can be investigated in more detail. We present near-IR light curves and colours obtained by the GASP from 2007 to 2013, and discuss the infrared-optical SEDs. These are analysed with the aim of understanding the interplay among different emission components. BL Lac SEDs are accounted for by synchrotron emission plus an important contribution from the host galaxy in the closest objects, and dust signatures in 3C 66A and Mkn 421. FSRQ SEDs require synchrotron emission with the addition of a quasar-like contribution, which includes radiation from a generally bright accretion disc, broad line region, and a relatively weak dust torus.
△ Less
Submitted 16 May, 2014;
originally announced May 2014.