Related papers: Field theory and structure-preserving geometric pa…
Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in…
Laboratory plasmas in open magnetic geometries can be found in many different applications such as (1) Scrape-Of-Layer (SOL) and divertor regions in toroidal confinement fusion devices (\approx1-10^2\hspace{1mm}\mathrm{eV}), (2) linear…
We propose a novel structure preserving discretization for viscous and resistive magnetohydrodynamics. We follow the recent line of work on discrete least action principle for fluid and plasma equation, incorporating the recent advances to…
The transport of heat and particles in the relatively collisional edge regions of magnetically confined plasmas is a scientifically challenging and technologically important problem. Understanding and predicting this transport requires the…
Variational symplectic algorithms have recently been developed for carrying out long-time simulation of charged particles in magnetic fields. As a direct consequence of their derivation from a discrete variational principle, these…
We consider the issue of strict, fully discrete \emph{local} energy conservation for a whole class of fully implicit local-charge- and global-energy-conserving particle-in-cell (PIC) algorithms. Earlier studies demonstrated these algorithms…
We present particle-in-cell (PIC) simulations of minority energetic protons in deuterium plasmas, which demonstrate a collective instability responsible for emission near the lower hybrid frequency and its harmonics. The simulations capture…
An axial-azimuthal two-dimensional Hall thruster discharge model was developed for analyzing gradient drift instability (GDI) and cross-field electron transport enhancement induced solely by the GDI. A hybrid particle-fluid model was used…
The dissipation mechanism of the magnetic energy in turbulent collisionless space and astrophysical plasmas is still not well understood. Its investigation requires efficient kinetic simulations of the energy transfer in collisionless…
Basic, local kinetic theory of ion temperature gradient driven (ITG) mode, with adiabatic electrons is reconsidered. Standard unstable, purely oscillating as well as damped solutions of the local dispersion relation are obtained using a…
Over the past decades, multiple gyrokinetic codes have shown to be able to simulate turbulence and associated transport in the core of Tokamak devices. However, their application to the edge and scrape-off layer (SOL) region presents…
Modeling of kinetic plasmas using electromagnetic particle in cell methods (EM-PIC) is a problem that is well worn, in that methods developed have been used extensively both understanding physics and exploiting them for device design.…
Shear flows, ubiquitous in space and astrophysical plasmas, can accelerate particles through turbulence excited by the Kelvin-Helmholtz instability. We present the first numerical study of particle acceleration in non-relativistic,…
Electromagnetic particle simulation model has been formulated and verified for nonlinear processes of lower hybrid (LH) waves in fusion plasmas. Electron dynamics is described by the drift kinetic equation using either kinetic momentum or…
Kinetic Particle In Cell (PIC) methods can extend greatly their range of applicability if implicit time differencing and spatial adaption are used to address the wide range of time and length scales typical of plasmas. For implicit…
The paper presents a 1D2V electrostatic PIC model with a drift-kinetic description of all particle types aiming at simulating classical longitudinal plasma transport in axially symmetric open traps. The model generalizes the semi-implicit…
A geometrical method is used for the analysis of stochastic processes in plasma turbulence. Distances between thermodynamic states can be computed according the thermodynamic length methodology which allows the use of a Riemannian metric on…
The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence…
The present paper is devoted to the convergence analysis of an asymptotic preserving particle scheme designed to serve as a particle pusher in a Particle-In-Cell (PIC) method for the Vlasov equation with a strong inhomogeneous magnetic…
We discuss structure-preserving numerical discretizations for repulsive and attractive Euler-Poisson equations that find applications in fluid-plasma and self-gravitation modeling. The scheme is fully discrete and structure preserving in…