Abstract:Explicit structure-preserving geometric Particle-in-Cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm [1-12], achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km/s, agreeing with experimental observations [13, 14] and analytical calculations [15, 16]. Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It shows that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times ...

Discretization form of the continuity condition at the polar axis, with application to the gyrokinetic simulation in a magnetic fusion torus

Electromagnetic Formulation of Global Gyrokinetic Particle Simulation in Toroidal Geometry

Gyrokinetic Particle Simulation of Microturbulence for General Magnetic Geometry and Experimental Profiles

Gyro-average Method for Global Gyrokinetic Particle Simulation in Realistic Tokamak Geometry

Gyrokinetic Simulation of Dissipative Trapped Electron Mode in Tokamak Edge

Turbulent Transport of Toroidal Angular Momentum in Fusion Plasmas

Global Gyrokinetic Simulation of Neoclassical Ambipolar Electric Field and Its Effects on Microturbulence in W7-X Stellarator

TH/P2-32 Gyrokinetic Simulation of Tokamak Edge Plasmas

High-order Discretization of a Gyrokinetic Vlasov Model in Edge Plasma Geometry

Asymptotically stable Particle-in-Cell methods for the magnetized Vlasov--Poisson equations in orthogonal curvilinear coordinates

Theory of Gyrokinetic Velocity Moment and Its Application for Zonal Flows in a Tokamak Plasma

An accurate SUPG-stabilized continuous Galerkin discretization for anisotropic heat flux in magnetic confinement fusion

A Proper Discretization of Hydrodynamic Equations in Cylindrical Coordinates for Astrophysical Simulations

Explicit Structure-Preserving Geometric Particle-in-Cell Algorithm in Curvilinear Orthogonal Coordinate Systems and Its Applications to Whole-Device 6D Kinetic Simulations of Tokamak Physics

A gyrokinetic model for the plasma periphery of tokamak devices

Gyrokinetic simulation of ITG turbulence with toroidal geometry including the magnetic axis by using field-aligned coordinates

Gyrokinetic Electromagnetic Particle Simulations in Triangular Meshes with C1 Finite Elements

On the Gyrokinetic Equilibrium

Linear Gyrokinetic Theory for Kinetic Magnetohydrodynamic Eigenmodes in Tokamak Plasmas

Linear Gyrokinetic Theory of Two-Dimensional Ion-Temperature-gradient Mode in Tokamaks

Nonlinear Gyrokinetic Theory with Polarization Drift