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

Title: A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma

Journal Article · · Journal of Computational Physics
ORCiD logo [1];  [1];  [1];  [2];  [3]
  1. Princeton Univ., Princeton, NJ (United States). Princeton Plasma Physics Lab.
  2. National Fusion Research Institute (Republic of Korea)
  3. Univ. of Colorado, Boulder, CO (United States)

In order to enable kinetic simulation of non-thermal edge plasmas at a reduced computational cost, a new hybrid-Lagrangian δf scheme has been developed that utilizes the phase space grid in addition to the usual marker particles, taking advantage of the computational strengths from both sides. The new scheme splits the particle distribution function of a kinetic equation into two parts. Marker particles contain the fast space-time varying, δf, part of the distribution function and the coarse-grained phase-space grid contains the slow space-time varying part. The coarse-grained phase-space grid reduces the memory-requirement and the computing cost, while the marker particles provide scalable computing ability for the fine-grained physics. Weights of the marker particles are determined by a direct weight evolution equation instead of the differential form weight evolution equations that the conventional delta-f schemes use. The particle weight can be slowly transferred to the phase space grid, thereby reducing the growth of the particle weights. The non-Lagrangian part of the kinetic equation – e.g., collision operation, ionization, charge exchange, heat-source, radiative cooling, and others – can be operated directly on the phase space grid. Deviation of the particle distribution function on the velocity grid from a Maxwellian distribution function – driven by ionization, charge exchange and wall loss – is allowed to be arbitrarily large. In conclusion, the numerical scheme is implemented in the gyrokinetic particle code XGC1, which specializes in simulating the tokamak edge plasma that crosses the magnetic separatrix and is in contact with the material wall.

Research Organization:
Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
Sponsoring Organization:
USDOE Office of Science (SC)
Contributing Organization:
Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); National Fusion Research Institute, Republic of Korea; University of Colorado Boulder, USA
Grant/Contract Number:
AC02-09CH11466; SC000801; AC05-00OR22725
OSTI ID:
1254687
Alternate ID(s):
OSTI ID: 1325288
Report Number(s):
PPPL-5212; PII: S0021999116300274
Journal Information:
Journal of Computational Physics, Vol. 315, Issue C; ISSN 0021-9991
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 58 works
Citation information provided by
Web of Science

References (15)

Gyrokinetic particle simulation model journal September 1987
Large-scale gyrokinetic turbulence simulations: Effects of profile variation journal May 1999
Electron temperature gradient driven turbulence journal May 2000
Electromagnetic gyrokinetic δf particle-in-cell turbulence simulation with realistic equilibrium profiles and geometry journal January 2007
A fluid–kinetic hybrid electron model for electromagnetic simulations journal May 2001
The Semi-Lagrangian Method for the Numerical Resolution of the Vlasov Equation journal March 1999
Beyond scale separation in gyrokinetic turbulence journal August 2007
Numerical simulation of plasmas with periodic smoothing in phase space journal February 1972
The particle-continuum method: an algorithmic unification of particle-in-cell and continuum methods journal December 2004
Full-f gyrokinetic particle simulation of centrally heated global ITG turbulence from magnetic axis to edge pedestal top in a realistic tokamak geometry journal September 2009
The role of dissipation in the theory and simulations of homogeneous plasma turbulence, and resolution of the entropy paradox journal October 1994
Generalized weighting scheme for δ f particle‐simulation method journal April 1994
Collisional δf method journal October 1997
A new bdelta f method for neoclassical transport studies journal August 1999
A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma journal June 2016

Cited By (27)

Pedestal and edge electrostatic turbulence characteristics from an XGC1 gyrokinetic simulation journal August 2017
A tight-coupling scheme sharing minimum information across a spatial interface between gyrokinetic turbulence codes journal July 2018
Gyrokinetic projection of the divertor heat-flux width from present tokamaks to ITER journal August 2017
Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code journal May 2017
Invited Review Article: Gas puff imaging diagnostics of edge plasma turbulence in magnetic fusion devices journal April 2017
XGC software June 2018
Cross-verification of the global gyrokinetic codes GENE and XGC journal June 2018
Gyrokinetic simulation study of magnetic island effects on neoclassical physics and micro-instabilities in a realistic KSTAR plasma journal May 2018
Gyroaveraging operations using adaptive matrix operators journal May 2018
Analysis Of Equilibrium And Turbulent Fluxes Across The Separatrix In A Gyrokinetic Simulation [Supplementary Data]
  • Charidakos, Ioannis Keramidas; Myra, James; Parker, Scott
  • Zenodo-Supplementary information for journal article at DOI: 10.1063/1.5037723, 15 files (333.7 kB) https://doi.org/10.5281/zenodo.1230040
dataset April 2018
Pressure balance in a lower collisionality, attached tokamak scrape-off layer journal July 2019
Verification of the global gyrokinetic stellarator code XGC-S for linear ion temperature gradient driven modes journal August 2019
Full- f version of GENE for turbulence in open-field-line systems journal June 2018
Integrated Tokamak modeling: When physics informs engineering and research planning journal May 2018
Comparative collisionless alpha particle confinement in stellarator reactors with the XGC gyrokinetic code journal March 2019
Improved linearized model collision operator for the highly collisional regime journal October 2019
Development of a Gyrokinetic Particle-in-Cell Code for Whole-Volume Modeling of Stellarators journal May 2019
A fast low-to-high confinement mode bifurcation dynamics in the boundary-plasma gyrokinetic code XGC1 journal May 2018
Blob-hole correlation model for edge turbulence and comparisons with NSTX gas puff imaging data journal June 2018
Analysis of equilibrium and turbulent fluxes across the separatrix in a gyrokinetic simulation journal July 2018
Analysis Of Equilibrium And Turbulent Fluxes Across The Separatrix In A Gyrokinetic Simulation [Supplementary Data]
  • Charidakos, Ioannis Keramidas; Myra, James; Parker, Scott
  • Zenodo-Supplementary information for journal article at DOI: 10.1063/1.5037723, 15 files (333.7 kB) https://doi.org/10.5281/zenodo.1230039
dataset April 2018
Gyrokinetic projection of the divertor heat-flux width from present tokamaks to ITER text January 2017
Gyrokinetic continuum simulation of turbulence in a straight open-field-line plasma text January 2017
Gyroaveraging operations using adaptive matrix operators text January 2018
Analysis of equilibrium and turbulent fluxes across the separatrix in a gyrokinetic simulation text January 2018
A tight-coupling scheme sharing minimum information across a spatial interface between gyrokinetic turbulence codes text January 2018
Electromagnetic full-$f$ gyrokinetics in the tokamak edge with discontinuous Galerkin methods text January 2019

Similar Records

A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma
Journal Article · Wed Jun 15 00:00:00 EDT 2016 · Journal of Computational Physics · OSTI ID:1254687

Partnership for Edge Physics Simulation (EPSI)
Technical Report · Wed Feb 11 00:00:00 EST 2015 · OSTI ID:1254687

A Fokker-Planck-Landau collision equation solver on two-dimensional velocity grid and its application to particle-in-cell simulation
Journal Article · Sat Mar 15 00:00:00 EDT 2014 · Physics of Plasmas · OSTI ID:1254687