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

Title: A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons

Abstract

A conservative scheme has been formulated and confirmed for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from the drift kinetic equation is used to time advance the electron density perturbation by using the perturbed mechanical flow calculated from the parallel vector potential, and the parallel vector potential is solved by using the perturbed canonical flow from the perturbed distribution function. In gyrokinetic particle simulations using this new scheme, the shear Alfvén wave dispersion relation in the shearless slab and continuum damping in the sheared cylinder have been recovered. The new scheme overcomes the stringent requirement in the conventional perturbative simulation method that perpendicular grid size needs to be as small as electron collisionless skin depth even for the long wavelength Alfvén waves. The new scheme also circumvents the problem in the typical approach that an unphysically large parallel electric field arises due to the inconsistency between electrostatic potential calculated from the perturbed density and vector potential calculated from the perturbed canonical flow. Lastly, the gyrokinetic particle simulations of the Alfvén waves in sheared cylinder have superior numerical properties compared with the fluid simulations, which suffer from numerical difficulties associated with singularmore » mode structures.« less

Authors:
 [1]; ORCiD logo [2]; ORCiD logo [3]
  1. Univ. of California, Irvine, CA (United States); Peking Univ., Beijing (China)
  2. Univ. of California, Irvine, CA (United States)
  3. Peking Univ., Beijing (China); Max Planck Inst. für Plasmaphysik, Garching (Germany)
Publication Date:
Research Org.:
Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC); Univ. of California, Oakland, CA (United States); UT-Battelle LLC/ORNL, Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR). Scientific Discovery through Advanced Computing (SciDAC); China Scholarship Council
OSTI Identifier:
1543850
Alternate Identifier(s):
OSTI ID: 1422452
Grant/Contract Number:  
AC02-05CH11231; AC05-00OR22725
Resource Type:
Accepted Manuscript
Journal Name:
Physics of Plasmas
Additional Journal Information:
Journal Volume: 25; Journal Issue: 2; Journal ID: ISSN 1070-664X
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Citation Formats

Bao, J., Lin, Z., and Lu, Z. X. A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons. United States: N. p., 2018. Web. doi:10.1063/1.5016432.
Bao, J., Lin, Z., & Lu, Z. X. A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons. United States. https://doi.org/10.1063/1.5016432
Bao, J., Lin, Z., and Lu, Z. X. Fri . "A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons". United States. https://doi.org/10.1063/1.5016432. https://www.osti.gov/servlets/purl/1543850.
@article{osti_1543850,
title = {A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons},
author = {Bao, J. and Lin, Z. and Lu, Z. X.},
abstractNote = {A conservative scheme has been formulated and confirmed for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from the drift kinetic equation is used to time advance the electron density perturbation by using the perturbed mechanical flow calculated from the parallel vector potential, and the parallel vector potential is solved by using the perturbed canonical flow from the perturbed distribution function. In gyrokinetic particle simulations using this new scheme, the shear Alfvén wave dispersion relation in the shearless slab and continuum damping in the sheared cylinder have been recovered. The new scheme overcomes the stringent requirement in the conventional perturbative simulation method that perpendicular grid size needs to be as small as electron collisionless skin depth even for the long wavelength Alfvén waves. The new scheme also circumvents the problem in the typical approach that an unphysically large parallel electric field arises due to the inconsistency between electrostatic potential calculated from the perturbed density and vector potential calculated from the perturbed canonical flow. Lastly, the gyrokinetic particle simulations of the Alfvén waves in sheared cylinder have superior numerical properties compared with the fluid simulations, which suffer from numerical difficulties associated with singular mode structures.},
doi = {10.1063/1.5016432},
journal = {Physics of Plasmas},
number = 2,
volume = 25,
place = {United States},
year = {Fri Feb 23 00:00:00 EST 2018},
month = {Fri Feb 23 00:00:00 EST 2018}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 11 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

An Eulerian gyrokinetic-Maxwell solver
journal, April 2003


Finite-β simulation of microinstabilities
journal, February 2014

  • Startsev, Edward A.; Lee, W. W.
  • Physics of Plasmas, Vol. 21, Issue 2
  • DOI: 10.1063/1.4863847

Nonlinear electromagnetic formulation for particle-in-cell simulation of lower hybrid waves in toroidal geometry
journal, June 2016

  • Bao, J.; Lin, Z.; Kuley, A.
  • Physics of Plasmas, Vol. 23, Issue 6
  • DOI: 10.1063/1.4952773

Shear-Alfvén waves in gyrokinetic plasmas
journal, October 2001

  • Lee, W. W.; Lewandowski, J. L. V.; Hahm, T. S.
  • Physics of Plasmas, Vol. 8, Issue 10
  • DOI: 10.1063/1.1400124

Electromagnetic formulation of global gyrokinetic particle simulation in toroidal geometry
journal, December 2009

  • Holod, I.; Zhang, W. L.; Xiao, Y.
  • Physics of Plasmas, Vol. 16, Issue 12
  • DOI: 10.1063/1.3273070

Plasmas
journal, December 1969


A fluid–kinetic hybrid electron model for electromagnetic simulations
journal, May 2001


Decay of MHD waves by phase mixing: II. The theta-pinch in cylindrical geometry
journal, June 1973

  • Grossmann, W.; Tataronis, J.
  • Zeitschrift für Physik A Hadrons and nuclei, Vol. 261, Issue 3
  • DOI: 10.1007/BF01391914

Decay of MHD waves by phase mixing: I. The sheet-pinch in plane geometry
journal, June 1973

  • Tataronis, J.; Grossmann, W.
  • Zeitschrift für Physik A Hadrons and nuclei, Vol. 261, Issue 3
  • DOI: 10.1007/BF01391913

Electromagnetic gyrokinetic PIC simulation with an adjustable control variates method
journal, July 2007


Nonlinear gyrokinetic theory for finite-beta plasmas
journal, January 1988

  • Hahm, T. S.; Lee, W. W.; Brizard, A.
  • Physics of Fluids, Vol. 31, Issue 7
  • DOI: 10.1063/1.866641

Microtearing Modes and Anomalous Transport in Tokamaks
journal, April 1980


Foundations of nonlinear gyrokinetic theory
journal, April 2007


Verification of electromagnetic fluid-kinetic hybrid electron model in global gyrokinetic particle simulation
journal, March 2013


Fluid electrons with kinetic closure for long wavelength energetic particles driven modes
journal, May 2011

  • Chen, Yang; Parker, Scott E.
  • Physics of Plasmas, Vol. 18, Issue 5
  • DOI: 10.1063/1.3567023

A conservative scheme of drift kinetic electrons for gyrokinetic simulation of kinetic-MHD processes in toroidal plasmas
journal, October 2017

  • Bao, J.; Liu, D.; Lin, Z.
  • Physics of Plasmas, Vol. 24, Issue 10
  • DOI: 10.1063/1.4995455

A fully nonlinear characteristic method for gyrokinetic simulation
journal, January 1993

  • Parker, S. E.; Lee, W. W.
  • Physics of Fluids B: Plasma Physics, Vol. 5, Issue 1
  • DOI: 10.1063/1.860870

Gyrokinetic particle simulation model
journal, September 1987


Verification of gyrokinetic particle simulation of current-driven instability in fusion plasmas. III. Collisionless tearing mode
journal, February 2016

  • Liu, Dongjian; Bao, Jian; Han, Tao
  • Physics of Plasmas, Vol. 23, Issue 2
  • DOI: 10.1063/1.4941094

Gyrokinetic simulation model for kinetic magnetohydrodynamic processes in magnetized plasmas
journal, January 2012


Turbulent Transport Reduction by Zonal Flows: Massively Parallel Simulations
journal, September 1998


Kinetic processes in plasma heating by resonant mode conversion of Alfvén wave
journal, January 1976

  • Hasegawa, Akira; Chen, Liu
  • Physics of Fluids, Vol. 19, Issue 12
  • DOI: 10.1063/1.861427

Partially Linearized Algorithms in Gyrokinetic Particle Simulation
journal, August 1993


A δf particle method for gyrokinetic simulations with kinetic electrons and electromagnetic perturbations
journal, August 2003


New variables for gyrokinetic electromagnetic simulations
journal, May 2014

  • Mishchenko, Alexey; Cole, Michael; Kleiber, Ralf
  • Physics of Plasmas, Vol. 21, Issue 5
  • DOI: 10.1063/1.4880560

Works referencing / citing this record:

Verification and validation of integrated simulation of energetic particles in fusion plasmas
journal, April 2019


Gyrokinetic simulations of nonlinear interactions between magnetic islands and microturbulence
journal, September 2019

  • Fang, Kaisheng; Bao, Jian; Lin, Zhihong
  • Plasma Science and Technology, Vol. 21, Issue 11
  • DOI: 10.1088/2058-6272/ab3c83