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

Explicit Non-Canonical Symplectic Algorithms for Charged Particle Dynamics

The High-Order Symplectic Finite-Difference Time-Domain Scheme

SYMPLECTIC PSEUDOSPECTRAL TIME-DOMAIN SCHEME FOR SOLVING TIME-DEPENDENT SCHRODINGER EQUATION

Simulating Maxwell–Schrödinger Equations by High-Order Symplectic FDTD Algorithm

A novel class of explicit energy-preserving splitting methods for charged-particle dynamics

A Symplectic Algorithm For The Stability Analysis Of Nonlinear Parametric Excited Systems

Splitting Methods for Semi-Classical Hamiltonian Dynamics of Charge Transfer in Nonlinear Lattices

Continuous-stage Adapted Exponential Methods for Charged-Particle Dynamics with Arbitrary Magnetic Fields

Novel semi-explicit symplectic schemes for nonseparable stochastic Hamiltonian systems

Novel structure-preserving schemes for stochastic Klein–Gordon–Schrödinger equations with additive noise

Novel structure-preserving schemes for stochastic Klein--Gordon--Schrödinger equations with additive noise

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

High order explicit Lorentz invariant volume-preserving algorithms for relativistic dynamics of charged particles

Explicit symplectic integrators with adaptive time steps in curved spacetimes

Structure-preserving quantum algorithms for linear and nonlinear Hamiltonian systems

Hamiltonian formulation and symplectic split-operator schemes for time-dependent density-functional-theory equations of electron dynamics in molecules

High order explicit symplectic integrators for the Discrete Non Linear Schrödinger equation

Applying explicit symplectic integrator to study chaos of charged particles around magnetized Kerr black hole

Local Structure-Preserving Relaxation Method for Charged Systems on Unstructured Meshes

Advanced fuel fusion, phase space engineering, and structure-preserving geometric algorithms

A Symplectic Approximation with Nonlinear Stability and Convergence Analysis for Efficiently Solving Semi-Linear Klein–Gordon Equations