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

Continuum electromagnetic gyrokinetic simulations of turbulence in the tokamak scrape-off layer and laboratory devices

Journal Article · · Physics of Plasmas
DOI:https://doi.org/10.1063/1.5141157· OSTI ID:1668050
 [1];  [2];  [3];  [4];  [1];  [5]
  1. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
  2. Princeton Univ., NJ (United States)
  3. General Atomics, San Diego, CA (United States); Univ. of Texas, Austin, TX (United States)
  4. MIT Plasma Science and Fusion Center, Cambridge, MA (United States); Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
  5. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Here, we present algorithms and results from Gkeyll, a full-f continuum, electromagnetic gyrokinetic code, designed to study turbulence in the edge region of fusion devices. The edge is computationally very challenging, requiring robust algorithms that can handle large-amplitude fluctuations and stable interactions with plasma sheaths. We present an energy-conserving high-order discontinuous Galerkin scheme that solves gyrokinetic equations in Hamiltonian form. Efficiency is improved by a careful choice of basis functions and automatically generated computation kernels. Previous verification tests were performed in the straight-field-line large plasma device [Shi et al., J. Plasma Phys. 83, 905830304 (2017)] and the Texas Helimak, a simple magnetized torus [Bernard et al., Phys. Plasmas 26, 042301 (2019)], including the effect of end-plate biasing on turbulence. Results for the scrape-off layer for NSTX parameters with a model helical magnetic geometry with bad curvature have been obtained [Shi et al., Phys. Plasmas 26, 012307 (2019)]. In this paper, we present algorithms for the two formulations of electromagnetic gyrokinetics: the Hamiltonian and the symplectic. We describe each formulation and show results of benchmark tests. Although our scheme works for the Hamiltonian formulation, the presence of spurious numerical modes for high-β and large $$k^2_⊥ρ^2_s$$ regimes shows that the symplectic formulation is more robust. We then review our recent algorithm for the symplectic formulation [Mandell et al., J. Plasma Phys. 86, 905860109 (2020)], along with example application of this new capability. Maintaining positivity of the distribution function can be challenging, and we describe a new and novel exponential recovery based algorithm to address this.
Research Organization:
Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
Sponsoring Organization:
USDOE
Grant/Contract Number:
AC02-09CH11466; FC02-08ER54966; FG02-04ER54742; FG02-95ER54309; FG02-97ER25308
OSTI ID:
1668050
Alternate ID(s):
OSTI ID: 1608504
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 (44)

Electron Physics in 3-D Two-Fluid 10-Moment Modeling of Ganymede's Magnetosphere journal April 2018
Nonlinear Dynamics in the Tokamak Edge journal May 2010
The Runge–Kutta Discontinuous Galerkin Method for Conservation Laws V journal April 1998
A High-Order Discontinuous Galerkin Method for 2D Incompressible Flows journal May 2000
The Serendipity Family of Finite Elements journal March 2011
Comparison of initial value and eigenvalue codes for kinetic toroidal plasma instabilities journal August 1995
Computational design for long-term numerical integration of the equations of fluid motion: Two-dimensional incompressible flow. Part I journal August 1966
A δf particle method for gyrokinetic simulations with kinetic electrons and electromagnetic perturbations journal August 2003
A numerical instability in an ADI algorithm for gyrokinetics journal November 2005
A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm journal December 2015
A high resolution wave propagation scheme for ideal Two-Fluid plasma equations journal November 2006
Electromagnetic gyrokinetic PIC simulation with an adjustable control variates method journal July 2007
Discontinuous Galerkin algorithms for fully kinetic plasmas journal January 2018
Plasma-material boundary conditions for discontinuous Galerkin continuum-kinetic simulations, with a focus on secondary electron emission journal April 2020
Gyrokinetic continuum simulation of turbulence in a straight open-field-line plasma journal May 2017
Effect of a weak ion collisionality on the dynamics of kinetic electrostatic shocks journal February 2019
Electromagnetic full- gyrokinetics in the tokamak edge with discontinuous Galerkin methods journal February 2020
Runge-Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems journal September 2001
Global Ten‐Moment Multifluid Simulations of the Solar Wind Interaction with Mercury: From the Planetary Conducting Core to the Dynamic Magnetosphere journal November 2019
An Extended MHD Study of the 16 October 2015 MMS Diffusion Region Crossing journal November 2019
Gyrokinetic turbulence simulations with kinetic electrons journal May 2001
Conventional δf-particle simulations of electromagnetic perturbations with finite elements journal December 2004
Model Fokker-Planck Equation for a Plasma and Its Solution journal January 1964
Intermittent convective transport carried by propagating electromagnetic filamentary structures in nonuniformly magnetized plasma journal February 2010
Finite-β simulation of microinstabilities journal February 2014
Pullback transformation in gyrokinetic electromagnetic simulations journal September 2014
Comparison of multi-fluid moment models with particle-in-cell simulations of collisionless magnetic reconnection journal January 2015
A gyrokinetic one-dimensional scrape-off layer model of an edge-localized mode heat pulse journal February 2015
Continuum kinetic and multi-fluid simulations of classical sheaths journal February 2017
Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code journal May 2017
Simulations of anti-parallel reconnection using a nonlocal heat flux closure journal August 2017
Nonlinear saturation of the Weibel instability journal November 2017
A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons journal February 2018
Full- f gyrokinetic simulation of turbulence in a helical open-field-line plasma journal January 2019
Gyrokinetic continuum simulations of plasma turbulence in the Texas Helimak journal April 2019
Collisional effects on kinetic electromagnetic modes and associated quasilinear transport journal January 1987
Gyrokinetic field theory journal February 2000
Three-dimensional computation of drift Alfvén turbulence journal October 1997
Low Mach-number collisionless electrostatic shocks and associated ion acceleration journal January 2018
Maximum-principle-satisfying and positivity-preserving high-order schemes for conservation laws: survey and new developments journal May 2011
Plasma Oscillations with Diffusion in Velocity Space journal December 1958
Foundations of nonlinear gyrokinetic theory journal April 2007
Hamiltonian theory of guiding-center motion journal May 2009
Temperature-dependent Saturation of Weibel-type Instabilities in Counter-streaming Plasmas journal February 2019

Similar Records

Fluid and gyrokinetic turbulence in open field-line, helical plasmas
Journal Article · Sun Aug 02 20:00:00 EDT 2020 · Physics of Plasmas · OSTI ID:1647927

Electromagnetic total-f algorithm for gyrokinetic particle-in-cell simulations of boundary plasma in XGC
Journal Article · Sun Nov 20 19:00:00 EST 2022 · Physics of Plasmas · OSTI ID:1905772

Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems
Journal Article · Sun Nov 20 19:00:00 EST 2016 · Physics of Plasmas · OSTI ID:1535287