-
Properties of an interplanetary shock observed at 0.07 and 0.7 Astronomical Units by Parker Solar Probe and Solar Orbiter
Authors:
D. Trotta,
A. Larosa,
G. Nicolaou,
T. S. Horbury,
L. Matteini,
H. Hietala,
X. Blanco-Cano,
L. Franci,
C. H. K. Chen,
L. Zhao,
G. P. Zank,
C. M. S. Cohen,
S. D. Bale,
R. Laker,
N. Fargette,
F. Valentini,
Y. Khotyaintsev,
R. Kieokaew,
N. Raouafi,
E. Davies,
R. Vainio,
N. Dresing,
E. Kilpua,
T. Karlsson,
C. J. Owen
, et al. (1 additional authors not shown)
Abstract:
The Parker Solar Probe (PSP) and Solar Orbiter (SolO) missions opened a new observational window in the inner heliosphere, which is finally accessible to direct measurements. On September 05, 2022, a coronal mass ejection (CME)-driven interplanetary (IP) shock has been observed as close as 0.07 au by PSP. The CME then reached SolO, which was well radially-aligned at 0.7 au, thus providing us with…
▽ More
The Parker Solar Probe (PSP) and Solar Orbiter (SolO) missions opened a new observational window in the inner heliosphere, which is finally accessible to direct measurements. On September 05, 2022, a coronal mass ejection (CME)-driven interplanetary (IP) shock has been observed as close as 0.07 au by PSP. The CME then reached SolO, which was well radially-aligned at 0.7 au, thus providing us with the opportunity to study the shock properties at so different heliocentric distances. We characterize the shock, investigate its typical parameters and compare its small-scale features at both locations. Using the PSP observations, we investigate how magnetic switchbacks and ion cyclotron waves are processed upon shock crossing. We find that switchbacks preserve their V--B correlation while compressed upon the shock passage, and that the signature of ion cyclotron waves disappears downstream of the shock. By contrast, the SolO observations reveal a very structured shock transition, with a population of shock-accelerated protons of up to about 2 MeV, showing irregularities in the shock downstream, which we correlate with solar wind structures propagating across the shock. At SolO, we also report the presence of low-energy ($\sim$ 100 eV) electrons scattering due to upstream shocklets. This study elucidates how the local features of IP shocks and their environments can be very different as they propagate through the heliosphere.
△ Less
Submitted 10 December, 2023;
originally announced December 2023.
-
The essential role of multi-point measurements in investigations of turbulence, three-dimensional structure, and dynamics: the solar wind beyond single scale and the Taylor Hypothesis
Authors:
W. H. Matthaeus,
S. Adhikari,
R. Bandyopadhyay,
M. R. Brown,
R. Bruno,
J. Borovsky,
V. Carbone,
D. Caprioli,
A. Chasapis,
R. Chhiber,
S. Dasso,
P. Dmitruk,
L. Del Zanna,
P. A. Dmitruk,
Luca Franci,
S. P. Gary,
M. L. Goldstein,
D. Gomez,
A. Greco,
T. S. Horbury,
Hantao Ji,
J. C. Kasper,
K. G. Klein,
S. Landi,
Hui Li
, et al. (27 additional authors not shown)
Abstract:
Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer h…
▽ More
Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer heliosphere. This shortcoming leaves numerous central questions incompletely answered. Deficiencies remain in at least two important subjects, Space Weather and fundamental plasma turbulence theory, due to a lack of a more complete understanding of the space-time structure of dynamic plasmas. Only with multispacecraft measurements over suitable spans of spatial separation and temporal duration can these ambiguities be resolved. We note that these characterizations apply to turbulence across a wide range of scales, and also equally well to shocks, flux ropes, magnetic clouds, current sheets, stream interactions, etc. In the following, we will describe the basic requirements for resolving space-time structure in general, using turbulence' as both an example and a principal target or study. Several types of missions are suggested to resolve space-time structure throughout the Heliosphere.
△ Less
Submitted 26 November, 2022; v1 submitted 22 November, 2022;
originally announced November 2022.
-
Ion-scale transition of plasma turbulence: Pressure-strain effect
Authors:
Petr Hellinger,
Victor Montagud-Camps,
Luca Franci,
Lorenzo Matteini,
Emanuele Papini,
Andrea Verdini,
Simone Landi
Abstract:
We investigate properties of solar wind-like plasma turbulence using direct numerical simulations. We analyze the transition from large, magnetohydrodynamic (MHD) scales to the ion characteristic ones using two-dimensional hybrid (fluid electrons, kinetic ions) simulations. To capture and quantify turbulence properties, we apply the Karman-Howarth-Monin (KHM) equation for compressible Hall MHD (ex…
▽ More
We investigate properties of solar wind-like plasma turbulence using direct numerical simulations. We analyze the transition from large, magnetohydrodynamic (MHD) scales to the ion characteristic ones using two-dimensional hybrid (fluid electrons, kinetic ions) simulations. To capture and quantify turbulence properties, we apply the Karman-Howarth-Monin (KHM) equation for compressible Hall MHD (extended by considering the plasma pressure as a tensor quantity) to the numerical results. The KHM analysis indicates that the transition from MHD to ion scales (the so called ion break in the power spectrum) results from a combination of an onset of Hall physics and of an effective dissipation owing to the pressure-strain energy-exchange channel and resistivity. We discuss the simulation results in the context of the solar wind.
△ Less
Submitted 23 March, 2022;
originally announced March 2022.
-
Spectral transfer and Kármán-Howarth-Monin equations for compressible Hall magnetohydrodynamics
Authors:
Petr Hellinger,
Emanuele Papini,
Andrea Verdini,
Simone Landi,
Luca Franci,
Lorenzo Matteini,
Victor Montagud-Camps
Abstract:
We derive two new forms of the Kármán-Howarth-Monin equation for decaying compressible Hall magnetohydrodynamic (MHD) turbulence. We test them on results of a weakly-compressible, two-dimensional, moderate-Reynolds-number Hall MHD simulation and compare them with an isotropic spectral transfer (ST) equation. The KHM and ST equations are automatically satisfied during the whole simulation owing to…
▽ More
We derive two new forms of the Kármán-Howarth-Monin equation for decaying compressible Hall magnetohydrodynamic (MHD) turbulence. We test them on results of a weakly-compressible, two-dimensional, moderate-Reynolds-number Hall MHD simulation and compare them with an isotropic spectral transfer (ST) equation. The KHM and ST equations are automatically satisfied during the whole simulation owing to the periodic boundary conditions and have complementary cumulative behavior. They are used here to analyze the onset of turbulence and its properties when it is fully developed. These approaches give equivalent results characterizing: the decay of the kinetic + magnetic energy at large scales, the MHD and Hall cross-scale energy transfer/cascade, the pressure dilatation, and the dissipation. The Hall cascade appears when the MHD one brings the energy close to the ion inertial range and is related to the formation of reconnecting current sheets. At later times, the pressure-dilation energy-exchange rate oscillates around zero with no net effect on the cross-scale energy transfer when averaged over a period of its oscillations. A reduced one-dimensional analysis suggests that all three methods may be useful to estimate the energy cascade rate from in situ observations.
△ Less
Submitted 14 April, 2021;
originally announced April 2021.
-
Modeling Kelvin-Helmholtz instability-driven turbulence with hybrid simulations of Alfvénic turbulence
Authors:
Luca Franci,
Julia E. Stawarz,
Emanuele Papini,
Petr Hellinger,
Takuma Nakamura,
David Burgess,
Simone Landi,
Andrea Verdini,
Lorenzo Matteini,
Robert Ergun,
Olivier Le Contel,
Per-Arne Lindqvist
Abstract:
Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfvénic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and…
▽ More
Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfvénic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and electron betas, turbulence strength, and injection scale) exhibits a quantitative agreement on the spectral, intermittency, and cascade-rate properties with in situ observations, despite the different driving mechanisms. Such agreement demonstrates a certain universality of the turbulent cascade from magnetohydrodynamic (MHD) to sub-ion scales, whose properties are mainly determined by the selected parameters, also indicating that the KH instability-driven turbulence has a quasi-2D nature. The validity of the Taylor hypothesis in the sub-ion spatial range suggests that the fluctuations at sub-ion scales have predominantly low frequencies, consistent with a kinetic Alfvén wave-like nature or with quasi-static structures. Finally, the third-order structure function analysis indicates that the cascade rate of the turbulence generated by a KH event in the magnetopause is an order of magnitude larger than in the ambient magnetosheath.
△ Less
Submitted 17 November, 2019;
originally announced November 2019.
-
Fast acceleration of transrelativistic electrons in astrophysical turbulence
Authors:
Domenico Trotta,
Luca Franci,
David Burgess,
Petr Hellinger
Abstract:
Highly energetic, relativistic electrons are commonly present in many astrophysical systems, from solar flares to the intra-cluster medium, as indicated by observed electromagnetic radiation. However, open questions remain about the mechanisms responsible for their acceleration, and possible re-acceleration. Ubiquitous plasma turbulence is one of the possible universal mechanisms. We study the ene…
▽ More
Highly energetic, relativistic electrons are commonly present in many astrophysical systems, from solar flares to the intra-cluster medium, as indicated by observed electromagnetic radiation. However, open questions remain about the mechanisms responsible for their acceleration, and possible re-acceleration. Ubiquitous plasma turbulence is one of the possible universal mechanisms. We study the energization of transrelativistic electrons in turbulence using hybrid particle-in-cell, which provide a realistic model of Alfvénic turbulence from MHD to sub-ion scales, and test particle simulations for electrons. We find that, depending on the electron initial energy and turbulence strength, electrons may undergo a fast and efficient phase of energization due to the magnetic curvature drift during the time they are trapped in dynamic magnetic structures. In addition, electrons are accelerated stochastically which is a slower process that yields lower maximum energies. The combined effect of these two processes determines the overall electron acceleration. With appropriate turbulence parameters, we find that superthermal electrons can be accelerated up to relativistic energies. For example, with heliospheric parameters and a relatively high turbulence level, rapid energization to MeV energies is possible.
△ Less
Submitted 17 April, 2020; v1 submitted 21 October, 2019;
originally announced October 2019.
-
Kinetic plasma turbulence: recent insights and open questions from 3D3V simulations
Authors:
S. S. Cerri,
D. Grošelj,
L. Franci
Abstract:
Turbulence and kinetic processes in magnetized space plasmas have been extensively investigated over the past decades via \emph{in-situ} spacecraft measurements, theoretical models and numerical simulations. In particular, multi-point high-resolution measurements from the \emph{Cluster} and \emph{MMS} space missions brought to light an entire new world of processes, taking place at the plasma kine…
▽ More
Turbulence and kinetic processes in magnetized space plasmas have been extensively investigated over the past decades via \emph{in-situ} spacecraft measurements, theoretical models and numerical simulations. In particular, multi-point high-resolution measurements from the \emph{Cluster} and \emph{MMS} space missions brought to light an entire new world of processes, taking place at the plasma kinetic scales, and exposed new challenges for their theoretical interpretation. A long-lasting debate concerns the nature of ion and electron scale fluctuations in solar-wind turbulence and their dissipation via collisionless plasma mechanisms. Alongside observations, numerical simulations have always played a central role in providing a test ground for existing theories and models. In this Perspective, we discuss the advances achieved with our 3D3V (reduced and fully) kinetic simulations, as well as the main questions left open (or raised) by these studies. To this end, we combine data from our recent kinetic simulations of both freely decaying and continuously driven fluctuations to assess the similarities and/or differences in the properties of plasma turbulence in the sub-ion range. Finally, we discuss possible future directions in the field and highlight the need to combine different types of numerical and observational approaches to improve the understanding of turbulent space plasmas.
△ Less
Submitted 11 October, 2019; v1 submitted 25 September, 2019;
originally announced September 2019.
-
Turbulence vs. fire hose instabilities: 3-D hybrid expanding box simulations
Authors:
P. Hellinger,
L. Matteini,
S. Landi,
L. Franci,
A. Verdini,
E. Papini
Abstract:
The relationship between a decaying plasma turbulence and proton fire hose instabilities in a slowly expanding plasma is investigated using three-dimensional (3-D) hybrid expanding box simulations. We impose an initial ambient magnetic field along the radial direction, and we start with an isotropic spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations with zero cross-hel…
▽ More
The relationship between a decaying plasma turbulence and proton fire hose instabilities in a slowly expanding plasma is investigated using three-dimensional (3-D) hybrid expanding box simulations. We impose an initial ambient magnetic field along the radial direction, and we start with an isotropic spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations with zero cross-helicity. A turbulent cascade rapidly develops and leads to a weak proton heating that is not sufficient to overcome the expansion-driven perpendicular cooling. The plasma system eventually drives the parallel and oblique fire hose instabilities that generate quasi-monochromatic wave packets that reduce the proton temperature anisotropy. The fire hose wave activity has a low amplitude with wave vectors quasi-parallel/oblique with respect to the ambient magnetic field outside of the region dominated by the turbulent cascade and is discernible in one-dimensional power spectra taken only in the direction quasi-parallel/oblique with respect to the ambient magnetic field; at quasi-perpendicular angles the wave activity is hidden by the turbulent background. These waves are partly reabsorbed by protons and partly couple to and participate in the turbulent cascade. Their presence reduces kurtosis, a measure of intermittency, and the Shannon entropy but increases the Jensen-Shannon complexity of magnetic fluctuations; these changes are weak and anisotropic with respect to the ambient magnetic field and it's not clear if they can be used to indirectly discern the presence of instability-driven waves.
△ Less
Submitted 21 August, 2019;
originally announced August 2019.
-
A Case for Electron-Astrophysics
Authors:
Daniel Verscharen,
Robert T. Wicks,
Olga Alexandrova,
Roberto Bruno,
David Burgess,
Christopher H. K. Chen,
Raffaella D'Amicis,
Johan De Keyser,
Thierry Dudok de Wit,
Luca Franci,
Jiansen He,
Pierre Henri,
Satoshi Kasahara,
Yuri Khotyaintsev,
Kristopher G. Klein,
Benoit Lavraud,
Bennett A. Maruca,
Milan Maksimovic,
Ferdinand Plaschke,
Stefaan Poedts,
Chirstopher S. Reynolds,
Owen Roberts,
Fouad Sahraoui,
Shinji Saito,
Chadi S. Salem
, et al. (5 additional authors not shown)
Abstract:
A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour,…
▽ More
A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.
△ Less
Submitted 6 August, 2019;
originally announced August 2019.
-
Three-dimensional local anisotropy of velocity fluctuations in the solar wind
Authors:
Andrea Verdini,
Roland Grappin,
Olga Alexandrova,
Luca Franci,
Simone Landi,
Lorenzo Matteini,
Emanuele Papini
Abstract:
We analyse velocity fluctuations in the solar wind at magneto-fluid scales in two datasets, extracted from Wind data in the period 2005-2015, that are characterised by strong or weak expansion. Expansion affects measurements of anisotropy because it breaks axisymmetry around the mean magnetic field. Indeed, the small-scale three-dimensional local anisotropy of magnetic fluctuations (δB) as measure…
▽ More
We analyse velocity fluctuations in the solar wind at magneto-fluid scales in two datasets, extracted from Wind data in the period 2005-2015, that are characterised by strong or weak expansion. Expansion affects measurements of anisotropy because it breaks axisymmetry around the mean magnetic field. Indeed, the small-scale three-dimensional local anisotropy of magnetic fluctuations (δB) as measured by structure functions (SF_B) is consistent with tube-like structures for strong expansion. When passing to weak expansion, structures become ribbon-like because of the flattening of SFB along one of the two perpendicular directions. The power-law index that is consistent with a spectral slope -5/3 for strong expansion now becomes closer to -3/2. This index is also characteristic of velocity fluctuations in the solar wind. We study velocity fluctuations (δV) to understand if the anisotropy of their structure functions (SF_V ) also changes with the strength of expansion and if the difference with the magnetic spectral index is washed out once anisotropy is accounted for. We find that SF_V is generally flatter than SF_B. When expansion passes from strong to weak, a further flattening of the perpendicular SF_V occurs and the small-scale anisotropy switches from tube-like to ribbon-like structures. These two types of anisotropy, common to SF_V and SF_B, are associated to distinct large-scale variance anisotropies of δB in the strong- and weak-expansion datasets. We conclude that SF_V shows anisotropic three-dimensional scaling similar to SF_B, with however systematic flatter scalings, reflecting the difference between global spectral slopes.
△ Less
Submitted 8 April, 2019;
originally announced April 2019.
-
Spectral anisotropies and intermittency of plasma turbulence at ion kinetic scales
Authors:
Simone Landi,
Luca Franci,
Emanuele Papini,
Andrea Verdini,
Lorenzo Matteini,
Petr Hellinger
Abstract:
By means of three dimensional high-resolution hybrid simulations we study the properties of the magnetic field spectral anisotropies near and beyond ion kinetic scales. By using both a Fourier analysis and a local analysis based on multi-point 2nd-order structure function techniques, we show that the anisotropy observed is less than what expected by standard wave normal modes turbulence theories a…
▽ More
By means of three dimensional high-resolution hybrid simulations we study the properties of the magnetic field spectral anisotropies near and beyond ion kinetic scales. By using both a Fourier analysis and a local analysis based on multi-point 2nd-order structure function techniques, we show that the anisotropy observed is less than what expected by standard wave normal modes turbulence theories although the non linear energy transfer is still in the perpendicular direction, only advected in the parallel direction as expected balancing the non-linear energy transfer time and the decorrelation time. Such result can be explained by a phenomenological model based on the formation of strong intermittent two-dimensional structures in the plane perpendicular to the local mean field that have some prescribed aspect ratio eventually depending on the scale. This model support the idea that small scales structures, such as reconnecting current sheets, contribute significantly to the formation of the turbulent cascade at kinetic scales.
△ Less
Submitted 17 April, 2019; v1 submitted 8 April, 2019;
originally announced April 2019.
-
[Plasma 2020 Decadal] The essential role of multi-point measurements in turbulence investigations: the solar wind beyond single scale and beyond the Taylor Hypothesis
Authors:
W. H. Matthaeus,
R. Bandyopadhyay,
M. R. Brown,
J. Borovsky,
V. Carbone,
D. Caprioli,
A. Chasapis,
R. Chhiber,
S. Dasso,
P. Dmitruk,
L. Del Zanna,
P. A. Dmitruk,
Luca Franci,
S. P. Gary,
M. L. Goldstein,
D. Gomez,
A. Greco,
T. S. Horbury,
Hantao Ji,
J. C. Kasper,
K. G. Klein,
S. Landi,
Hui Li,
F. Malara,
B. A. Maruca
, et al. (24 additional authors not shown)
Abstract:
This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for su…
▽ More
This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.
△ Less
Submitted 16 March, 2019;
originally announced March 2019.
-
Von Karman-Howarth equation for Hall magnetohydrodynamics: Hybrid simulations
Authors:
Petr Hellinger,
Andrea Verdini,
Simone Landi,
Luca Franci,
Lorenzo Matteini
Abstract:
A dynamical vectorial equation for homogeneous incompressible Hall-MHD turbulence together with the exact scaling law for third-order correlation tensors, analogous to that for the incompressible MHD, is rederived and applied to the results of two-dimensional hybrid simulations of plasma turbulence. At large (MHD) scales the simulations exhibits a clear inertial range where the MHD dynamic law is…
▽ More
A dynamical vectorial equation for homogeneous incompressible Hall-MHD turbulence together with the exact scaling law for third-order correlation tensors, analogous to that for the incompressible MHD, is rederived and applied to the results of two-dimensional hybrid simulations of plasma turbulence. At large (MHD) scales the simulations exhibits a clear inertial range where the MHD dynamic law is valid. In the sub-ion range the cascade continues via the Hall term but the dynamic law derived in the framework of incompressible Hall MHD equations is obtained only in a low plasma beta simulation. For a higher beta plasma the cascade rate decreases in the sub-ion range and the change becomes more pronounced as the plasma beta increases. This break in the cascade flux can be ascribed to non thermal (kinetic) features or to others terms in the dynamical equation that are not included in the Hall-MHD incompressible approximation.
△ Less
Submitted 26 March, 2018;
originally announced March 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.
-
Solar wind turbulent cascade from MHD to sub-ion scales: large-size 3D hybrid particle-in-cell simulations
Authors:
Luca Franci,
Simone Landi,
Andrea Verdini,
Lorenzo Matteini,
Petr Hellinger
Abstract:
Spectral properties of the turbulent cascade from fluid to kinetic scales in collisionless plasmas are investigated by means of large-size three-dimensional (3D) hybrid (fluid electrons, kinetic protons) particle-in-cell simulations. Initially isotropic Alfvènic fluctuations rapidly develop a strongly anisotropic turbulent cascade, mainly in the direction perpendicular to the ambient magnetic fiel…
▽ More
Spectral properties of the turbulent cascade from fluid to kinetic scales in collisionless plasmas are investigated by means of large-size three-dimensional (3D) hybrid (fluid electrons, kinetic protons) particle-in-cell simulations. Initially isotropic Alfvènic fluctuations rapidly develop a strongly anisotropic turbulent cascade, mainly in the direction perpendicular to the ambient magnetic field. The omnidirectional magnetic field spectrum shows a double power-law behavior over almost two decades in wavenumber, with a Kolmogorov-like index at large scales, a spectral break around ion scales, and a steepening at sub-ion scales. Power laws are also observed in the spectra of the ion bulk velocity, density, and electric field, both at magnetohydrodynamic (MHD) and at kinetic scales. Despite the complex structure, the omnidirectional spectra of all fields at ion and sub-ion scales are in remarkable quantitative agreement with those of a two-dimensional (2D) simulation with similar physical parameters. This provides a partial, a-posteriori validation of the 2D approximation at kinetic scales. Conversely, at MHD scales, the spectra of the density and of the velocity (and, consequently, of the electric field) exhibit differences between the 2D and 3D cases. Although they can be partly ascribed to the lower spatial resolution, the main reason is likely the larger importance of compressible effects in a full geometry. Our findings are also in remarkable quantitative agreement with solar wind observations.
△ Less
Submitted 6 November, 2017;
originally announced November 2017.
-
Magnetic reconnection as a driver for a sub-ion scale cascade in plasma turbulence
Authors:
L. Franci,
S. S. Cerri,
F. Califano,
S. Landi,
E. Papini,
A. Verdini,
L. Matteini,
F. Jenko,
P. Hellinger
Abstract:
A new path for the generation of a sub-ion scale cascade in collisionless plasma turbulence, triggered by magnetic reconnection, is uncovered by means of high-resolution two-dimensional hybrid-kinetic simulations employing two complementary approaches, Lagrangian and Eulerian, and different driving mechanisms. The simulation results provide clear numerical evidences that the development of power-l…
▽ More
A new path for the generation of a sub-ion scale cascade in collisionless plasma turbulence, triggered by magnetic reconnection, is uncovered by means of high-resolution two-dimensional hybrid-kinetic simulations employing two complementary approaches, Lagrangian and Eulerian, and different driving mechanisms. The simulation results provide clear numerical evidences that the development of power-law energy spectra below the so-called ion break occurs as soon as the first magnetic reconnection events take place, regardless of the actual state of the turbulent cascade at MHD scales. In both simulations, the reconnection-mediated small-scale energy spectrum of parallel magnetic fluctuations exhibits a very stable spectral slope of ~-2.8, whether or not a large-scale turbulent cascade has already fully developed. Once a quasi-stationary turbulent state is achieved, the spectrum of the total magnetic fluctuations settles towards a spectral index of -5/3 in the MHD range and of ~-3 at sub-ion scales.
△ Less
Submitted 19 July, 2017;
originally announced July 2017.
-
Mirror instability in the turbulent solar wind
Authors:
P. Hellinger,
S. Landi,
L. Matteini,
A. Verdini,
L. Franci
Abstract:
The relationship between a decaying strong turbulence and the mirror instability in a slowly expanding plasma is investigated using two-dimensional hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations which have energy equipartition betwe…
▽ More
The relationship between a decaying strong turbulence and the mirror instability in a slowly expanding plasma is investigated using two-dimensional hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at sub-ion scales. The imposed expansion (taking a strictly transverse ambient magnetic field) leads to generation of an important perpendicular proton temperature anisotropy that eventually drives the mirror instability. This instability generates large-amplitude, nonpropagating, compressible, pressure-balanced magnetic structures in a form of magnetic enhancements/humps that reduce the perpendicular temperature anisotropy.
△ Less
Submitted 17 March, 2017;
originally announced March 2017.
-
Plasma turbulence at ion scales: a comparison between PIC and Eulerian hybrid-kinetic approaches
Authors:
S. S. Cerri,
L. Franci,
F. Califano,
S. Landi,
P. Hellinger
Abstract:
Kinetic-range turbulence in magnetized plasmas and, in particular, in the context of solar-wind turbulence has been extensively investigated over the past decades via numerical simulations. Among others, one of the widely adopted reduced plasma model is the so-called hybrid-kinetic model, where the ions are fully kinetic and the electrons are treated as a neutralizing (inertial or massless) fluid.…
▽ More
Kinetic-range turbulence in magnetized plasmas and, in particular, in the context of solar-wind turbulence has been extensively investigated over the past decades via numerical simulations. Among others, one of the widely adopted reduced plasma model is the so-called hybrid-kinetic model, where the ions are fully kinetic and the electrons are treated as a neutralizing (inertial or massless) fluid. Within the same model, different numerical methods and/or approaches to turbulence development have been employed. In the present work, we present a comparison between two-dimensional hybrid-kinetic simulations of plasma turbulence obtained with two complementary approaches spanning about two decades in wavenumber - from MHD inertial range to scales well below the ion gyroradius - with a state-of-the-art accuracy. One approach employs hybrid particle-in-cell (HPIC) simulations of freely-decaying Alfvénic turbulence, whereas the other consists of Eulerian hybrid Vlasov-Maxwell (HVM) simulations of turbulence continuously driven with partially-compressible large-scale fluctuations. Despite the completely different initialization and injection/drive at large scales, the same properties of turbulent fluctuations at $k_\perpρ_i\gtrsim1$ are observed. The system indeed self-consistently "reprocesses" the turbulent fluctuations while they are cascading towards smaller and smaller scales, in a way which actually depends on the plasma beta parameter. Small-scale turbulence has been found to be mainly populated by kinetic Alfvén wave (KAW) fluctuations for $β\geq1$, whereas KAW fluctuations are only sub-dominant for low-$β$.
△ Less
Submitted 7 March, 2017;
originally announced March 2017.
-
Plasma beta dependence of the ion-scale spectral break of solar wind turbulence: high-resolution 2D hybrid simulations
Authors:
Luca Franci,
Simone Landi,
Lorenzo Matteini,
Andrea Verdini,
Petr Hellinger
Abstract:
We investigate properties of the ion-scale spectral break of solar wind turbulence by means of two-dimensional high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box and add a spectrum of in-plane, large-scale, magnetic and kinetic fluctuations. We perform a set of simulations with different values of the plasma beta, di…
▽ More
We investigate properties of the ion-scale spectral break of solar wind turbulence by means of two-dimensional high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box and add a spectrum of in-plane, large-scale, magnetic and kinetic fluctuations. We perform a set of simulations with different values of the plasma beta, distributed over three orders of magnitude, from 0.01 to 10. In all the cases, once turbulence is fully developed, we observe a power-law spectrum of the fluctuating magnetic field on large scales (in the inertial range) with a spectral index close to -5/3, while in the sub-ion range we observe another power-law spectrum with a spectral index systematically varying with $β$ (from around -3.6 for small values to around -2.9 for large ones). The two ranges are separated by a spectral break around ion scales. The length scale at which this transition occurs is found to be proportional to the ion inertial length, $d_i$, for $β\ll 1$ and to the ion gyroradius, $ρ_i = d_i\sqrtβ$, for $β\gg 1$, i.e., to the larger between the two scales in both the extreme regimes. For intermediate cases, i.e., $β\sim 1$, a combination of the two scales is involved. We infer an empiric relation for the dependency of the spectral break on $β$ that provides a good fit over the whole range of values. We compare our results with in situ observations in the solar wind and suggest possible explanations for such a behavior.
△ Less
Submitted 17 October, 2016;
originally announced October 2016.
-
Two-dimensional Hybrid Simulations of Kinetic Plasma Turbulence: Current and Vorticity vs Proton Temperature
Authors:
Luca Franci,
Petr Hellinger,
Lorenzo Matteini,
Andrea Verdini,
Simone Landi
Abstract:
Proton temperature anisotropies between the directions parallel and perpendicular to the mean magnetic field are usually observed in the solar wind plasma. Here, we employ a high-resolution hybrid particle-in-cell simulation in order to investigate the relation between spatial properties of the proton temperature and the peaks in the current density and in the flow vorticity. Our results indicate…
▽ More
Proton temperature anisotropies between the directions parallel and perpendicular to the mean magnetic field are usually observed in the solar wind plasma. Here, we employ a high-resolution hybrid particle-in-cell simulation in order to investigate the relation between spatial properties of the proton temperature and the peaks in the current density and in the flow vorticity. Our results indicate that, although regions where the proton temperature is enhanced and temperature anisotropies are larger correspond approximately to regions where many thin current sheets form, no firm quantitative evidence supports the idea of a direct causality between the two phenomena. On the other hand, quite a clear correlation between the behavior of the proton temperature and the out-of-plane vorticity is obtained.
△ Less
Submitted 4 April, 2016;
originally announced April 2016.
-
Plasma turbulence and kinetic instabilities at ion scales in the expanding solar wind
Authors:
Petr Hellinger,
Lorenzo Matteini,
Simone Landi,
Andrea Verdini,
Luca Franci,
Pavel M. Travnicek
Abstract:
The relationship between a decaying strong turbulence and kinetic instabilities in a slowly expanding plasma is investigated using two-dimensional (2-D) hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvénic fluctuations which have energy equipartition…
▽ More
The relationship between a decaying strong turbulence and kinetic instabilities in a slowly expanding plasma is investigated using two-dimensional (2-D) hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvénic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at ion scales. The turbulent cascade leads to an overall anisotropic proton heating, protons are heated in the perpendicular direction, and, initially, also in the parallel direction. The imposed expansion leads to generation of a large parallel proton temperature anisotropy which is at later stages partly reduced by turbulence. The turbulent heating is not sufficient to overcome the expansion-driven perpendicular cooling and the system eventually drives the oblique firehose instability in a form of localized nonlinear wave packets which efficiently reduce the parallel temperature anisotropy. This work demonstrates that kinetic instabilities may coexist with strong plasma turbulence even in a constrained 2-D regime.
△ Less
Submitted 29 September, 2015; v1 submitted 13 August, 2015;
originally announced August 2015.
-
High-resolution hybrid simulations of kinetic plasma turbulence at proton scales
Authors:
Luca Franci,
Simone Landi,
Lorenzo Matteini,
Andrea Verdini,
Petr Hellinger
Abstract:
We investigate properties of plasma turbulence from magneto-hydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field, perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations, with energy equipartition and vanishing correlation. Once the turbu…
▽ More
We investigate properties of plasma turbulence from magneto-hydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field, perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations, with energy equipartition and vanishing correlation. Once the turbulence is fully developed, we observe a MHD inertial range, where the spectra of the perpendicular magnetic field and the perpendicular proton bulk velocity fluctuations exhibit power-law scaling with spectral indices of -5/3 and -3/2, respectively. This behavior is extended over a full decade in wavevectors and is very stable in time. A transition is observed around proton scales. At sub-ion scales, both spectra steepen, with the former still following a power law with a spectral index of ~-3. A -2.8 slope is observed in the density and parallel magnetic fluctuations, highlighting the presence of compressive effects at kinetic scales. The spectrum of the perpendicular electric fluctuations follows that of the proton bulk velocity at MHD scales, and flattens at small scales. All these features, which we carefully tested against variations of many parameters, are in good agreement with solar wind observations. The turbulent cascade leads to on overall proton energization with similar heating rates in the parallel and perpendicular directions. While the parallel proton heating is found to be independent on the resistivity, the number of particles per cell and the resolution employed, the perpendicular proton temperature strongly depends on these parameters.
△ Less
Submitted 19 June, 2015;
originally announced June 2015.
-
Solar wind turbulence from MHD to sub-ion scales: high-resolution hybrid simulations
Authors:
Luca Franci,
Andrea Verdini,
Lorenzo Matteini,
Simone Landi,
Petr Hellinger
Abstract:
We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several propert…
▽ More
We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several properties of the observed solar wind fluctuations: spectral indices of the magnetic, kinetic, and residual energy spectra in the magneto-hydrodynamic (MHD) inertial range along with a flattening of the electric field spectrum, an increase in magnetic compressibility, and a strong coupling of the cascade with the density and the parallel component of the magnetic fluctuations at sub-proton scales. Our findings support the interpretation that in the solar wind large-scale MHD fluctuations naturally evolve beyond proton scales into a turbulent regime that is governed by the generalized Ohm's law.
△ Less
Submitted 18 April, 2015; v1 submitted 18 March, 2015;
originally announced March 2015.
-
Neutron Star instabilities in full General Relativity using a $Γ=2.75$ ideal fluid
Authors:
Roberto De Pietri,
Alessandra Feo,
Luca Franci,
Frank Löffler
Abstract:
We present results about the effect of the use of a stiffer equation of state, namely the ideal-fluid $Γ=2.75$ ones, on the dynamical bar-mode instability in rapidly rotating polytropic models of neutron stars in full General Relativity. We determine the change on the critical value of the instability parameter $β$ for the emergence of the instability when the adiabatic index $Γ$ is changed from 2…
▽ More
We present results about the effect of the use of a stiffer equation of state, namely the ideal-fluid $Γ=2.75$ ones, on the dynamical bar-mode instability in rapidly rotating polytropic models of neutron stars in full General Relativity. We determine the change on the critical value of the instability parameter $β$ for the emergence of the instability when the adiabatic index $Γ$ is changed from 2 to 2.75 in order to mimic the behavior of a realistic equation of state. In particular, we show that the threshold for the onset of the bar-mode instability is reduced by this change in the stiffness and give a precise quantification of the change in value of the critical parameter $β_c$. We also extend the analysis to lower values of $β$ and show that low-beta shear instabilities are present also in the case of matter described by a simple polytropic equation of state.
△ Less
Submitted 15 June, 2014; v1 submitted 31 March, 2014;
originally announced March 2014.
-
Bar-mode instability suppression in magnetized relativistic stars
Authors:
Luca Franci,
Roberto De Pietri,
Kyriaki Dionysopoulou,
Luciano Rezzolla
Abstract:
We show that magnetic fields stronger than about $10^{15}$ G are able to suppress the development of the hydrodynamical bar-mode instability in relativistic stars. The suppression is due to a change in the rest-mass density and angular velocity profiles due to the formation and to the linear growth of a toroidal component that rapidly overcomes the original poloidal one, leading to an amplificatio…
▽ More
We show that magnetic fields stronger than about $10^{15}$ G are able to suppress the development of the hydrodynamical bar-mode instability in relativistic stars. The suppression is due to a change in the rest-mass density and angular velocity profiles due to the formation and to the linear growth of a toroidal component that rapidly overcomes the original poloidal one, leading to an amplification of the total magnetic energy. The study is carried out performing three-dimensional ideal-magnetohydrodynamics simulations in full general relativity, superimposing to the initial (matter) equilibrium configurations a purely poloidal magnetic field in the range $10^{14}-10^{16}$ G. When the seed field is a few parts in $10^{15}$ G or above, all the evolved models show the formation of a low-density envelope surrounding the star. For much weaker fields, no effect on the matter evolution is observed, while magnetic fields which are just below the suppression threshold are observed to slow down the growth-rate of the instability.
△ Less
Submitted 25 September, 2013;
originally announced September 2013.
-
Dynamical bar-mode instability in rotating and magnetized relativistic stars
Authors:
Luca Franci,
Roberto De Pietri,
Kyriaki Dionysopoulou,
Luciano Rezzolla
Abstract:
We present three-dimensional simulations of the dynamical bar-mode instability in magnetized and differentially rotating stars in full general relativity. Our focus is on the effects that magnetic fields have on the dynamics and the onset of the instability. In particular, we perform ideal-magnetohydrodynamics simulations of neutron stars that are known to be either stable or unstable against the…
▽ More
We present three-dimensional simulations of the dynamical bar-mode instability in magnetized and differentially rotating stars in full general relativity. Our focus is on the effects that magnetic fields have on the dynamics and the onset of the instability. In particular, we perform ideal-magnetohydrodynamics simulations of neutron stars that are known to be either stable or unstable against the purely hydrodynamical instability, but to which a poloidal magnetic field in the range of $10^{14}$--$10^{16}$ G is superimposed initially. As expected, the differential rotation is responsible for the shearing of the poloidal field and the consequent linear growth in time of the toroidal magnetic field. The latter rapidly exceeds in strength the original poloidal one, leading to a magnetic-field amplification in the the stars. Weak initial magnetic fields, i.e. $ \lesssim 10^{15}$ G, have negligible effects on the development of the dynamical bar-mode instability, simply braking the stellar configuration via magnetic-field shearing, and over a timescale for which we derived a simple algebraic expression. On the other hand, strong magnetic fields, i.e. $\gtrsim 10^{16}$ G, can suppress the instability completely, with the precise threshold being dependent also on the amount of rotation. As a result, it is unlikely that very highly magnetized neutron stars can be considered as sources of gravitational waves via the dynamical bar-mode instability.
△ Less
Submitted 19 August, 2013;
originally announced August 2013.