Source author record

D. Told

D. Told 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

11works
3topics
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

11 published item(s)

preprint2026arXiv

Plasmoid formation via competing lower-hybrid drift and Kelvin-Helmholtz instabilities: A hybrid kinetic-gyrokinetic simulation study

We investigate the nonlinear formation of plasmoids in 2D low-beta current sheets through the interplay between the Kelvin-Helmholtz instability (KHI) and the lower-hybrid drift instability (LHDI). Using a hybrid kinetic-gyrokinetic model-based Super Simple Vlasov (ssV) code with fully kinetic ions and drift-kinetic electrons, we simulate Harris-type current sheets and velocity shear layers with strong cross-field density gradients. Our central hypothesis is that steep density gradients drive LHDI, which can grow faster than KHI and initiate an inverse cascade from kinetic to fluid scales, potentially suppressing KHI. Our simulations confirm that, in thin current sheets, LHDI develops rapidly at the sheet edges and nonlinearly merges into larger-scale magnetic islands before KHI can evolve. These LHDI-driven structures distort the velocity shear and suppress classical KH vortices. In contrast, for thicker current sheets or weaker density gradients, KHI dominates and produces the expected rolled-up vortices and associated plasmoids. These findings demonstrate that LHDI-induced turbulence can act as both a seed and a regulator of plasmoid-generating instabilities, mediating cross-scale energy transfer. This mechanism is relevant to thin boundary layers in space plasmas, such as the solar wind magnetosphere interface, and suggests that microturbulence can govern large-scale magnetic topology during collisionless reconnection.

preprint2025arXiv

Verification of a hybrid gyrokinetic model using the advanced semi-Lagrange code ssV

The super simple Vlasov (ssV) code was developed to study instabilities, turbulence, and reconnection in weakly magnetized plasmas, such as the solar wind in the dissipation range and the edge of fusion plasmas. The ssV code overcomes the limitations of standard gyrokinetic theory by using a hybrid model that incorporates fully kinetic ions and gyrokinetic electrons. This hybrid gyrokinetic model enables accurate modeling in regimes characterized by steep gradients and high-frequency dynamics. To achieve this, ssV implements a set of semi-Lagrangian numerical schemes, including Positive Flux Conservative (PFC), Flux Conservative fifth-order (FCV), FCV with Umeda limiters, and a Semi-Lagrangian Monotonicity-Preserving fifth-order scheme (SLMP5). Benchmark problems such as Landau damping, ion-acoustic waves, ion Bernstein waves, and kinetic Alfven waves were employed to evaluate the schemes. The SLMP5 scheme consistently delivered the best overall accuracy and numerical stability performance. The code also addresses well-known electromagnetic gyrokinetic simulation issues, such as the Ampere cancellation problem, using carefully chosen velocity-space resolutions and accurate integral evaluation.

preprint2016arXiv

Subproton-scale cascades in solar wind turbulence: driven hybrid-kinetic simulations

A long-lasting debate in space plasma physics concerns the nature of subproton-scale fluctuations in solar wind (SW) turbulence. Over the past decade, a series of theoretical and observational studies were presented in favor of either kinetic Alfvén wave (KAW) or whistler turbulence. Here, we investigate numerically the nature of the subproton-scale turbulent cascade for typical SW parameters by means of unprecedented high-resolution simulations of forced hybrid-kinetic turbulence in two real-space and three velocity-space dimensions. Our analysis suggests that small-scale turbulence in this model is dominated by KAWs at $β\gtrsim1$ and by magnetosonic/whistler fluctuations at lower $β$. The spectral properties of the turbulence appear to be in good agreement with theoretical predictions. A tentative interpretation of this result in terms of relative changes in the damping rates of the different waves is also presented. Overall, the results raise interesting new questions about the properties and variability of subproton-scale turbulence in the SW, including its possible dependence on the plasma $β$, and call for detailed and extensive parametric explorations of driven kinetic turbulence in three dimensions.

preprint2016arXiv

The Structure of Plasma Heating in Gyrokinetic Alfvénic Turbulence

We analyze plasma heating in weakly collisional kinetic Alfvén wave (KAW) turbulence using high resolution gyrokinetic simulations spanning the range of scales between the ion and the electron gyroradii. Real space structures that have a higher than average heating rate are shown not to be confined to current sheets. This novel result is at odds with previous studies, which use the electromagnetic work in the local electron fluid frame, i.e. $\mathbf{J} \!\cdot\! (\mathbf{E} + \mathbf{v}_e\times\mathbf{B})$, as a proxy for turbulent dissipation to argue that heating follows the intermittent spatial structure of the electric current. Furthermore, we show that electrons are dominated by parallel heating while the ions prefer the perpendicular heating route. We comment on the implications of the results presented here.

preprint2015arXiv

Gyrokinetic and kinetic particle-in-cell simulations of guide-field reconnection. I: Macroscopic effects of the electron flows

In this work, we compare gyrokinetic (GK) and fully kinetic Particle-in-Cell (PIC) simulations of magnetic reconnection in the limit of strong guide field. In particular, we analyze the limits of applicability of the GK plasma model compared to a fully kinetic description of force free current sheets for finite guide fields ($b_g$). Here we report the first part of an extended comparison, focusing on the macroscopic effects of the electron flows. For a low beta plasma ($β_i=0.01$), it is shown that both plasma models develop magnetic reconnection with similar features in the secondary magnetic islands if a sufficiently high guide field ($b_g\gtrsim 30$) is imposed in the kinetic PIC simulations. Outside of these regions, in the separatrices close to the X points, the convergence between both plasma descriptions is less restrictive ($b_g\gtrsim 5$). Kinetic PIC simulations using guide fields $b_g \lesssim 30$ reveal secondary magnetic islands with a core magnetic field and less energetic flows inside of them in comparison to the GK or kinetic PIC runs with stronger guide fields. We find that these processes are mostly due to an initial shear flow absent in the GK initialization and negligible in the kinetic PIC high guide field regime, in addition to fast outflows on the order of the ion thermal speed that violate the GK ordering. Since secondary magnetic islands appear after the reconnection peak time, a kinetic PIC/GK comparison is more accurate in the linear phase of magnetic reconnection. For a high beta plasma ($β_i=1.0$) where reconnection rates and fluctuations levels are reduced, similar processes happen in the secondary magnetic islands in the fully kinetic description, but requiring much lower guide fields ($b_g\lesssim 3$).

preprint2015arXiv

Gyrokinetic studies of core turbulence features in ASDEX Upgrade H-mode plasmas

Gyrokinetic validation studies are crucial in developing confidence in the model incorporated in numerical simulations and thus improving their predictive capabilities. As one step in this direction, we simulate an ASDEX Upgrade discharge with the GENE code, and analyze various fluctuating quantities and compare them to experimental measurements. The approach taken is the following. First, linear simulations are performed in order to determine the turbulence regime. Second, the heat fluxes in nonlinear simulations are matched to experimental fluxes by varying the logarithmic ion temperature gradient within the expected experimental error bars. Finally, the dependence of various quantities with respect to the ion temperature gradient is analyzed in detail. It is found that density and temperature fluctuations can vary significantly with small changes in this parameter, thus making comparisons with experiments very sensitive to uncertainties in the experimental profiles. However, cross-phases are more robust, indicating that they are better observables for comparisons between gyrokinetic simulations and experimental measurements.

preprint2015arXiv

Multiscale nature of the dissipation range in gyrokinetic simulations of Alfvénic turbulence

Nonlinear energy transfer and dissipation in Alfvén wave turbulence are analyzed in the first gyrokinetic simulation spanning all scales from the tail of the MHD range to the electron gyroradius scale. For typical solar wind parameters at 1 AU, about 30% of the nonlinear energy transfer close to the electron gyroradius scale is mediated by modes in the tail of the MHD cascade. Collisional dissipation occurs across the entire kinetic range $k_\perpρ_i\gtrsim 1$. Both mechanisms thus act on multiple coupled scales, which have to be retained for a comprehensive picture of the dissipation range in Alfvénic turbulence.

preprint2014arXiv

Collision-dependent power law scalings in 2D gyrokinetic turbulence

Nonlinear gyrokinetics provides a suitable framework to describe short-wavelength turbulence in magnetized laboratory and astrophysical plasmas. In the electrostatic limit, this system is known to exhibit a free energy cascade towards small scales in (perpendicular) real and/or velocity space. The dissipation of free energy is always due to collisions (no matter how weak the collisionality), but may be spread out across a wide range of scales. Here, we focus on freely-decaying 2D electrostatic turbulence on sub-ion-gyroradius scales. An existing scaling theory for the turbulent cascade in the weakly collisional limit is generalized to the moderately collisional regime. In this context, non-universal power law scalings due to multiscale dissipation are predicted, and this prediction is confirmed by means of direct numerical simulations.

preprint2014arXiv

Electromagnetic stabilization of tokamak microturbulence in a high-$β$ regime

The impact of electromagnetic stabilization and flow shear stabilization on ITG turbulence is investigated. Analysis of a low-$β$ JET L-mode discharge illustrates the relation between ITG stabilization, and proximity to the electromagnetic instability threshold. This threshold is reduced by suprathermal pressure gradients, highlighting the effectiveness of fast ions in ITG stabilization. Extensive linear and nonlinear gyrokinetic simulations are then carried out for the high-$β$ JET hybrid discharge 75225, at two separate locations at inner and outer radii. It is found that at the inner radius, nonlinear electromagnetic stabilization is dominant, and is critical for achieving simulated heat fluxes in agreement with the experiment. The enhancement of this effect by suprathermal pressure also remains significant. It is also found that flow shear stabilization is not effective at the inner radii. However, at outer radii the situation is reversed. Electromagnetic stabilization is negligible while the flow shear stabilization is significant. These results constitute the high-$β$ generalization of comparable observations found at low-$β$ at JET. This is encouraging for the extrapolation of electromagnetic ITG stabilization to future devices. An estimation of the impact of this effect on the ITER hybrid scenario leads to a 20% fusion power improvement.

preprint2013arXiv

Ion temperature profile stiffness: non-linear gyrokinetic simulations and comparison with experiment

Recent experimental observations at JET show evidence of reduced ion temperature profile stiffness, hypothesised to be due to concomitant low magnetic shear (s) and significant toroidal rotational flow shear. Non-linear gyrokinetic simulations are performed, aiming to investigate the physical mechanism behind the observations. A comprehensive set of simulations are carried out, comparing the impact on the ion heat flux of various parameters that differ within the data-set. These parameters include q, s, rotation, effect of rotation on the magnetohydrodynamic (MHD) equilibrium, R/L_n, beta_e, Z_eff, and the fast particle content. The effect of toroidal flow shear itself is not predicted by the simulations to lead to a significant reduction in ion heat flux, due both to an insufficient magnitude of flow shear and significant parallel velocity gradient destabilisation. It is however found that non-linear electromagnetic effects due to both thermal and fast-particle pressure gradients, even at low beta_e, can significantly reduce the profile stiffness. A total of five discharges are examined, at both inner and outer radii. For all cases studied, the simulated and experimental ion heat flux values agree within reasonable variations of input parameters around the experimental uncertainties.

preprint2013arXiv

Nonlinear stabilization of tokamak microturbulence by fast ions

Nonlinear electromagnetic stabilization by suprathermal pressure gradients found in specific regimes is shown to be a key factor in reducing tokamak microturbulence, augmenting significantly the thermal pressure electromagnetic stabilization. Based on nonlinear gyrokinetic simulations investigating a set of ion heat transport experiments on the JET tokamak, described by Mantica et al. [Phys. Rev. Lett. 107 135004 (2011)], this result explains the experimentally observed ion heat flux and stiffness reduction. These findings are expected to improve the extrapolation of advanced tokamak scenarios to reactor relevant regimes.