DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE

Abstract

we report a multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.

Authors:
ORCiD logo [1];  [2];  [2];  [3];  [4];  [2]; ORCiD logo [2];  [2];  [2];  [5]
  1. Univ. of California, Los Angeles, CA (United States)
  2. Max Planck Institute for Plasma Physics, Garching (Germany)
  3. Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States)
  4. Univ. of Texas, Austin, TX (United States)
  5. Max Planck Institute for Plasma Physics, Garching (Germany); Univ. of Texas, Austin, TX (United States)
Publication Date:
Research Org.:
Univ. of California, Los Angeles, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC); Euratom Research and Training Programme
OSTI Identifier:
1852090
Grant/Contract Number:  
SC0016073; 633053
Resource Type:
Accepted Manuscript
Journal Name:
Computer Physics Communications
Additional Journal Information:
Journal Volume: 255; Journal Issue: C; Journal ID: ISSN 0010-4655
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; gyrokinetics; collisions; GENE

Citation Formats

Crandall, P., Jarema, D., Doerk, H., Pan, Q., Merlo, G., Görler, T., Navarro, A. Bañón, Told, D., Maurer, M., and Jenko, F. Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE. United States: N. p., 2020. Web. doi:10.1016/j.cpc.2020.107360.
Crandall, P., Jarema, D., Doerk, H., Pan, Q., Merlo, G., Görler, T., Navarro, A. Bañón, Told, D., Maurer, M., & Jenko, F. Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE. United States. https://doi.org/10.1016/j.cpc.2020.107360
Crandall, P., Jarema, D., Doerk, H., Pan, Q., Merlo, G., Görler, T., Navarro, A. Bañón, Told, D., Maurer, M., and Jenko, F. Thu . "Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE". United States. https://doi.org/10.1016/j.cpc.2020.107360. https://www.osti.gov/servlets/purl/1852090.
@article{osti_1852090,
title = {Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE},
author = {Crandall, P. and Jarema, D. and Doerk, H. and Pan, Q. and Merlo, G. and Görler, T. and Navarro, A. Bañón and Told, D. and Maurer, M. and Jenko, F.},
abstractNote = {we report a multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.},
doi = {10.1016/j.cpc.2020.107360},
journal = {Computer Physics Communications},
number = C,
volume = 255,
place = {United States},
year = {Thu May 14 00:00:00 EDT 2020},
month = {Thu May 14 00:00:00 EDT 2020}
}

Works referenced in this record:

Approximate Fokker–Planck collision operator for transport theory applications
journal, January 1976

  • Hirshman, S. P.; Sigmar, D. J.
  • Physics of Fluids, Vol. 19, Issue 10
  • DOI: 10.1063/1.861356

Implications of advanced collision operators for gyrokinetic simulation
journal, February 2017


Tokamak profile prediction using direct gyrokinetic and neoclassical simulation
journal, June 2009

  • Candy, J.; Holland, C.; Waltz, R. E.
  • Physics of Plasmas, Vol. 16, Issue 6
  • DOI: 10.1063/1.3167820

SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems
journal, September 2005

  • Hernandez, Vicente; Roman, Jose E.; Vidal, Vicente
  • ACM Transactions on Mathematical Software, Vol. 31, Issue 3
  • DOI: 10.1145/1089014.1089019

The global version of the gyrokinetic turbulence code GENE
journal, August 2011

  • Görler, T.; Lapillonne, X.; Brunner, S.
  • Journal of Computational Physics, Vol. 230, Issue 18
  • DOI: 10.1016/j.jcp.2011.05.034

Progress with the COGENT Edge Kinetic Code: Collision Operator Options
journal, June 2012

  • Dorf, M. A.; Cohen, R. H.; Compton, J. C.
  • Contributions to Plasma Physics, Vol. 52, Issue 5-6
  • DOI: 10.1002/ctpp.201210042

Towards optimal explicit time-stepping schemes for the gyrokinetic equations
journal, July 2014


Block-structured grids for Eulerian gyrokinetic simulations
journal, January 2016


Energetically consistent collisional gyrokinetics
journal, October 2015

  • Burby, J. W.; Brizard, A. J.; Qin, H.
  • Physics of Plasmas, Vol. 22, Issue 10
  • DOI: 10.1063/1.4935124

Implementation of a gyrokinetic collision operator with an implicit time integration scheme and its computational performance
journal, February 2019


Improved collision operator for plasma kinetic simulations with multi-species ions and electrons
journal, December 2015

  • Nakata, Motoki; Nunami, Masanori; Watanabe, Tomo-Hiko
  • Computer Physics Communications, Vol. 197
  • DOI: 10.1016/j.cpc.2015.08.007

Electron temperature gradient driven turbulence
journal, May 2000

  • Jenko, F.; Dorland, W.; Kotschenreuther, M.
  • Physics of Plasmas, Vol. 7, Issue 5
  • DOI: 10.1063/1.874014

Linearized Coulomb Collision Operator for Simulation of Interpenetrating Plasma Streams
journal, May 2019

  • Dimits, A. M.; Joseph, I.; Banks, J. W.
  • IEEE Transactions on Plasma Science, Vol. 47, Issue 5
  • DOI: 10.1109/TPS.2019.2897790

Gyrokinetic linearized Landau collision operator
journal, January 2013


A high-accuracy Eulerian gyrokinetic solver for collisional plasmas
journal, November 2016


Conservative discretization of the Landau collision integral
journal, March 2017

  • Hirvijoki, E.; Adams, M. F.
  • Physics of Plasmas, Vol. 24, Issue 3
  • DOI: 10.1063/1.4979122

A multi-species collisional operator for full-F global gyrokinetics codes: Numerical aspects and verification with the GYSELA code
journal, January 2019


Differential formulation of the gyrokinetic Landau operator
journal, January 2017

  • Hirvijoki, Eero; Brizard, Alain J.; Pfefferlé, David
  • Journal of Plasma Physics, Vol. 83, Issue 1
  • DOI: 10.1017/S0022377816001203

Synergy between ion temperature gradient turbulence and neoclassical processes in global gyrokinetic particle-in-cell simulations
journal, April 2012

  • Vernay, T.; Brunner, S.; Villard, L.
  • Physics of Plasmas, Vol. 19, Issue 4
  • DOI: 10.1063/1.3699189

Gyrokinetic Landau collision operator in conservative form
journal, February 2019


A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma
journal, June 2016


Enhanced stabilisation of trapped electron modes by collisional energy scattering in tokamaks
journal, June 2015

  • Manas, P.; Camenen, Y.; Benkadda, S.
  • Physics of Plasmas, Vol. 22, Issue 6
  • DOI: 10.1063/1.4922754

A multi-species collisional operator for full-F gyrokinetics
journal, December 2015

  • Estève, D.; Garbet, X.; Sarazin, Y.
  • Physics of Plasmas, Vol. 22, Issue 12
  • DOI: 10.1063/1.4937373

Kinetic calculation of neoclassical transport including self-consistent electron and impurity dynamics
journal, July 2008


A 5D gyrokinetic full- f global semi-Lagrangian code for flux-driven ion turbulence simulations
journal, October 2016


Noncircular, finite aspect ratio, local equilibrium model
journal, April 1998

  • Miller, R. L.; Chu, M. S.; Greene, J. M.
  • Physics of Plasmas, Vol. 5, Issue 4
  • DOI: 10.1063/1.872666

Linearized model Fokker–Planck collision operators for gyrokinetic simulations. I. Theory
journal, December 2008

  • Abel, I. G.; Barnes, M.; Cowley, S. C.
  • Physics of Plasmas, Vol. 15, Issue 12
  • DOI: 10.1063/1.3046067

Theory of plasma transport in toroidal confinement systems
journal, April 1976


Block-structured grids in full velocity space for Eulerian gyrokinetic simulations
journal, June 2017


Full linearized Fokker–Planck collisions in neoclassical transport simulations
journal, December 2011


Linearized model Fokker–Planck collision operators for gyrokinetic simulations. II. Numerical implementation and tests
journal, July 2009

  • Barnes, M.; Abel, I. G.; Dorland, W.
  • Physics of Plasmas, Vol. 16, Issue 7
  • DOI: 10.1063/1.3155085