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 (39)

• An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm
• GEMPIC: Geometric ElectroMagnetic Particle-In-Cell Methods
• Structure and structure-preserving algorithms for plasma physics
• Variational formulation of particle algorithms for kinetic plasma simulations
• Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme
• A Hamiltonian splitting for the Vlasov-Maxwell system
• Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems
• Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov-Maxwell system
• Exact Charge-Conserving Scatter-Gather Algorithm for Particle-in-Cell Simulations on Unstructured Grids: A Geometric Perspective
• Structure-Preserving Geometric Particle-in-Cell Methods for Vlasov-Maxwell Systems
• Explicit high-order symplectic integrators for charged particles in general electromagnetic fields
• Global Gyrokinetic Simulation of Tokamak Edge Pedestal Instabilities
• Explicit symplectic algorithms based on generating functions for charged particle dynamics
• Hamiltonian integration methods for Vlasov-Maxwell equations
• Energy-conserving time propagation for a geometric particle-in-cell Vlasov--Maxwell solver
• Variational integration for ideal magnetohydrodynamics with built-in advection equations
• Hamiltonian particle-in-cell methods for Vlasov-Maxwell equations
• Hamiltonian formulation of the gyrokinetic Vlasov-Maxwell equations
• Finite-dimensional collisionless kinetic theory
• Variational approach to low-frequency kinetic-MHD in the current coupling scheme
• Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems
• Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm
• Toroidal regularization of the guiding center Lagrangian
• Formation of current singularity in a topologically constrained plasma
• Degenerate Variational Integrators for Magnetic Field Line Flow and Guiding Center Trajectories
• Simulations of relativistic-quantum plasmas using real-time lattice scalar QED
• Gauge-free electromagnetic gyrokinetic theory
• Canonical symplectic structure and structure-preserving geometric algorithms for Schrödinger-Maxwell systems
• Comment on "Symplectic integration of magnetic systems": a proof that the Boris algorithm is not variational
• The geometric theory of charge conservation in particle-in-cell simulations
• Effective action approach to wave propagation in scalar QED plasmas
• Comment on "Hamiltonian splitting for the Vlasov-Maxwell equations"
• Field theory and structure-preserving geometric particle-in-cell algorithm for drift wave instability and turbulence
• Explicit High-Order Gauge-Independent Symplectic Algorithms for Relativistic Charged Particle Dynamics
• A Spectral Canonical Electrostatic Algorithm
• Local Energy Conservation Law for Spatially-Discretized Hamiltonian Vlasov-Maxwell System
• Constructing current singularity in a 3D line-tied plasma
• Subcycling of particle orbits in variational, geometric electromagnetic particle-in-cell methods
• Gauge invariant canonical symplectic algorithms for real-time lattice strong-field quantum electrodynamics

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