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

Variational Symplectic Algorithm for Guiding Center Dynamics and Its Application in Tokamak Geometry

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

Symplectic Simulations of Radial Diffusion of Fast Alpha Particles in the Presence of Low-Frequency Modes in Rippled Tokamaks

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

Dynamic Optimization Of Tokamak Plasmas Via Control Parameterization And The Time-Scaling Transformation

Explicit Runge–Kutta Integrator with Hamiltonian Correction for Long-Time Simulations of Guiding-Center Orbit in Tokamak Configurations

Design and Simulation of Extremum-Seeking Open-Loop Optimal Control of Current Profile in the DIII-D Tokamak

Dynamic Optimization of Open-Loop Input Signals for Ramp-Up Current Profiles in Tokamak Plasmas

Finite Element Approximation and Input Parameterization for the Optimal Control of Current Profiles in Tokamak Plasmas

Optimal Tracking for a Divergent-Type Parabolic PDE System in Current Profile Control

Simserver Simulation of a Model-Based Current Profile Controller in the DIII-D Plasma Control System

High-Precision Nonsingular Integrator of Guiding-Center Orbit Close to Magnetic Axis in Tokamak Equilibrium Configuration

Numerical simulation of magnetohydrodynamic instability with variational principle in the tokamak

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

Canonical Hamiltonian Guiding Center Dynamics and Its Intrinsic Magnetic Moment

Targeting realistic geometry in Tokamak code Gysela

A mimetic finite difference based quasi-static magnetohydrodynamic solver for force-free plasmas in tokamak disruptions

Variational formulation of macro-particle plasma simulation algorithms

Emulation Techniques for Scenario and Classical Control Design of Tokamak Plasmas

Linear Gyrokinetic Theory and Computation of the Gyrocenter Motion Based on the Exact Canonical Variables for Axisymmetric Tokamaks

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