-
ViDA: a Vlasov-DArwin solver for plasma physics at electron scales
Authors:
Oreste Pezzi,
Giulia Cozzani,
Francesco Califano,
Francesco Valentini,
Massimiliano Guarrasi,
Enrico Camporeale,
Gianfranco Brunetti,
Alessandro Retinò,
Pierluigi Veltri
Abstract:
We present a Vlasov-DArwin numerical code (ViDA) specifically designed to address plasma physics problems, where small-scale high accuracy is requested even during the non linear regime to guarantee a clean description of the plasma dynamics at fine spatial scales. The algorithm provides a low-noise description of proton and electron kinetic dynamics, by splitting in time the multi-advection Vlaso…
▽ More
We present a Vlasov-DArwin numerical code (ViDA) specifically designed to address plasma physics problems, where small-scale high accuracy is requested even during the non linear regime to guarantee a clean description of the plasma dynamics at fine spatial scales. The algorithm provides a low-noise description of proton and electron kinetic dynamics, by splitting in time the multi-advection Vlasov equation in phase space. Maxwell equations for the electric and magnetic fields are reorganized according to Darwin approximation to remove light waves. Several numerical tests show that ViDA successfully reproduces the propagation of linear and nonlinear waves and captures the physics of magnetic reconnection. We also discuss preliminary tests of the parallelization algorithm efficiency, performed at CINECA on the Marconi-KNL cluster. ViDA will allow to run Eulerian simulations of a non-relativistic fully-kinetic collisionless plasma and it is expected to provide relevant insights on important problems of plasma astrophysics such as, for instance, the development of the turbulent cascade at electron scales and the structure and dynamics of electron-scale magnetic reconnection, such as the electron diffusion region.
△ Less
Submitted 24 October, 2019; v1 submitted 8 May, 2019;
originally announced May 2019.
-
3D MHD modeling of the expanding remnant of SN 1987A. Role of magnetic field and non-thermal radio emission
Authors:
S. Orlando,
M. Miceli,
O. Petruk,
M. Ono,
S. Nagataki,
M. A. Aloy,
P. Mimica,
S. -H. Lee,
F. Bocchino,
G. Peres,
M. Guarrasi
Abstract:
Aims. We investigate the role played by a pre-supernova (SN) ambient magnetic field on the dynamics of the expanding remnant of SN 1987A and the origin and evolution of the radio emission from the remnant, in particular, during the interaction of the blast wave with the nebula surrounding the SN. Methods. We model the evolution of SN 1987A from the breakout of the shock wave at the stellar surface…
▽ More
Aims. We investigate the role played by a pre-supernova (SN) ambient magnetic field on the dynamics of the expanding remnant of SN 1987A and the origin and evolution of the radio emission from the remnant, in particular, during the interaction of the blast wave with the nebula surrounding the SN. Methods. We model the evolution of SN 1987A from the breakout of the shock wave at the stellar surface to the expansion of its remnant through the surrounding nebula by 3D MHD simulations. The model considers the radiative cooling, the deviations from equilibrium of ionization, the deviation from temperature-equilibration between electrons and ions, and a plausible configuration of the pre-SN ambient magnetic field. From the simulations, we synthesize both thermal X-ray and non-thermal radio emission and compare the model results with observations. Results. The presence of an ambient magnetic field reduces the erosion and fragmentation of the dense equatorial ring after the impact of the SN blast wave. As a result, the ring survives the passage of the blast, at least, during the time covered by the simulations (40 years). Our model is able to reproduce the morphology and lightcurves of SN 1987A in both X-ray and radio bands. The model reproduces the observed radio emission if the flux originating from the reverse shock is heavily suppressed. In this case, the radio emission originates mostly from the forward shock traveling through the H II region. Possible mechanisms for the suppression of emission from the reverse shock are investigated. We find that synchrotron self-absorption and free-free absorption have negligible effects on the emission during the interaction with the nebula. We suggest that the emission from the reverse shock at radio frequencies might be limited by highly magnetized ejecta.
△ Less
Submitted 30 November, 2018;
originally announced December 2018.
-
Three-dimensional simulations of solar wind turbulence with the hybrid code CAMELIA
Authors:
L. Franci,
P. Hellinger,
M. Guarrasi,
C. H. K. Chen,
E. Papini,
A. Verdini,
L. Matteini,
S. Landi
Abstract:
We investigate the spectral properties of plasma turbulence from fluid to sub-ion scales by means of high-resolution three-dimensional (3D) numerical simulations performed with the hybrid particle-in-cell (HPIC) code CAMELIA. We produce extended turbulent spectra with well-defined power laws for the magnetic, ion bulk velocity, density, and electric fluctuations. The present results are in good ag…
▽ More
We investigate the spectral properties of plasma turbulence from fluid to sub-ion scales by means of high-resolution three-dimensional (3D) numerical simulations performed with the hybrid particle-in-cell (HPIC) code CAMELIA. We produce extended turbulent spectra with well-defined power laws for the magnetic, ion bulk velocity, density, and electric fluctuations. The present results are in good agreement with previous two-dimensional (2D) HPIC simulations, especially in the kinetic range of scales, and reproduce several features observed in solar wind spectra. By providing scaling tests on many different architectures and convergence studies, we prove CAMELIA to represent a very efficient, accurate and reliable tool for investigating the develpoment of the turbulent cascade in the solar wind, being able to cover simultaneously several decades in wavenumber, also in 3D.
△ Less
Submitted 11 December, 2017;
originally announced December 2017.
-
3D MHD modeling of twisted coronal loops
Authors:
F. Reale,
S. Orlando,
M. Guarrasi,
A. Mignone,
G. Peres,
A. W. Hood,
E. R. Priest
Abstract:
We perform MHD modeling of a single bright coronal loop to include the interaction with a non-uniform magnetic field. The field is stressed by random footpoint rotation in the central region and its energy is dissipated into heating by growing currents through anomalous magnetic diffusivity that switches on in the corona above a current density threshold. We model an entire single magnetic flux tu…
▽ More
We perform MHD modeling of a single bright coronal loop to include the interaction with a non-uniform magnetic field. The field is stressed by random footpoint rotation in the central region and its energy is dissipated into heating by growing currents through anomalous magnetic diffusivity that switches on in the corona above a current density threshold. We model an entire single magnetic flux tube, in the solar atmosphere extending from the high-beta chromosphere to the low-beta corona through the steep transition region. The magnetic field expands from the chromosphere to the corona. The maximum resolution is ~30 km. We obtain an overall evolution typical of loop models and realistic loop emission in the EUV and X-ray bands. The plasma confined in the flux tube is heated to active region temperatures (~3 MK) after ~2/3 hr. Upflows from the chromosphere up to ~100 km/s fill the core of the flux tube to densities above 10^9 cm^-3. More heating is released in the low corona than the high corona and is finely structured both in space and time.
△ Less
Submitted 19 July, 2016;
originally announced July 2016.
-
MHD modeling of coronal loops: the transition region throat
Authors:
M. Guarrasi,
F. Reale,
S. Orlando,
A. Mignone,
J. A. Klimchuk
Abstract:
The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in…
▽ More
The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in terms of steady heating models. We study the area response with a time-dependent 2D MHD loop model, including the description of the expanding magnetic field, coronal heating and losses by thermal conduction and radiation from optically thin plasma. We run a simulation for a loop 50 Mm long and quasi-statically heated to about 4 MK. We find that the area can change substantially with the quasi-steady heating rate, e.g. by ~40% at 0.5 MK as the loop temperature varies between 1 and 4 MK, and, therefore, affects the interpretation of DEM(T) curves.
△ Less
Submitted 3 February, 2014;
originally announced February 2014.
-
Solar Dynamics Observatory discovers thin high temperature strands in coronal active regions
Authors:
Fabio Reale,
Massimiliano Guarrasi,
Paola Testa,
Edward E. DeLuca,
Giovanni Peres,
Leon Golub
Abstract:
One scenario proposed to explain the million degrees solar corona is a finely-stranded corona where each strand is heated by a rapid pulse. However, such fine structure has neither been resolved through direct imaging observations nor conclusively shown through indirect observations of extended superhot plasma. Recently it has been shown that the observed difference in appearance of cool and warm…
▽ More
One scenario proposed to explain the million degrees solar corona is a finely-stranded corona where each strand is heated by a rapid pulse. However, such fine structure has neither been resolved through direct imaging observations nor conclusively shown through indirect observations of extended superhot plasma. Recently it has been shown that the observed difference in appearance of cool and warm coronal loops (~1 MK, ~2-3 MK, respectively) -- warm loops appearing "fuzzier" than cool loops -- can be explained by models of loops composed of subarcsecond strands, which are impulsively heated up to ~10 MK. That work predicts that images of hot coronal loops (>~6 MK) should again show fine structure. Here we show that the predicted effect is indeed widely observed in an active region with the Solar Dynamics Observatory, thus supporting a scenario where impulsive heating of fine loop strands plays an important role in powering the active corona.
△ Less
Submitted 8 June, 2011;
originally announced June 2011.
-
Coronal fuzziness modelled with pulse-heated multistranded loop systems
Authors:
Massimiliano Guarrasi,
Fabio Reale,
Giovanni Peres
Abstract:
Coronal active regions are observed to get fuzzier and fuzzier (i.e. more and more confused and uniform) in harder and harder energy bands or lines. We explain this evidence as due to the fine multi-temperature structure of coronal loops. To this end, we model bundles of loops made of thin strands, each heated by short and intense heat pulses. For simplicity, we assume that the heat pulses are all…
▽ More
Coronal active regions are observed to get fuzzier and fuzzier (i.e. more and more confused and uniform) in harder and harder energy bands or lines. We explain this evidence as due to the fine multi-temperature structure of coronal loops. To this end, we model bundles of loops made of thin strands, each heated by short and intense heat pulses. For simplicity, we assume that the heat pulses are all equal and triggered only once in each strand at a random time. The pulse intensity and cadence are selected so as to have steady active region loops ($\sim 3$ MK), on the average. We compute the evolution of the confined heated plasma with a hydrodynamic loop model. We then compute the emission along each strand in several spectral lines, from cool ($\leq 1$ MK), to warm ($2-3$ MK) lines, detectable with Hinode/EIS, to hot X-ray lines. The strands are then put side-by-side to construct an active region loop bundle. We find that in the warm lines ($2-3$ MK) the loop emission fills all the available image surface. Therefore the emission appears quite uniform and it is difficult to resolve the single loops, while in the cool lines the loops are considerably more contrasted and the region is less fuzzy. The main reasons for this effect are that, during their evolution, i.e. pulse heating and slow cooling, each strand spends a relatively long time at temperatures around $2-3$ MK, and that it has a high emission measure during that phase, so the whole region appears more uniform or smudged. We make the prediction that the fuzziness should be reduced in the hot UV and X-ray lines.
△ Less
Submitted 17 June, 2010;
originally announced June 2010.