Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

First implementation of gyrokinetic exact linearized Landau collision operator and comparison with models

Journal Article · · Physics of Plasmas
DOI:https://doi.org/10.1063/1.5143374· OSTI ID:1799694
 [1];  [2];  [3]
  1. Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States). Plasma Science and Fusion Center; Massachusetts Institute of Technology
  2. Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States). Plasma Science and Fusion Center
  3. Univ. of California, Los Angeles, CA (United States)

Gyrokinetic simulations are fundamental to understanding and predicting turbulent transport in magnetically confined fusion plasmas. Previous simulations have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius (FLR) effects. We have implemented the linearized Fokker–Planck collision operator with exact field-particle terms and full FLR effects in a gyrokinetic code (GENE). The new operator, referred to as “exact” in this paper, allows the accuracy of model collision operators to be assessed. The conservative Landau form is implemented because its symmetry underlies the conservation laws and the H-theorem, and enables numerical methods to preserve this conservation, independent of resolution. The implementation utilizes the finite-volume method recently employed to discretize the Sugama collision model in GENE, allowing direct comparison between the two operators. Results show that the Sugama model appears accurate for the growth rates of trapped electron modes (TEMs) driven only by density gradients, but appreciably underestimates the growth rates as the collisionality and electron temperature gradient increase. The TEM turbulent fluxes near the nonlinear threshold using the exact operator are similar to the Sugama model for the ne = d ln Te/d ln ne = 0 case, but substantially larger than the Sugama model for the ne = 1 case. The FLR effects reduce the growth rates increasingly with wavenumber deepening a “valley” at the intermediate binormal wavenumber as the unstable mode extends from the TEM regime to the electron temperature gradient instability regime. Application to the Hinton–Rosenbluth problem shows that zonal flows decay faster as the radial wavenumber increases and the exact operator yields weaker decay rates.

Research Organization:
Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States); Univ. of California, Oakland, CA (United States)
Sponsoring Organization:
USDOE Office of Science (SC)
Grant/Contract Number:
FC02-08ER54966; AC02-05CH11231
OSTI ID:
1799694
Alternate ID(s):
OSTI ID: 1615074
Journal Information:
Physics of Plasmas, Journal Name: Physics of Plasmas Journal Issue: 4 Vol. 27; ISSN 1070-664X
Publisher:
American Institute of Physics (AIP)Copyright Statement
Country of Publication:
United States
Language:
English

References (41)

Comparison of initial value and eigenvalue codes for kinetic toroidal plasma instabilities journal August 1995
Remarks on algorithm 006: An adaptive algorithm for numerical integration over an N-dimensional rectangular region journal December 1980
Towards optimal explicit time-stepping schemes for the gyrokinetic equations journal July 2014
Improved collision operator for plasma kinetic simulations with multi-species ions and electrons journal December 2015
New conservative gyrokinetic full-f Vlasov code and its comparison to gyrokinetic δf particle-in-cell code journal September 2007
New velocity-space discretization for continuum kinetic calculations and Fokker–Planck collisions journal June 2013
A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma journal June 2016
Nonlinear gyrokinetic Coulomb collision operator journal November 2019
Role of trapped electron mode turbulence in internal transport barrier control in the Alcator C-Mod Tokamak journal May 2004
A guiding-center Fokker–Planck collision operator for nonuniform magnetic fields journal September 2004
Gyrokinetic simulation of collisionless trapped-electron mode turbulence journal July 2005
Linearized model Fokker–Planck collision operators for gyrokinetic simulations. I. Theory journal December 2008
Role of zonal flows in trapped electron mode turbulence through nonlinear gyrokinetic particle and continuum simulation journal May 2009
Linearized model Fokker–Planck collision operators for gyrokinetic simulations. II. Numerical implementation and tests journal July 2009
Linearized model collision operators for multiple ion species plasmas and gyrokinetic entropy balance equations journal November 2009
Drift-kinetic simulation of neoclassical transport with impurities in tokamaks journal February 2010
A Fokker-Planck-Landau collision equation solver on two-dimensional velocity grid and its application to particle-in-cell simulation journal March 2014
Energetically consistent collisional gyrokinetics journal October 2015
Role of density gradient driven trapped electron mode turbulence in the H-mode inner core with electron heating journal May 2016
Gyrokinetic energy conservation and Poisson‐bracket formulation journal July 1989
Numerical simulation of ion‐temperature‐gradient‐driven modes journal March 1991
Approximate Fokker–Planck collision operator for transport theory applications journal January 1976
Linearized gyro-kinetic equation with collisions journal January 1977
Nonlinear gyrokinetic equations for low-frequency electromagnetic waves in general plasma equilibria journal January 1982
Nonlinear gyrokinetic equations journal January 1983
Collisional effects on kinetic electromagnetic modes and associated quasilinear transport journal January 1987
Gyrokinetic particle simulation of neoclassical transport journal August 1995
Electron temperature gradient driven turbulence journal May 2000
Dynamics of axisymmetric and poloidal flows in tokamaks journal January 1999
Zonal flows in plasma—a review journal April 2005
Alternate form of model like particle collision operator journal May 2009
Full linearized Fokker–Planck collisions in neoclassical transport simulations journal December 2011
Implications of advanced collision operators for gyrokinetic simulation journal February 2017
Uncertainty estimation and a stopping rule in nonlinear gyrokinetic simulations journal November 2016
Fokker-Planck Equation for an Inverse-Square Force journal July 1957
Collision operators for partially linearized particle simulation codes journal January 1994
Gyrokinetic linearized Landau collision operator journal January 2013
Gyrokinetic Landau collision operator in conservative form journal February 2019
Gyrokinetic Fokker-Planck Collision Operator journal May 2011
Foundations of nonlinear gyrokinetic theory journal April 2007
The Exponentially Convergent Trapezoidal Rule journal January 2014