Related papers: Second order Gyrokinetic theory for Particle-In-Ce…
Recent development of structure-preserving geometric particle-in-cell (PIC) algorithms for Vlasov-Maxwell systems is summarized. With the arriving of 100 petaflop and exaflop computing power, it is now possible to carry out direct…
Gyrokinetic simulations have greatly improved our theoretical understanding of turbulent transport in fusion devices. Most gyrokinetic models in use are delta-f simulations in which the slowly varying radial profiles of density and…
The Gyrokinetic Toroidal Code at Princeton (GTC-P) is a highly scalable and portable particle-in-cell (PIC) code. It solves the 5D Vlasov-Poisson equation featuring efficient utilization of modern parallel computer architectures at the…
Two-dimensional gyrokinetics is a simple paradigm for the study of kinetic magnetised plasma turbulence. In this paper, we present a comprehensive theoretical framework for this turbulence. We study both the inverse and direct cascades (the…
We reanalyse an arbitrary-wavelength gyrokinetic formalism [A. M. Dimits, Phys. Plasmas $\bf17$, 055901 (2010)], which orders only the vorticity to be small and allows strong, time-varying flows on medium and long wavelengths. We obtain a…
We present an energy-conserving numerical scheme to solve the Vlasov-Maxwell (VM) system based on the regularized moment method proposed in [Z. Cai, Y. Fan, and R. Li. CPAM, 2014]. The globally hyperbolic moment system is deduced for the…
For more than half a century, most of the plasma scientists have encountered a violation of the conservation laws of charge, momentum, and energy whenever they have numerically solve the first-principle equations of kinetic plasmas, such as…
This paper describes a new code for simulating astrophysical plasmas that solves a hybrid model composed of gyrokinetic ions (GKI) and an isothermal electron fluid (ITEF) [A. Schekochihin et al., Astrophys. J. Suppl. \textbf{182}, 310…
A basic prerequisite for the investigation of relativistic astrophysical magnetoplasmas, occurring typically in the vicinity of massive stellar objects (black holes, neutron stars, active galactic nuclei, etc.), is the accurate description…
We present the standard electromagnetic Particle-in-Cell method, starting from the discrete approximation of derivatives on a uniform grid. The application to second-order, centered, finite-difference discretization of the equations of…
The nonlinear gyrokinetic equations describe plasma turbulence in laboratory and astrophysical plasmas. To solve these equations, massively parallel codes have been developed and run on present-day supercomputers. This paper describes…
We describe Pegasus, a new hybrid-kinetic particle-in-cell code tailored for the study of astrophysical plasma dynamics. The code incorporates an energy-conserving particle integrator into a stable, second-order--accurate, three-stage…
The gyro-kinetic model is an approximation of the Vlasov-Maxwell system in a strongly magnetized magnetic field. We propose a new algorithm for solving it combining the Semi-Lagrangian (SL) method and the Arakawa (AKW) scheme with a…
The global total-$f$ gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas, implements a continuum grid to perform the dissipative operations, such as plasma collisions. To transfer the distribution…
Predicting the dynamics of a thermonuclear plasma during a magnetic confinement experiment is fundamental in order to make nuclear fusion a reliable source of energy. The development of a set of equations describing the plasma evolution on…
We explore the possibilities of applying structure-preserving numerical methods to a plasma hybrid model with kinetic ions and mass-less fluid electrons satisfying the quasi-neutrality relation. The numerical schemes are derived by finite…
The global and electromagnetic gyrokinetic particle-in-cell code ORB5 is employed to investigate the self-consistent interactions between Alfv\'en modes (AM) and ion temperature gradient (ITG) turbulence in a magnetically confined plasma.…
We introduce a simplified model of the electron-beam/plasma system to model the electrical breakdown caused by the inductive electric field created by a rapidly rising electron beam current. The rigid-beam model is a reduction to the…
Plasmas are highly nonlinear and multi-scale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD). Galerkin models,…
In this paper, we extend the geometric Particle in Cell framework on dual grids to a gauge-free drift-kinetic Vlasov--Maxwell model and its coupling with the fully kinetic model. We derive a discrete action principle on dual grids for our…