Source author record

Simone Landi

Simone Landi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2026arXiv

Self-organization of local streamline structures and energy transfer rate in compressible plasma turbulence

We examine how local streamline topology and energy cascade rate self-organize in plasma turbulence for both compressible and incompressible regimes. Using a fully-compressible Hall-magnetohydrodynamic simulation, we quantify the subgrid-scale energy transfer and analyze its relationship to streamline structures by means of grandient tensor geometric invariants of the velocity field. Our results highlight how streamline topology is crucial for diagnosing turbulence: for nearly-incompressible fluctuations the energy is primarily transferred to smaller scales through strain-dominated and stable-vortical structures, while is back-transferred towards larger scales through unstable-vortical structures. Compressible fluctuations, on the contrary, do not show a clear topological selection of the energy transfer since the overall direction of the local cascade rate is found to be determined by the sign of $-\nabla\cdot u$ (plasma volumetric compression or expansion).

preprint2024arXiv

A GPU-Accelerated Modern Fortran Version of the ECHO Code for Relativistic Magnetohydrodynamics

The numerical study of relativistic magnetohydrodynamics (MHD) plays a crucial role in high-energy astrophysics, but unfortunately is computationally demanding, given the complex physics involved (high Lorentz factor flows, extreme magnetization, curved spacetimes near compact objects) and the large variety of spatial scales needed to resolve turbulent motions. A great benefit comes from the porting of existing codes running on standard processors to GPU-based platforms. However, this usually requires a drastic rewriting of the original code, the use of specific languages like CUDA, and a complex analysis of data management and optimization of parallel processes. Here we describe the porting of the ECHO code for special and general relativistic MHD to accelerated devices, simply based on native Fortran language built-in constructs, especially 'do concurrent' loops, few OpenACC directives, and the straightforward data management provided by the Unified Memory option of NVIDIA compilers.Thanks to these very minor modifications to the original code, the new version of ECHO runs at least 16 times faster on GPU platforms compared to CPU-based ones. The chosen benchmark is the 3D propagation of a relativistic MHD Alfvén wave, for which strong and weak scaling tests performed on the LEONARDO pre-exascale supercomputer at CINECA are provided (using up to 256 nodes corresponding to 1024 GPUs, and over 14 billion cells). Finally, an example of high-resolution relativistic MHD Alfvénic turbulence simulation is shown, demonstrating the potential for astrophysical plasmas of the new GPU-based version of ECHO.

preprint2022arXiv

Anisotropic electron heating in turbulence-driven magnetic reconnection in the near-Sun solar wind

We perform a high-resolution two-dimensional fully-kinetic numerical simulation of a turbulent plasma system with observation-driven conditions, in order to investigate the interplay between turbulence, magnetic reconnection, and particle heating from ion to sub-electron scales in the near-Sun solar wind. We find that the power spectra of the turbulent plasma and electromagnetic fluctuations show multiple power-law intervals down to scales smaller than the electron gyroradius. Magnetic reconnection is observed to occur in correspondence of current sheets with a thickness of the order of the electron inertial length, which form and shrink due to interacting ion-scale vortexes. In some cases, both ion and electron outflows are observed (the classic reconnection scenario), while in others -- typically for the shortest current sheets -- only electron jets are presents ("electron-only reconnection"). At the onset of reconnection, the electron temperature starts to increase and a strong parallel temperature anisotropy develops. This suggests that in strong turbulence electron-scale coherent structures may play a significant role for electron heating, as impulsive and localized phenomena such as magnetic reconnection may transfer energy from the electromagnetic fields to particles more efficiently than damping mechanisms related to interactions with wave-like fluctuations.

preprint2022arXiv

Ion-scale transition of plasma turbulence: Pressure-strain effect

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.

preprint2020arXiv

Coronal Electron Temperature inferred from the Strahl Electrons in the Inner Heliosphere: Parker Solar Probe and Helios observations

The shape of the electron velocity distribution function plays an important role in the dynamics of the solar wind acceleration. Electrons are normally modelled with three components, the core, the halo, and the strahl. We investigate how well the fast strahl electrons in the inner heliosphere preserve the information about the coronal electron temperature at their origin. We analysed the data obtained by two missions, Helios spanning the distances between 65 and 215 R$_S$, and Parker Solar Probe (PSP) reaching down to 35 R$_S$ during its first two orbits around the Sun. The electron strahl was characterised with two parameters, pitch-angle width (PAW), and the strahl parallel temperature (T$_{s\parallel}$). PSP observations confirm the already reported dependence of strahl PAW on core parallel plasma beta ($β_{ec\parallel}$)\citep{Bercic2019}. Most of the strahl measured by PSP appear narrow with PAW reaching down to 30$^o$. The portion of the strahl velocity distribution function aligned with the magnetic field is for the measured energy range well described by a Maxwellian distribution function. T$_{s\parallel}$ was found to be anti-correlated with the solar wind velocity, and independent of radial distance. These observations imply that T$_{s\parallel}$ carries the information about the coronal electron temperature. The obtained values are in agreement with coronal temperatures measured using spectroscopy (David et al. 2998), and the inferred solar wind source regions during the first orbit of PSP agree with the predictions using a PFSS model (Bale et al. 2019, Badman et al. 2019).

preprint2020arXiv

In-situ observation of Hall Magnetohydrodynamic Cascade in Space Plasma

We present estimates of the turbulent energy cascade rate, derived from a Hall-MHD third-order law. We compute the contribution from the Hall term and the MHD term to the energy flux. We use MMS data accumulated in the magnetosheath and the solar wind, and compare the results with previously established simulation results. We find that in observation, the MHD contribution is dominant at inertial scales, as in the simulations, but the Hall term becomes significant in observations at larger scales than in the simulations. Possible reasons are offered for this unanticipated result.

preprint2020arXiv

On cascade of kinetic energy in compressible hydrodynamic turbulence

Properties of the turbulent cascade of kinetic energy are studied using direct numerical simulations of three-dimensional hydrodynamic decaying turbulence with a moderate Reynolds number and the initial Mach number $M=1$. Compressible and incompressible versions of the Karman-Howarth-Monin (KHM) and low-pass filtering/coarse-graining approaches are compared. In the simulation the total energy is well conserved; the scale dependent KHM and coarse-grained energy equations are also well conserved; the two approaches show similar results, the system does not have an inertial range for the cascade of kinetic energy, the region where this cascade dominates also have a non-negligible contribution of the kinetic-energy decay, dissipation, and pressure-dilatation effects. While the two approaches give semi-quantitatively similar results for the kinetic energy cascade, dissipation and pressure-dilatation rates, they differ in the increment separation and filtering scales; these scales are not simply related. The two approaches may be used to find the inertial range and to determine the cascade/dissipation rate of the kinetic energy.

preprint2019arXiv

Modeling Kelvin-Helmholtz instability-driven turbulence with hybrid simulations of Alfvénic turbulence

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.

preprint2016arXiv

"Ideally" unstable current sheets and the triggering of fast magnetic reconnection

Magnetic reconnection is thought to be the dynamical mechanism underlying many explosive phenomena observed both in space and in the laboratory, though the question of how fast magnetic reconnection is triggered in such high Lundquist ($S$) number plasmas has remained elusive. It has been well established that reconnection can develop over timescales faster than those predicted traditionally once kinetic scales are reached. It has also been shown that, within the framework of resistive Magnetohydrodynamics (MHD), fast reconnection is achieved for thin enough sheets via the onset of the so-called plasmoid instability. The latter was discovered in studies specifically devoted to the Sweet-Parker current sheet, either as an initial condition or an apparent transient state developing in nonlinear studies. On the other hand, a fast tearing instability can grow on an ideal, i.e., $S$-independent, timescale (dubbed "ideal" tearing) within current sheets whose aspect ratio scales with the macroscopic Lundquist number as $L/a\sim S^{1/3}$ -- much smaller than the Sweet-Parker one -- suggesting a new way to approach to the initiation of fast reconnection in collapsing current configurations. Here we present an overview of what we have called "ideal" tearing in resistive MHD, and discuss how the same reasoning can be extended to other plasma models commonly used that include electron inertia and kinetic effects. We then discuss a scenario for the onset of "ideal" fast reconnection via collapsing current sheets and describe a quantitative model for the interpretation of the nonlinear evolution of "ideally" unstable sheets in two dimensions.

preprint2016arXiv

Plasma beta dependence of the ion-scale spectral break of solar wind turbulence: high-resolution 2D hybrid simulations

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.

preprint2016arXiv

Two-dimensional Hybrid Simulations of Kinetic Plasma Turbulence: Current and Vorticity vs Proton Temperature

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.

preprint2015arXiv

Anisotropy of third-order structure functions in MHD turbulence

The measure of the third-order structure function, Y, is employed in the solar wind to compute the cascade rate of turbulence. In the absence of a mean field B0=0, Y is expected to be isotropic (radial) and independent of the direction of increments, so its measure yields directly the cascade rate. For turbulence with mean field, as in the solar wind, Y is expected to become more two dimensional (2D), that is, to have larger perpendicular components, loosing the above simple symmetry. To get the cascade rate one should compute the flux of Y, which is not feasible with single-spacecraft data, thus measurements rely upon assumptions about the unknown symmetry. We use direct numerical simulations (DNS) of magneto-hydrodynamic (MHD) turbulence to characterize the anisotropy of Y. We find that for strong guide field B0=5 the degree of two-dimensionalization depends on the relative importance of shear and pseudo polarizations (the two components of an Alfvén mode in incompressible MHD). The anisotropy also shows up in the inertial range. The more Y is 2D, the more the inertial range extent differs along parallel and perpendicular directions. We finally test the two methods employed in observations and find that the so-obtained cascade rate may depend on the angle between B0 and the direction of increments. Both methods yield a vanishing cascade rate along the parallel direction, contrary to observations, suggesting a weaker anisotropy of solar wind turbulence compared to our DNS. This could be due to a weaker mean field and/or to solar wind expansion.

preprint2015arXiv

Fire Hose instability driven by alpha particle temperature anisotropy

We investigate properties of a solar wind-like plasma including a secondary alpha particle population exhibiting a parallel temperature anisotropy with respect to the background magnetic field, using linear and quasi-linear predictions and by means of one-dimensional hybrid simulations. We show that anisotropic alpha particles can drive a parallel fire hose instability analogous to that generated by protons, but that, remarkably, the instability can be triggered also when the parallel plasma beta of alpha particles is below unity. The wave activity generated by the alpha anisotropy affects the evolution of the more abundant protons, leading to their anisotropic heating. When both ion species have sufficient parallel anisotropies both of them can drive the instability, and we observe generation of two distinct peaks in the spectra of the fluctuations, with longer wavelengths associated to alphas and shorter ones to protons. If a non-zero relative drift is present, the unstable modes propagate preferentially in the direction of the drift associated with the unstable species. The generated waves scatter particles and reduce their temperature anisotropy to marginally stable state, and, moreover, they significantly reduce the relative drift between the two ion populations. The coexistence of modes excited by both species leads to saturation of the plasma in distinct regions of the beta/anisotropy parameter space for protons and alpha particles, in good agreement with in situ solar wind observations. Our results confirm that fire hose instabilities are likely at work in the solar wind and limit the anisotropy of different ion species in the plasma.

preprint2015arXiv

High-resolution hybrid simulations of kinetic plasma turbulence at proton scales

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.

preprint2015arXiv

Plasma turbulence and kinetic instabilities at ion scales in the expanding solar wind

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.

preprint2015arXiv

Resistive magnetohydrodynamic simulations of the ideal tearing mode

We study the linear and nonlinear evolution of the tearing instability on thin current sheets by means of two-dimensional numerical simulations, within the framework of compressible, resistive magnetohydrodynamics. In particular we analyze the behavior of current sheets whose inverse aspect ratio scales with the Lundquist number $S$ as $S^{-1/3}$. This scaling has been recently recognized to yield the threshold separating fast, ideal reconnection, with an evolution and growth which are independent of $S$ provided this is high enough, as it should be natural having the ideal case as a limit for $S\to\infty$. Our simulations confirm that the tearing instability growth rate can be as fast as $γ\approx 0.6\,{τ_A}^{-1}$, where $τ_A$ is the ideal Alfvénic time set by the macroscopic scales, for our least diffusive case with $S=10^7$. The expected instability dispersion relation and eigenmodes are also retrieved in the linear regime, for the values of $S$ explored here. Moreover, in the nonlinear stage of the simulations we observe secondary events obeying the same critical scaling with $S$, here calculated on the \emph{local}, much smaller lengths, leading to increasingly faster reconnection. These findings strongly support the idea that in a fully dynamic regime, as soon as current sheets develop, thin and reach this critical threshold in their aspect ratio, the tearing mode is able to trigger plasmoid formation and reconnection on the local (ideal) Alfvénic timescales, as required to explain the explosive flaring activity often observed in solar and astrophysical plasmas.

preprint2015arXiv

Solar wind turbulence from MHD to sub-ion scales: high-resolution hybrid simulations

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.

preprint2012arXiv

On the unconstrained expansion of a spherical plasma cloud turning collisionless : case of a cloud generated by a nanometer dust grain impact on an uncharged target in space

Nano and micro meter sized dust particles travelling through the heliosphere at several hundreds of km/s have been repeatedly detected by interplanetary spacecraft. When such fast moving dust particles hit a solid target in space, an expanding plasma cloud is formed through the vaporisation and ionisation of the dust particles itself and part of the target material at and near the impact point. Immediately after the impact the small and dense cloud is dominated by collisions and the expansion can be described by fluid equations. However, once the cloud has reached micro-m dimensions, the plasma may turn collisionless and a kinetic description is required to describe the subsequent expansion. In this paper we explore the late and possibly collisionless spherically symmetric unconstrained expansion of a single ionized ion-electron plasma using N-body simulations. Given the strong uncertainties concerning the early hydrodynamic expansion, we assume that at the time of the transition to the collisionless regime the cloud density and temperature are spatially uniform. We do also neglect the role of the ambient plasma. This is a reasonable assumption as long as the cloud density is substantially higher than the ambient plasma density. In the case of clouds generated by fast interplanetary dust grains hitting a solid target some 10^7 electrons and ions are liberated and the in vacuum approximation is acceptable up to meter order cloud dimensions. ...

preprint2009arXiv

Three-dimensional evolution of magnetic and velocity shear driven instabilities in a compressible magnetized jet

The problem of three-dimensional combined magnetic and velocity shear driven instabilities of a compressible magnetized jet modeled with a plane neutral/current double vortex sheet in the framework of the resistive magnetohydrodynamics is addressed. The resulting dynamics given by the stream+current sheet interaction is analyzed and the effects of a variable geometry of the basic fields are considered. Depending on the basic asymptotic magnetic field configuration, a selection rule of the linear instability modes can be obtained. Hence, the system follows a two-stage path developing either through a fully three-dimensional dynamics with a rapid evolution of kink modes leading to a final turbulent state, or rather through a driving two-dimensional instability pattern that develops on parallel planes on which a reconnection+coalescence process takes place.