skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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

Journal Article · · Plasma Science and Technology
 [1];  [2]
  1. University of Science and Technology of China, Hefei (China). School of Physical Sciences
  2. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)

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 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-1, agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown 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 at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

Research Organization:
Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Sponsoring Organization:
USDOE; National MCF Energy R&D Program; National Key Research and Development Program; National Natural Science Foundation of China (NSFC)
Grant/Contract Number:
AC02-09CH11466; 2018YFE0304100; 2016YFA0400600; 2016YFA0400601; 016YFA0400602; 11905220; 11805273
OSTI ID:
1809092
Journal Information:
Plasma Science and Technology, Vol. 23, Issue 5; ISSN 1009-0630
Publisher:
IOPScienceCopyright Statement
Country of Publication:
United States
Language:
English

References (128)

Degenerate variational integrators for magnetic field line flow and guiding center trajectories journal May 2018
Gyrocenter-gauge kinetic theory journal November 2000
Symplectic integrators with adaptive time step applied to runaway electron dynamics journal January 2019
A general concurrent algorithm for plasma particle-in-cell simulation codes journal December 1989
Volume-preserving maps, source-free systems and their local structures journal April 2006
Linear gyrokinetic analysis of a DIII-D H-mode pedestal near the ideal ballooning threshold journal September 2012
Gyrokinetic theory for arbitrary wavelength electromagnetic modes in tokamaks journal April 1998
Comment on “Hamiltonian splitting for the Vlasov–Maxwell equations” journal September 2015
Hamiltonian particle-in-cell methods for Vlasov-Maxwell equations journal September 2016
Explicit symplectic methods for solving charged particle trajectories journal May 2017
Simulation of plasma turbulence in scrape-off layer conditions: the GBS code, simulation results and code validation journal November 2012
Hamiltonian gyrokinetic Vlasov–Maxwell system journal September 2015
The symplecticity of multi-step methods journal February 1993
Nonphysical noises and instabilities in plasma simulation due to a spatial grid journal December 1972
Volume-preserving algorithm for secular relativistic dynamics of charged particles journal April 2015
Gauge invariant canonical symplectic algorithms for real-time lattice strong-field quantum electrodynamics journal February 2021
Implicit time integration for plasma simulation journal August 1982
Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov–Maxwell equations journal December 2015
Global gyrokinetic simulations of the H-mode tokamak edge pedestal journal May 2013
Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media journal May 1966
Compressed ion temperature gradient turbulence in diverted tokamak edge journal May 2009
Numerical simulations of one laser-plasma model based on Poisson structure journal March 2020
Machine learning and serving of discrete field theories journal November 2020
KAM theorem of symplectic algorithms for Hamiltonian systems journal September 1999
Formation of current singularity in a topologically constrained plasma journal February 2016
A δf particle method for gyrokinetic simulations with kinetic electrons and electromagnetic perturbations journal August 2003
Nonlinear reduced Braginskii equations with ion thermal dynamics in toroidal plasma journal June 1997
Structure and structure-preserving algorithms for plasma physics journal April 2017
Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems journal November 2015
3D multispecies nonlinear perturbative particle simulations of collective processes in intense particle beams for heavy ion fusion
  • Qin, Hong; Davidson, Ronald C.; Lee, W. Wei-li
  • Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, Vol. 464, Issue 1-3 https://doi.org/10.1016/S0168-9002(01)00713-6
journal May 2001
Volume-preserving algorithms for charged particle dynamics journal January 2015
Exact charge conservation scheme for Particle-in-Cell simulation with an arbitrary form-factor journal April 2001
Structure-preserving electromagnetic–kinetic simulations of lower hybrid-wave injection and current drive journal November 2020
LIGKA: A linear gyrokinetic code for the description of background kinetic and fast particle effects on the MHD stability in tokamaks journal September 2007
Variational symplectic algorithm for guiding center dynamics and its application in tokamak geometry journal April 2009
Toroidal regularization of the guiding center Lagrangian journal November 2017
First results from Alcator‐C‐MOD* journal May 1994
Rigorous charge conservation for local electromagnetic field solvers journal March 1992
Geometric Particle-in-Cell Simulations of the Vlasov--Maxwell System in Curvilinear Coordinates journal January 2021
Electromagnetic gyrokinetic δf particle-in-cell turbulence simulation with realistic equilibrium profiles and geometry journal January 2007
Gyrokinetic turbulence simulations at high plasma beta journal October 2008
Can time be a discrete dynamical variable? journal March 1983
Discrete mechanics and variational integrators journal May 2001
Exact charge-conserving scatter–gather algorithm for particle-in-cell simulations on unstructured grids: A geometric perspective journal September 2015
Runge-kutta schemes for Hamiltonian systems journal December 1988
Conservative semi-Lagrangian schemes for Vlasov equations journal March 2010
Difference Discrete Variational Principles, Euler–Lagrange Cohomology and Symplectic, Multisymplectic Structures III: Application to Symplectic and Multisymplectic Algorithms journal March 2002
GTC Simulation of Ideal Ballooning Mode in Tokamak Plasmas journal June 2013
Multisymplecticity of the centred box discretization for hamiltonian PDEs with m ≥ 2 space dimensions journal November 2002
Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme journal August 2012
A spectral canonical electrostatic algorithm journal February 2016
Higher order volume-preserving schemes for charged particle dynamics journal January 2016
Explicit K -symplectic algorithms for charged particle dynamics journal February 2017
Canonical symplectic structure and structure-preserving geometric algorithms for Schrödinger–Maxwell systems journal November 2017
Canonicalization and symplectic simulation of the gyrocenter dynamics in time-independent magnetic fields journal March 2014
Development of variational guiding center algorithms for parallel calculations in experimental magnetic equilibria journal April 2015
Nonlinear gyrokinetic equations for tokamak microturbulence journal September 1988
Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems journal November 2016
Hamiltonian time integrators for Vlasov-Maxwell equations journal December 2015
Observation of Self-Generated Flows in Tokamak Plasmas with Lower-Hybrid-Driven Current journal January 2009
The Sunway TaihuLight supercomputer: system and applications journal June 2016
The geometric theory of charge conservation in particle-in-cell simulations journal May 2020
Variational Symplectic Integrator for Long-Time Simulations of the Guiding-Center Motion of Charged Particles in General Magnetic Fields journal January 2008
Energy-conserving time propagation for a structure-preserving particle-in-cell Vlasov–Maxwell solver journal January 2021
Particle-in-cell charged-particle simulations, plus Monte Carlo collisions with neutral atoms, PIC-MCC journal April 1991
Linear gyrokinetic theory for kinetic magnetohydrodynamic eigenmodes in tokamak plasmas journal June 1999
Explicit high-order gauge-independent symplectic algorithms for relativistic charged particle dynamics journal August 2019
A symplectic integration algorithm for separable Hamiltonian functions journal January 1991
The virtual particle electromagnetic particle-mesh method journal May 1991
Comment on “Symplectic integration of magnetic systems”: A proof that the Boris algorithm is not variational journal November 2015
Performance and optimization of direct implicit particle simulation journal March 1989
Subcycling of particle orbits in variational, geometric electromagnetic particle-in-cell methods journal September 2020
Mesh refinement for particle-in-cell plasma simulations: Applications to and benefits for heavy ion fusion journal October 2002
A study of the numerical heating in electrostatic particle simulations journal April 1994
Symplectic integration of Hamiltonian systems journal May 1990
Gauge-free electromagnetic gyrokinetic theory journal June 2019
QUICKPIC: A highly efficient particle-in-cell code for modeling wakefield acceleration in plasmas journal September 2006
Global Gyrokinetic Simulation of Tokamak Edge Pedestal Instabilities journal November 2012
Three-dimensional multispecies nonlinear perturbative particle simulations of collective processes in intense particle beams journal August 2000
On plasma rotation induced by waves in tokamaks journal October 2013
Variational integration for ideal magnetohydrodynamics with built-in advection equations journal October 2014
Gyrokinetic perpendicular dynamics journal May 1999
Field theory and a structure-preserving geometric particle-in-cell algorithm for drift wave instability and turbulence journal September 2019
A lattice Maxwell system with discrete space–time symmetry and local energy–momentum conservation journal February 2019
Difference equations and conservation laws journal March 1987
The Hamiltonian structure of the Maxwell-Vlasov equations journal March 1982
Particle simulation of plasmas journal April 1983
Finite-dimensional collisionless kinetic theory journal March 2017
Nonlinear Simulations of Peeling-Ballooning Modes with Anomalous Electron Viscosity and their Role in Edge Localized Mode Crashes journal October 2010
Integrable discrete-time systems and difference operators journal January 1988
GEMPIC: geometric electromagnetic particle-in-cell methods journal July 2017
Local energy conservation law for a spatially-discretized Hamiltonian Vlasov-Maxwell system journal June 2017
Variational approach to low-frequency kinetic-MHD in the current coupling scheme journal March 2017
Multi-scale particle-in-cell plasma simulation journal September 1991
Volume-preserving algorithms for source-free dynamical systems journal October 1995
Hamiltonian splitting for the Vlasov–Maxwell equations journal February 2015
A variational multi-symplectic particle-in-cell algorithm with smoothing functions for the Vlasov-Maxwell system journal October 2013
Comments on: The Maxwell-Vlasov equations as a continuous hamiltonian system journal November 1981
A charge- and energy-conserving implicit, electrostatic particle-in-cell algorithm on mapped computational meshes journal January 2013
An energy- and charge-conserving, implicit, electrostatic particle-in-cell algorithm journal August 2011
The Maxwell-Vlasov equations as a continuous hamiltonian system journal December 1980
Constructing Current Singularity in a 3D Line-tied Plasma journal December 2017
Observations of counter-current toroidal rotation in Alcator C-Mod LHCD plasmas journal January 2009
An Explicit Symplectic Integration Scheme for Plasma Simulations journal July 1993
Direct implicit large time-step particle simulation of plasmas journal July 1983
VORPAL: a versatile plasma simulation code journal May 2004
A Can0nical Integrati0n Technique journal August 1983
Variational formulation of macro-particle plasma simulation algorithms journal May 2014
Simulations of relativistic quantum plasmas using real-time lattice scalar QED journal May 2018
Why is Boris algorithm so good? journal August 2013
Pullback transformations in gyrokinetic theory journal March 2004
Structure-preserving geometric particle-in-cell methods for Vlasov-Maxwell systems journal September 2018
Nonlinear saturation of kinetic ballooning modes by zonal fields in toroidal plasmas journal January 2019
Kinetic-ballooning-mode theory in general geometry journal November 1980
Effective-action approach to wave propagation in scalar QED plasmas journal July 2016
Explicit high-order symplectic integrators for charged particles in general electromagnetic fields journal December 2016
Turbulent Transport Reduction by Zonal Flows: Massively Parallel Simulations journal September 1998
A family of new explicit, revertible, volume-preserving numerical schemes for the system of Lorentz force journal December 2016
Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm journal October 2016
Nonlinear gyrokinetic Maxwell-Vlasov equations using magnetic co-ordinates journal June 1989
Construction of higher order symplectic integrators journal November 1990
Variational formulation of particle algorithms for kinetic plasma simulations journal July 2013
Solving the Vlasov–Maxwell equations using Hamiltonian splitting journal November 2019
Geometric gyrokinetic theory for edge plasmas journal May 2007
H mode confinement in Alcator C-Mod journal June 1997
Gyrokinetic field theory journal February 2000
Structure-preserving vs. standard particle-in-cell methods: The case of an electron hybrid model journal February 2020
3D nonlinear perturbative particle simulations of two-stream collective processes in intense particle beams journal August 2000