Source author record

Fulvio Militello

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

ResearcherUnclaimed source record

physics.plasm-ph

Catalog footprint

What is connected

5works

1topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2020arXiv

Blob interaction in 2D scrape off layer simulations

Interaction of coherent structures known as blobs in the scrape-off layer of magnetic confinement fusion devices is investigated. Isolated and interacting seeded blobs as well as full plasma turbulence are studied with a two dimensional fluid code. The features of the blobs (size, amplitude, position) are determined with a blob tracking algorithm, which identifies them as coherent structures above a chosen density threshold and compared to a conventional center of mass approach. The agreement of these two methods is shown to be affected by the parameters of the blob tracking algorithm. The benchmarked approach is then extended to a population of interacting plasma blobs with statistically distributed amplitudes, sizes and initial positions for different levels of intermittency. As expected, for decreasing intermittency, we observe an increasing number of blobs deviating from size-velocity scaling laws of perfectly isolated blobs. This is found to be caused by the interaction of blobs with the electrostatic potential of one another, leading to higher average blob velocities. The degree of variation from the picture of perfectly isolated blobs is quantified as a function of the average waiting time of the seeded blobs.

preprint2019arXiv

Dependence on ion temperature of shallow-angle magnetic presheaths with adiabatic electrons

The magnetic presheath is a boundary layer occurring when magnetized plasma is in contact with a wall and the angle $α$ between the wall and the magnetic field $\vec{B}$ is oblique. Here, we consider the fusion-relevant case of a shallow-angle, $α\ll 1$, electron-repelling sheath, with the electron density given by a Boltzmann distribution, valid for $α/ \sqrt{τ+1} \gg \sqrt{m_{\text{e}}/m_{\text{i}}}$, where $m_{\text{e}}$ is the electron mass, $m_{\text{i}}$ is the ion mass, $τ= T_{\text{i}}/ZT_{\text{e}}$, $T_{\text{e}}$ is the electron temperature, $T_{\text{i}}$ is the ion temperature, and $Z$ is the ionic charge state. The thickness of the magnetic presheath is of the order of a few ion sound Larmor radii $ρ_{\text{s}} = \sqrt{m_{\text{i}} \left(ZT_{\text{e}} + T_{\text{i}} \right) } / ZeB$, where $e$ is the proton charge and $B = |\vec{B}|$ is the magnitude of the magnetic field. We study the dependence on $τ$ of the electrostatic potential and ion distribution function in the magnetic presheath by using a set of prescribed ion distribution functions at the magnetic presheath entrance, parameterized by $τ$. The kinetic model is shown to be asymptotically equivalent to Chodura's fluid model at small ion temperature, $τ\ll 1$, for $|\ln α| > 3|\ln τ| \gg 1$. In this limit, despite the fact that fluid equations give a reasonable approximation to the potential, ion gyro-orbits acquire a spatial extent that occupies a large portion of the magnetic presheath. At large ion temperature, $τ\gg 1$, relevant because $T_{\text{i}}$ is measured to be a few times larger than $T_{\text{e}}$ near divertor targets of fusion devices, ions reach the Debye sheath entrance (and subsequently the wall) at a shallow angle whose size is given by $\sqrtα$ or $1/\sqrtτ$, depending on which is largest.

preprint2018arXiv

Solution to a collisionless shallow-angle magnetic presheath with kinetic ions

Using a kinetic model for the ions and adiabatic electrons, we solve a steady state, electron-repelling magnetic presheath in which a uniform magnetic field makes a small angle $α\ll 1$ (in radians) with the wall. The presheath characteristic thickness is the typical ion gyroradius $ρ_{\text{i}}$. The Debye length $λ_{\text{D}}$ and the collisional mean free path of an ion $λ_{\text{mfp}}$ satisfy the ordering $λ_{\text{D}} \ll ρ_{\text{i}} \ll αλ_{\text{mfp}}$, so a quasineutral and collisionless model is used. We assume that the electrostatic potential is a function only of distance from the wall, and it varies over the scale $ρ_{\text{i}}$. Using the expansion in $α\ll 1$, we derive an analytical expression for the ion density that only depends on the ion distribution function at the entrance of the magnetic presheath and the electrostatic potential profile. Importantly, we have added the crucial contribution of the orbits in the region near the wall. By imposing the quasineutrality equation, we derive a condition that the ion distribution function must satisfy at the magnetic presheath entrance --- the kinetic equivalent of the Chodura condition. Using an ion distribution function at the entrance of the magnetic presheath that satisfies the kinetic Chodura condition, we find numerical solutions for the self-consistent electrostatic potential, ion density and flow across the magnetic presheath for several values of $α$. Our numerical results also include the distribution of ion velocities at the Debye sheath entrance. We find that at small values of $α$ there are substantially fewer ions travelling with a large normal component of the velocity into the wall.

preprint2015arXiv

Investigation of the Effect of Resistivity on Scrape Off Layer Filaments using Three Dimensional Simulations

The propagation of filaments in the Scrape Off Layer (SOL) of tokamaks largely determine the plasma profiles in the region. In a conduction limited SOL, parallel temperature gradients are expected, such that the resistance to parallel currents is greater at the target than further upstream. Since the perpendicular motion of an isolated filament is largely determined by balance of currents that flow through it, this may be expected to affect filament transport. 3D simulations have thus been used to study the influence of enhanced parallel resistivity on the dynamics of filaments. Filaments with the smallest perpendicular length scales, which were inertially limited at low resistivity (meaning that polarization rather than parallel currents determine their radial velocities), were unaffected by resistivity. For larger filaments, faster velocities were produced at higher resistivities, due to two mechanisms. Firstly parallel currents were reduced and polarization currents were enhanced, meaning that the inertial regime extended to larger filaments, and secondly a potential difference formed along the parallel direction so that higher potentials were produced in the region of the filament for the same amount of current to flow into the sheath. These results indicate that broader SOL profiles could be produced at higher resistivities.

preprint2014arXiv

3D Simulations of Plasma Filaments in the Scrape Off Layer: A Comparison with Models of Reduced Dimensionality

This paper presents simulations of isolated 3D filaments in a slab geometry obtained using a 3D reduced fluid code. First, systematic scans were performed to investigate how the dynamics of a filament are affected by its amplitude, perpendicular size and parallel extent. The perpendicular size of the filament was found to have a strong influence on its motions, as it determined the relative importance of parallel currents to polarisation and viscous currents, whilst drift-wave instabilities were observed if the initial amplitude of the blob was increased sufficiently. Next, the 3D simulations were compared to 2D simulations using different parallel closures; namely, the sheath dissipation closure, which neglects parallel gradients, and the vorticity advection closure, which neglects the influence of parallel currents. The vorticity advection closure was found to not replicate the 3D perpendicular dynamics and overestimated the initial radial acceleration of all the filaments studied. In contrast, a more satisfactory comparison with the sheath dissipation closure was obtained, even in the presence of significant parallel gradients, where the closure is no longer valid. Specifically it captured the contrasting dynamics of filaments with different perpen- dicular sizes that were observed in the 3D simulations which the vorticity advection closure failed to replicate. However, neither closure successfully replicated the Boltzmann spinning effects and associated poloidal drift of the blob that was observed in the 3D simulations. Although the sheath dissipation closure was concluded to be more successful in replicating the 3D dynamics, it is emphasised that the vorticity closure may still be relevant for situations where the parallel current is inhibited from closing through the sheath due to effects such as increased resistivity near the targets.