Source author record

Emanuele Papini

Emanuele Papini 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

9works
7topics
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

9 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

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.

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

Simulating acoustic waves in spotted stars

Acoustic modes of oscillation are affected by stellar activity, however it is unclear how starspots contribute to these changes. Here we investigate the nonmagnetic effects of starspots on global modes with angular degree $\ell \leq 2$ in highly active stars, and characterize the spot seismic signature on synthetic light curves. We perform 3D time-domain simulations of linear acoustic waves to study their interaction with a model starspot. We model the spot as a 3D change in the sound speed stratification with respect to a convectively stable stellar background, built from solar Model S. We perform a parametric study by considering different depths and perturbation amplitudes. Exact numerical simulations allow the investigation of the wavefield-spot interaction beyond first order perturbation theory. The interaction of the axisymmetric modes with the starspot is strongly nonlinear. As mode frequency increases, the frequency shifts for radial modes exceed the value predicted by linear theory, while the shifts for the $\ell=2, m=0$ modes are smaller than predicted by linear theory, with avoided-crossing-like patterns forming between the $m=0$ and $m=1$ mode frequencies. The nonlinear behavior increases with increasing spot amplitude and/or decreasing depth. Linear theory still reproduces the correct shifts for nonaxisymmetric modes. In the nonlinear regime the mode eigenfunctions are not pure spherical harmonics, but rather a mixture of different spherical harmonics. This mode mixing, together with the frequency changes, may lead to misidentification of the modes in the observed acoustic power spectra.

preprint2013arXiv

Propagating Linear Waves in Convectively Unstable Stellar Models: a Perturbative Approach

Linear time-domain simulations of acoustic oscillations are unstable in the stellar convection zone. To overcome this problem it is customary to compute the oscillations of a stabilized background stellar model. The stabilization, however, affects the result. Here we propose to use a perturbative approach (running the simulation twice) to approximately recover the acoustic wave field, while preserving seismic reciprocity. To test the method we considered a 1D standard solar model. We found that the mode frequencies of the (unstable) standard solar model are well approximated by the perturbative approach within $1~μ$Hz for low-degree modes with frequencies near $3~μ$Hz. We also show that the perturbative approach is appropriate for correcting rotational-frequency kernels. Finally, we comment that the method can be generalized to wave propagation in 3D magnetized stellar interiors because the magnetic fields have stabilizing effects on convection.