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

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

References in corpus (6)

• Variational integration for ideal magnetohydrodynamics with built-in advection equations
• Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems
• Variational approach to low-frequency kinetic-MHD in the current coupling scheme
• Toroidal regularization of the guiding center Lagrangian
• Gauge-free electromagnetic gyrokinetic theory
• Subcycling of particle orbits in variational, geometric electromagnetic particle-in-cell methods

Cited by in corpus (9)

• A quantum-inspired method for solving the Vlasov-Poisson equations
• Rubrics for Charge Conserving Current Mapping in Finite Element Particle in Cell Methods
• Time Integrator Agnostic Charge Conserving Finite Element PIC
• Subcycling of particle orbits in variational, geometric electromagnetic particle-in-cell methods
• Advanced fuel fusion, phase space engineering, and structure-preserving geometric algorithms
• Geometric Electrostatic Particle-In-Cell Algorithm on Unstructured Meshes
• Multispecies structure-preserving particle discretization of the Landau collision operator
• Metriplectic foundations of gyrokinetic Vlasov-Maxwell-Landau theory
• A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus